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SPACECRAFT  AND  STELLAR  OCCULTATIONS  BY  TURBULENT  PLANETARY 
ATMOSPHERES 

A  theoretical  investigation  of  various  wave  propagation  effects  and  their  impact  on 
derived  profiles  of  refractivity ,  temperature  and  pressure 


ABSTRACT 

The  long  propagation  paths  involved  in  radio  and  stellar  occultations  by  turbulent 
planetary  atmospheres  require  that  the  classical,  weak  scattering  scintillation  theory 
be  expanded  to  account  for  the  inhomogeneous  ambient  atmosphere  upon  which  the 
turbulence  is  superimposed.  Such  coupling  between  the  turbulent  and  the  ambient 
components  of  refractivity  reduces  the  scintillation  index  by  less  than  a  factor  of  two 
in  shallow  radio  occultations.  For  stellar  occultations  the  reduction  varies  between 
this  value  for  very  small  projected  stellar  radii,  to  approximately  a  factor  of  3.6  when 
the  projected  stellar  radius  above  the  planetary  limb  greatly  exceeds  the  radius  of  the 
free -space  Fresnel  zone  at  this  distance.  More  profound  changes  are  found  in  the 
scintillation  power  spectrum,  the  shape  of  which  depends  strongly  on  both  occulta- 
tion  depth  and  geometry  when  coupling  to  the  inhomogeneous  background  is  pro¬ 
perly  accounted  for. 

Second-order,  systematic  propagation  effects  calculated  from  both  geometrical  optics 
and  a  wave-optical  formulation  show  that  the  average  phase  velocity  is  increased  in 
the  presence  of  turbulence.  The  finite  (X’l/6)  wavelength  dependence  of  the  phase 
path  bias  implies  that  an  initially  non-dispersive  medium  becomes  slightly  dispersive 
by  the  addition  of  turbulence. 

Atmospheric  profiles  derived  from  occultations  by  turbulent  planetary  atmospheres 
differ  only  slightly  from  those  of  the  corresponding  non-turbulent  atmosphere  when 
the  weak  scattering  condition  is  satisfied.  Even  so,  there  is  strong  indication  that 
profiles  obtained  from  radio  Doppler  measurements  are  much  less  affected  by  turbu 
lence  than  those  from  radio  or  stellar  intensity  measurements.  Systematic  errors  in 
refractivity.  temperature  and  pressure  arc  all  small  for  either  kind  of  measurement  in 
the  limit  of  weak  scattering. 

Radio  and  stellar  occultations,  when  properly  conducted  and  interpreted,  are  a  poten¬ 
tial  source  of  information  also  on  turbulence  strength,  bulk  flow  scales,  planetary 
rotation  and  atmospheric  winds. 


1  INTRODUCTION 
1.1  Historical  perspective 

I'hcorctical  interpretation  of  the  occultation  of  stars  by  planetary  atmospheres 
appears  to  have  originated  with  Pannekoek  (1903),  who  observed  the  occultation  of  a 
star  by  Jupiter  on  September  19  of  that  year.  He  recommended  that  such  observa¬ 
tions  be  made  carefully:  "Because  they  can  give  a  determination  of  the  temperature 
distribution  in  the  outermost  layers  of  the  planetary  atmosphere,  or  more  directly, 
the  way  in  which  the  horizontal  refraction  varies  with  height”.  The  theoretical  prob¬ 
lem  of  deriving  a  refractivity  profile  from  a  luminosity  record  of  the  occulted  star 
was  later  reviewed  by  Fabry  (1929),  by  Baum  and  Code  (1953),  and  finally  by  Link 
(1969),  who  was  the  first  to  employ  numerical  inversion  techniques  to  the  optical 
data. 

To  date  observations  of  stellar  occultations  have  yielded  information  on  the  atmo¬ 
sphere  of  Jupiter  through  the  occultation  of  a  Arictis  (Baum  and  Code  1953)  and  the 
(J-Scorpii  system  (eg  Hubbard  cl  at  1972,  Vevcrka  ct  at  19741)).  For  a  review  of  stellar 
occultations  by  Jupiter,  sec  Hunten  and  Vevcrka  (1976).  Structural  information  on 
the  upper  atmosphere  of  Neptune  was  similarly  obtained  by  the  occultation  of 
BI)-17°4388  on  7  April  1968  (Osawa  dal  1968,  Kovalevsky  and  Link  1968,  Freeman 
and  Lynga  1970,  Rages  ct  at  1 974,  Vevcrka  d  at  1 974a).  Stellar  occultations  by  Mars 
(Tcxas-Arizona  Occultation  Croup  1977,  Elliot  d  at  1977b,  Young  1977),  by  Venus 
(dc  Vaucoulcurs  and  Mcn/.c!  1960,  Vevcrka  and  Wasserman  1974)  and  by  Uranus 
(Elliot  d  at  1977a,  Millis  d  at  1977,  Hubbard  d  at  1977)  have  also  been  recorded. 


With  the  advent  of  spacecraft  missions  to  the  planets,  occultation  of  the  spacecraft 
itself  has  become  a  very  important  and  complementary  source  of  information  of  the 
structure  of  planetarv  atmospheres.  While  stellar  occulta! ions  are  limited  to  fairly  tenu¬ 
ous  atmospheric  regions,  spacccrait  occultations,  by  virtue  ol  their  much  higher  signal-to- 
noisc  ratio,  smaller  planet-to-spacccraft  distance,  and  the  Doppler  frequency  as  a  new 
data  source,  have  opened  up  for  study  much  deeper  atmospheric  regions  not  acces¬ 
sible  to  earth-based  optical  observations.  1'his  subject  was  one  of  several  related  bista¬ 
tic  radar  experiments  introduced  by  V  R  Fshleman  at  the  1962  Review  of  Space 
Research  conducted  at  the  State  University  of  Iowa  by  the  US  National  Academy  of 
Sciences  (Publication  1079),  and  carle  reports  be  Fjcldbo  (1904),  Fshleman  (1964), 
and  Kliore  ct  al  ( 1  964)  describe  details  of  its  scientific  and  engineering  aspects.  Since  the 
first  encounter  of  Mars  by  Mariner  4  on  If)  July  1965  (Kliore  ct  al  1 9 6 f> ,  Fjcldbo  ct 
al  1966a,  Fjcldbo  and  F.shleman  1968),  the  subsequent  occultations  of  Mariner  a  and 
10  be  Venus  (e  g  Kliore  ct  al  1967,  Mariner  Stanford  Group  1967.  Fjcldbo  and 
F.shleman  1969,  Fjcldbo  ct  al  1971,  Howard  ct  al  1974,  Fjcldbo  ct  al  197")),  of  Pio¬ 
neer  10  and  11  by  Jupiter  te  g  Kliore  ct  al  1974,  Kliore  ct  al  1975,  Kliore  and 
Woiceshyn  1976,  Fjcldbo  ct  al  i 975).  of  Pioneer  1  1  at  Saturn  (Kliore  ct  al  1980)  and 
the  occultations  of  Yovagcr  1  and  2  bv  Jupiter  ( F.shleman  <7  al  1979b)  have  yielded 
important  new  results  on  the  structure  of  both  the  atmosphere  and  ionosphere  of 
these  very  planets.  Additional  in  situ  information  on  temperature,  pressure  and  chemi¬ 
cal  composition  have  been  provided  bv  the  Venera  landers  descending  into  the  atmo¬ 
sphere  of  Venus  (e  g  X'akhnin  1968,  Vinogradov  ct  al  1968)  and  recently  also  by  the 
Pioneer  Venus  entry  probes  ( Seif f  ct  al  1979). 

As  the  potential  viability  of  the  occultation  tec  hnique  has  been  demonstrated  in  a 
number  of  both  optical  and  radio  occultations,  a  question  of  utmost  importance 
concerns  the  accuracy  of  the  information  on  atmospheric  structure  obtained  in  this 
way.  As  it  turns  out,  the  answer  to  this  question  is  far  from  simple:  sources  of  widely 
different  character  and  magnitude  contribute  to  the  ultimate  error.  It  is  remarkable 
that  not  until  the  early  analyses  of  the  Pioneer  10  and  II  occultations  by  Jupiter 
produced  grossly  unrealistic  results,  was  the  question  of  accuracy  raised  on  a  serious 
and  systematic  basis  leg  F.shleman  1975,  Hubbard  ct  al  1975). 

Prior  to  these  difficulties,  to  which  1  shall  return  shortly,  stellar  occultation  profiles 
ad  long  been  suspected  ol  containing  substantial  errors.  An  outstanding  and  contri¬ 
buting  fac  tor  has  here  been  the  large-  sfn'kcs  in  intensity  that  appear  to  be  a  regular 
feature  of  stellar  occultation  light  curves.  Without  doubt,  these  features,  caused  hv 
some  as  yet  unknown  vertical  microstructure  in  the  atmosphere,  have  hampered  the 
derivation  ol  a  refractivity  profile  from  the  occultation  data,  in  particular  in  the 
earlier  clays  before  the  numerical  inversion  technique  was  introduced  bv  Kovalevsky 
and  Link  (1968)  to  fac  ilitate  the  analysis  of  the  Bl)-17°4888  occultation  by  Neptune. 
Hut  it  is  also  clear  that  other  major  sources  of  error  are  present.  Perhaps  the  most 
dramatic  example  ol  the  difficulties  associated  with  this  technique  is  provided  by  the 
discrepant  refractivity  profiles  obtained  by  several  research  groups  for  the  (3-Scorpii 
occultation  by  Jupiter.  Ol  the  six  recorded  profiles  only  three  seem  to  approach  the 
factual  temperature  distribution  in  the  upper  Jovian  atmosphere  (Vapilion  dal  1978, 
Hubbard  ct  al  19/  2 .  Yevcrka  ct  al  1974b).  For  a  review  of  the  (5-Scorpii  occultations, 
see  Huntcn  and  Yevcrka  (1976).  It  is  clear  that  meaningful  results  demand  extreme 
•  are  in  both  the-  experimental  design  phase  and  also  during  the  recording  of  data.  It  is 
also  essential  that  due  regard  be  paid  to  the  different  potential  sourc  es  of  error  that 
may  arise  in  the  subsequent  data  reduction  to  obtain  profiles  of  refractivity.  tempera¬ 
ture  and  pressure  (sec  eg  Wasserman  and  Vcverka  1978,  Frenc  h  ct  al  1978). 

When  the  first  radio  oc  c  ultation  profiles  from  the-  Mariner  4  Mars  encounter  appeared 
(Kliore  dal  196),  l-jeldbo  dal  1966a, b,  Fjcldbo  and  F.shleman  1968),  there  was  no 
reason  to  doubt  their  general  validity,  except  that  the  upper  lew  scale  heights  of 
these  profiles  were  somewhat  indeterminate  because  of  partially  unknown  boundary 
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conditions,  a  problem  common  to  both  radio  and  optical  occultations.  The  sub¬ 
sequent  occultations  of  Mariner  5  and  10  by  Venus  on  19  October  1967  and  5  feb- 
ruary  1974,  respectively,  also  yielded  results  that  were  in  all  major  respects  reasonable 
(Kliore  et  al  1967,  Mariner  Stanford  Group  1967,  Howard  et  al  1974).  Indeed,  com¬ 
parisons  between  the  Mariner  5  profiles  and  that  of  the  Venera  4  landing  module 
(Eshleman  et  al  1968,  Eshleman  1970),  and  between  the  Mariner  5  and  10  profiles 
and  the  Venera  8,  9  and  10  data  (Nicholson  and  Muhleman  1978),  display  a  substan¬ 
tial  degree  of  consistence,  both  between  the  individual  radio  occultation  profiles  and 
between  these  profiles  and  those  of  the  landers  where  a  comparison  could  be  made. 
On  a  finer  level  it  is  nevertheless  clear  that  the  small  differences  that  still  do  exist 
between  the  different  Mariner  10  profiles  are  not  due  to  either  statistical  fluctuations 
in  the  data  processing,  or  in  differing  assumptions  regarding  the  initial  conditions  for 
the  profile  computations  (Lipa  and  Tyler  1979).  Differences  by  as  much  as  10  K 
between  the  profile  computations  of  Nicholson  and  Muhleman  (1978)  and  Lipa  and 
Tyler  (1979)  and  the  initial  profile  of  Howard  et  al  (1974)  probably  require  an 
explanation  in  terms  of  the  difference  in  computational  procedures  used  to  obtain 
these  profiles  (Lipa  and  Tyler  1979). 

1  he  high  confidence  in  the  radio  occultation  technique,  as  a  tool  for  probing  plane¬ 
tary  ionospheres  and  atmospheres  to  depths  not  accessible  to  optical  occultations,  was 
seriously  shaken  bv  the  publication  of  the  preliminary  temperature  -pressure  profiles 
obtained  from  the  Pioneer  10  and  11  encounters  with  Jupiter  on  3  December  1973 
and  4  December  1974,  respectively  (Kliore  ct  al  1974,  Kliore  ct  al  1975).  Analysis  of 
the  ionospheric  part  was  hampered  by  multipath  propagation  until  the  several  simul¬ 
taneous  signals  were  laboriously  sorted  out  by  hand  ( Ljeldbo  et  al  1975).  In  the  lower 
neutral  atmosphere,  temperatures  were  far  in  excess  of  those  obtained  from  radiative- 
convective  balance  calculations  and  on-board  radiometric  observations  (Trafton  1973, 
Wallace  ct  al  1974,  Orton  and  Ingersoli  1976).  The  radiometric  observations  fixed  the 
effective  atmospheric  temperature  at  125±2K  at  the  500  mb  level,  far  below  the 
5  400  K  computed  from  the  first  analyses  of  the  Pioneer  1  0  and  1  1  occultation  data, 
this  circumstance,  further  aggrevated  by  Gulkis’  (Bcrge  and  Gulkis  1976)  argument 
that  the  13  cm  S-band  radio  signal  could  not  propagate  below  the  280  K  level  corres¬ 
ponding  to  the  13  cm  atmospheric  brightness  temperature  of  Jupiter,  led  one  to 
conclude  that  the  Pioneer  10  and  11  radio  occultation  profiles  were  in  all  probability 
wrong. 

This  unexpected  situation  spurred  hectic  efforts  to  uncover  the  problem  with  the 
Pioneer  occultation  profiles.  A  major  source  of  error,  the  neglect  of  planetary  oblate¬ 
ness  in  the  inversion  calculation  to  obtain  a  refractivity  profile  from  the  atmospheric- 
phase  perturbations,  was  identified  by  Hubbard  ct  al  (1975);  sec  also  Hubbard  (1976). 
Almost  simultaneously  Eshleman  (1975)  performed  a  general  parametric  sensitivity 
analysis  of  radio  occultation  experiments,  showing  that  small  errors  in  cither  Doppler 
rate,  spacecraft-limb  distance  or  in  spacecraft  velocity,  were  magnified  by  a  poten¬ 
tially  large  factor  in  such  experiments.  In  the  Pioneer  10  entry  occulation  this  factor 
became  as  large  as  =70.  Eshleman  also  showed  that  there  was  more  than  one  source 
of  major  error  in  the  original  profiles,  although  correction  for  the  oblatcncss  or 
velocity  factor  did  make  it  possible  to  produce  several  temperature-pressure  profiles 
that  are  both  mutually  consistent  and  in  reasonable  agreement  with  the  Pioneer 
radiometer  observations,  in  particular  in  deeper  regions  of  the  atmosphere,  below 
about  10  mb,  where  the  effect  of  partially  unknown  boundary  conditions  is  rapidly 
decreasing  (Kliore  ct  al  1976,  Kliore  and  W'oiccshyn  1976). 

Simultaneously,  and  independently  of  the  events  described  above,  a  different  line  of 
development,  from  which  also  the  present  work  arose,  sought  the  explanation  for  the 
erroneous  Pioneer  profiles  in  terms  of  smaller  scale  refractive  anomalies,  such  as 
turbulence,  in  the  Jovian  atmosphere.  While  layers,  vertically  propagating  internal 
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gravity  waves  or  turbulence  had  long  been  suspected  of  causing  the  large  spikes  seen 
in  stellar  occultation  light  curves,  virtually  no  systematic  effort  was  expended  up  to 
the  year  of  1976  to  explore  these  possibilities.  For  radio  occultation  a  first  quantita¬ 
tive  step  was  taken  by  Hubbard  and  Jokipii  (1975),  who  calculated  the  second-order, 
systematic  effect  of  turbulence  on  radio  occultation  profiles  derived  from  Doppler 
measurements.  Hubbard  and  Jokipii  showed  that  if  exponential  turbulence  was 
imbedded  in  an  otherwise  isothermal  atmosphere,  the  average  radio  phase  would  be 
retarded,  progressively  with  increasing  atmospheric  depth,  and  yield  a  resulting  tem¬ 
perature  profile  that  is  slightly  offset  from  the  true  profile  of  the  quiescent  atmo¬ 
sphere.  The  Hubbard-Jokipii  calculation  was  later  critisized  by  Eshleman  and  Haug- 
stad  (1977,  1978)  and  Haugstad  (1976),  who  showed  that  the  average  phase  bias 
caused  by  the  turbulence  was  in  the  opposite  direction  (i  c,  the  average  phase  was 
advanced),  and  furthermore  that  the  analysis,  contrary  to  their  claim,  was  valid  only 
in  the  limit  of  geometrical  optics  which  is  never  realized  in  radio  occultation  experi¬ 
ments.  It  was  also  demonstrated  that  taking  proper  account  of  diffraction  effects 
both  yielded  a  wavelength  dependence  of  the  phase  path  bias  and  reduced  its  magni¬ 
tude  far  below  the  value  calculated  from  geometrical  optics  (Eshleman  and  Haugstad 
1978,  Haugstad  1978b, c,d,  Haugstad  and  Eshleman  1979).  While  the  precise  form  of  this 
systematic  effect  of  turbulence,  its  magnitude  and  the  question  of  wavelength  depen¬ 
dence,  is  still  under  debate,  both  parties  appear  to  agree  that  this  propagation  effect 
is  too  small  to  be  of  any  practical  concern  (Hubbard  and  Jokipii  1977a,b,  Haugstad 
and  Eshleman  1979,  Hubbard  and  Jokipii  1979,  Hubbard  i979,  Eshleman  and  Haug¬ 
stad  1979). 

Another,  and  in  retrospect  more  important  step  to  clarify  the  effect  of  turbulence, 
was  taken  by  A  T  Young  (1976).  Young  considered  the  larger,  first-order  random 
effects  of  turbulence  and  argued  forcefully  that  the  intensity  spikes  typical  of  stellar 
occultation  luminosity  curves  were  uncorrelated  on  a  planetary  scale  and  were  in  all 
probability  manifestations  of  turbulence  rather  than  layers  or  other  horizontally 
extended  structures.  Developing  a  hybrid  geometrical  optics  —  wave  optics  scintilla¬ 
tion  theory  that  accounted  for  "coupling”  between  the  turbulence  and  the  inhomo¬ 
geneous  background  upon  which  it  is  superimposed.  Young  argued  that  also  radio 
occultation  experiments  should  be  severely  contaminated  by  turbulence  effects. 
Although  Young’s  heuristic  approach  did  not  yield  the  correct  form  of  this  coupling 
effect,  which  appears  as  a  result  of  propagation  over  very  long  distances,  his  paper  did 
bring  into  focus  turbulence  as  a  major  concern  and  undoubtedly  stimulated  other, 
independent  attempts  to  develop  the  correct  theoretical  basis  for  the  first-order, 
random  effects  of  turbulence  (Ishimaru  1977,  Haugstad  1977,  Haugstad  1978b, c,d,e, 
Hubbard  ct  al  1978,  Haugstad  1979a,  Woo  et  al  1980).  These  analyses  showed  that 
Young’s  heuristic  model  seriously  overestimated  the  effect  of  turbulence  and  further¬ 
more  indicated  that  previous  spacecraft  occultation  profiles  may  not  have  been 
seriously  affected  by  turbulence  (Haugstad  1978d,e). 

Along  a  parallel  line  of  development,  interest  was  shifted  from  the  possible 
degrading  effect  of  turbulence  on  radio  occultation  profiles,  to  the  possibility  of 
extracting  from  the  radio  signal  characteristics  of  the  underlying  turbulence.  As 
stressed  by  Woo  et  al  (1974).  information  on  the  strength  and  spatial  extent  of 
turbulence  is  important  as  a  means  of  constraining  theoretical  circulation  models  of 
planetary  atmospheres.  Not  surprisingly,  developing  procedures  for  extracting  such 
information  turned  out  to  be  but  another  aspect  of  the  general  problem  of  electro¬ 
magnetic  wave  propagation  in  a  randomly  inhomogeneous  medium  with  inhomo¬ 
geneous  (and  possibly  anisotropic)  mean  characteristics.  The  pioneering  work  here  was 
done  by  R  Woo  and  A  Ishimaru,  who  applied  standard  weak  scintillation  theory  to 
conditions  typical  of  radio  occultations  of  planetary  spacecraft  (Woo  and  Ishimaru 
1973,  Woo  and  Ishimaru  1974).  The  theory,  not  yet  accounting  for  the  coupling 
effect  discussed  earlier,  was  subsequently  applied  to  the  Mariner  5  and  10  occultations 
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to  yield  rough  estimates  of  the  scintillation  index  and  structure  constant  in  regions 
where  turbulence  could  be  reliablv  identified  (Woo  et  al  1974,  Woo  1975a);  see  also 
Gurvich  ( 1969),  and  Golitsyn  and  Gurvich  ( 1 97 1 ).  The  theory  was  later  expanded  by 
Haugstad  ( 1978c, d,e),  Hubbard  et  al  ( 1978),  and  by  Woo  et  al  (1980)  to  include  effects 
of  coupling  to  the  inhomogeneous  ambient  atmosphere.  Power  spectra  of  phase 
(Haugstad  1979a)  and  intensity  scintillations  (Haugstad  1979a,  Woo  et  al  1980)  were 
calculated  for  a  range  of  different  occultation  geometries  and  discussed  in  terms  of 
their  potential  usefulness  for  deriving  turbulence  characteristics  from  the  radio  occul¬ 
tation  data. 

Even  though  the  past  5—6  years  have  seen  major  advances  in  the  theory  of  turbulence 
effects  in  radio  and  optical  occultation*,  fundamental  problems  still  remain  unsolved.  A 
common  feature  of  all  the  turbulence  theories  quoted  here,  for  example,  is  their 
limitation  to  weak  scattering  conditions.  While  strong  scattering  theories  suitable  for 
short-path,  terrestrial  propagation  conditions  have  been  successfully  formulated  during 
the  past  decade  (for  a  review,  see  eg  Fante  1975),  inclusion  of  the  inhomogeneous 
background  essential  to  radio  and  optical  occultation  conditions,  still  poses  an 
unsolved  problem  in  the  context  of  strong  scattering.  This  and  other  related  prob¬ 
lems,  to  be  considered  in  the  next  section,  undoubtedly  constitute  a  principal  chal¬ 
lenge  to  future  theoretical  work  in  this  field. 


1.2  Objective  and  outline 

The  general  problem  to  be  addressed  in  this  report  is  twofold: 

a)  To  establish,  by  theoretical  analysis,  the  effect  of  turbulence  on  radio  and  stellar 
occultation  temperature  and  pressure  profiles 

b)  To  provide  a  theoretical  basis  from  which  the  characteristics  of  turbulence  in 
planetary  atmospheres  can  be  denied  from  occultation  measurements 


In  order  to  uncover  the  various  facets  of  this  extensive  problem  and,  at  the  same 
time,  provide  a  justification  for  the  particular  aspects  of  this  problem  that  will  be 
subjected  to  detailed  examination  here,  it  is  convenient  to  consult  the  schematic 
illustration  in  Figure  1.1.  In  order  to  solve  the  problem  specilied  above,  an  essential 


first  step  will  be  to  formulate  a  realistic 


Figure  1.1  Schematic  illustration  of  the  effect 
of  turbulence  on  atmospheric  pro¬ 
files  derived  from  occultation 
measurements 


model  of  both  the  quiescent,  ambient 
atmosphere  and  the  turbulence. 
Once  established,  the  problem  con¬ 
sists  of  calculating  its  effect  on  the 
probing  electromagnetic  signal.  The 
resulting  signal,  contaminated  by 
effects  due  to  turbulence,  is  then 
subjected  to  the  standard  inversion 
algorithm  (i  c,  to  its  mathematical 
representation),  by  which  a  profile 
of  rcfractivity  versus  height  is  con¬ 
structed  from  the  probing  signal. 
Knowing  the  effect  of  the  turbulence 
on  this  signal  allows  one,  at  least  in 
principle,  to  calculate  also  the  corres¬ 
ponding  effect  of  the  turbulence  on 
the  rcfractivity  profile.  If,  further¬ 
more.  the  chemical  composition  of 
the  atmosphere  is  known,  also  pro¬ 
files  of  temperature  and  pressure 
may  be  computed  and  separated,  at 
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least  in  principle  and  in  an  average  sense,  into  an  ambient  part  and  a  fictitious 
component  caused  by  turbulence. 

While  this  procedure  is  in  principle  complete  and  will  solve  the  stated  problem, 
difficulties  arise  on  different  levels  when  practical  solutions  are  sought.  For  instance, 
when  applying  the  result  of  an  analysis  of  this  sort  to  am  specific  oceultat ion,  a 
principal  difficulty  arises  from  the  fact  that  the  characteristics  of  the  turbulence, 
which  will  affect  the  outcome  of  the  wave  propagation  analvsis,  is  usually  not  known 
u  priori,  l  he  problem  is  therefore  iterative  in  nature:  from  a  crude  model  of  the 
turbulence  its  effect  on  the  probing  signal  may  be  computed  and  compared  with  the 
actually  observed  signal.  From  this  comparison  a  refined  turbulence  model,  closer  to 
the  factual  distribution  of  turbulence,  may  in  principle  be  constructed,  etc. 

but  also  major  technical  problems  arise  when  implementation  of  the  scheme  in  Figure 
1.1  is  attempted.  As  in  terrestrial  propagation  experiments,  the  effect  of  turbulence 
or.  electromagnetic  radiation  frequently  falls  into  one  of  the  two  limiting  categories 
of  weak  or  strong  scattering.  In  the  former  case  the  fluctuations  in  signal  strength  arc 
small  in  terms  of  the  mean  signal,  whereas  in  strong  scattering  the  root-mean-square 
lluctuations  are  comparable  to  the  mean  signal  strength  (saturated  scintillations). 
While  weak  scattering  theories  were  fully  developed  by  the  mid-sixties,  only  more 
recently  have  strong  scattering  theories,  due  to  their  greater  mathematical  complexity, 
been  developed  into  explicit  forms  that  allow  confrontation  with  observations.  In 
radio  occupations  two  additional  factors  add  to  the  complexity  of  this  problem: 
ti  the  actual  inhomogeneous  nature  of  the  ambient  atmosphere  cannot  be  neglected, 
a-  is  routinely  done  in  propagation  theories  tailored  to  terrestrial  conditions  involving 
much  shorter  path  lengths,  and  (ii)  the  fluctuations  in  angle  of  arrival  of  the  wave- 
front,  a  key  quantity  for  assessing  the  effect  of  turbulence  on  atmospheric  profiles 
derived  from  Doppler  measurements,  requires  also  the  precise  spatial  extent  of  the 
turbulence  to  be  stipulated  (local  properties  are  not  sufficient)  and  leads  to  a  very 
•  omplicated  analysis.  (For  stellar  occupations  only  the  former  of  these  complications 
applies. I  We  also  note  that,  because  of  the  non-linear  relationship  between  the  turbu¬ 
lent  refructivity  fluctuations  and  the  amplitude  and  phase,  the  calculated  profiles  of 
temperature  and  pressure  contain  both  a  systematic  and  a  random  part.  During  weak 
-  titering  conditions  the  systematic  effect  will  be  of  second  order  in  small  quantities, 
and  hence  much  smaller  than  the  leading  first-order,  random  effect. 

Having  described  and  categorized  above  the  full  problem,  we  here  restrict  ourselves  to 
giving  a  detailed  treatment  of  certain,  however  important,  aspects  of  the  overall 
problem,  essentially  covering  and  expanding  on  the  material  developed  in  the  quoted 
publications  by  this  author  (Ilaugstad  1978b, c,d,e,  Ilaugstad  1979a,  Kshleman  and 
ilaugstad  1977.  Kshleman  and  Ilaugstad  1978,  Kshleman  and  Ilaugstad  1979.  Haug¬ 
stad  and  Kshleman  1979,  Ilaugstad  1980).  I  bus,  giving  an  expose  of  the  basic  theory 
of  radio  and  stellar  occultations  in  section  2.  we  proceed  to  develop  in  section  3  the 
theory  of  electromagnetic  wave  propagation  in  a  non-turbulcnt  atmosphere  as  the 
limiting  form  of  a  general  wave-optical  formulation.  From  this  approach  the  error 
incurred  by  the  routine  neglect  of  wave -optical  effects  in  the  conventional  occultation 
theories  is  readily  assessed.  In  section  4  the  iccak  scattering  theory  of  radiation  propa¬ 
gating  through  a  turbulent  planetary  atmosphere  is  developed  from  a  stationary  phase 
method  (Ilaugstad  1978c).  After  having  established  an  atmospheric  model,  the 
resulting  formulae  arc  used  to  derive  explicit  results  for  the  weak  scattering  scintilla¬ 
tion  index  in  both  radio  and  stellar  occultation  experiments.  We  establish  here  also 
the  form  of  the  power  spectra  of  the  fluctuations  in  phase  and  intensity,  and  com¬ 
pare  tfu:  predicted  spectra  with  experimental  data  from  the  encounters  of  Mariner  5 
and  10  with  Venus.  Section  4  concludes  with  a  general  discussion  of  the  validity  of 
the  weak  scattering  approximation,  which  is  a  fundamental  limitation  of  all  the  re¬ 
sults  derived  in  this  work.  In  section  we  investigate  to  what  extent  turbulence,  again 
within  the  realm  of  weak  scattering,  affects  the  derived  profiles  of  temperature  and 
pressure.  Answering  this  question  turns  out  to  be  relatively  simple  for  atmospheric 
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profiles  derived  from  stellar  or  radio  occultation  intensity  measurements.  However, 
for  atmospheric  profiles  derived  from  radio  Doppler  measurements  the  technical  prob¬ 
lems  are  formidable,  as  previously  indicated.  For  this  case  a  physical,  semi-quantita¬ 
tive  approach  avoiding  the  core  of  the  mathematical  difficulties  has  been  adopted. 
The  section  concludes  with  a  discussion  of  the  small  systematic  effects  of  turbulence 
on  occultation  profiles.  In  section  6  we  devote  attention  to  a  relatively  new  topic  in 
this  field:  the  effect  of  turbulence  in  deep  radio  occultation  measurements.  The  basic- 
theory  of  deep  radio  occultations  (Eshleman  et  al  1979a)  is  combined  here  with  the 
scintillation  theory  developed  in  section  4  to  yield  results  for  the  scintillation  index 
ami  p'>u<i  spectrum  close  to  a  refractive  focus  (Haugstad  1981).  The  theory  is 
applied  to  deep  radio  occultations  o|  planetars  atmospheres  (the  Voyager  1  Jupiter 
encounter),  and  to  t he  suggested  use  of  the  gravitational  lens  of  the  sun  for  eaves¬ 
dropping  and  communication  over  interstellar  distances  (Eshleman  1979).  We  con¬ 
clude  the  present  work  in  section  7  by  discussing  a  general  procedure  by  which 
information  on  the  strength  and  spatial  distribution  of  turbulence  may  be  extracted 
from  occultation  measurements,  exploring  both  the  potential  feasibility  and  the  inhe¬ 
rent  limitations  of  this  scheme.  Section  8  summarizes  the  results  of  the  present  work, 
and  expands  the  perspective  by  recommending  where  future  work  in  this  field  should 
be  directed. 


2  RADIO  AND  STELLAR  OCCULTATIONS;  FORMAL  THEORY 
2.1  Profiles  of  refractivity,  temperature  and  pressure 

1  he  derivation  of  vertical  profiles  of  temperature  and  pressure  from  occultation 
measurements  proceeds  in  two  stages.  In  the  first  a  refractivity  profile  is  determiner! 
from  the  atmospheric  perturbations  of  cither  phase  or  intensity  of  the  probing  signal. 
In  the  second  stage  the  refractivity  profile  is  used  together  with  possible  knowledge 
about  the  chemical  composition  of  the  atmosphere  to  calculate  profiles  of  temperature 
and  pressure. 

For  the  first  of  these  steps,  consider  Figure  2.1  which  illustrates  a  spacecraft  partially 
occulted  by  a  planetars'  atmosphere,  the  refractivity  of  which  is  assumed  to  be  a 


ligure  2  I  Hay  path  and  simplified  occultation  geometry'  defining  the  rax  I tending 
angle  a,  the  impact  parameter  a,  the  radial  distance  to  the  point  of  closest 
approach  (periapsis)  r0,  and  the  planet  to  spacecraft  distance  I)  as  measured 
from  the  planetary  centre 

The  spacecraft  velocity  is  resolved  into  a  component  parallel  to  the  incoming  rav 
asymptote,  and  a  component  vt  in  the  plane  of  the  sky. 

I 

j 
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function  of  the  radial  coordinate  only-  (For  stellar  occultations  the  earth  replaces  the 
spacecraft  and  1)  is  now  the  planet-to-earth  distance.)  It  is  assumed  throughout  that 
ro  >  r0c,  where  roc  is  the  level  of  superrefraction,  below  which  the  ray  encircles  the 
planet  and  fails  to  reach  the  spacecraft.  (For  a  discussion  of  this  interesting  possibi¬ 
lity,  see  Lipa  and  Croft  (1975),  who  also  calculate  the  refractive  attenuation  of  the 
critical  ray.) 

Let  the  refractive  index  n  =  1+-  v,  where  v  is  the  atmospheric  refractivity.  The  quanti¬ 
ties  n  and  a  are  related  by  an  Abelian  transform  pair  (eg,  Fjeldbo  et  al  1971, 
Phinncy  and  Anderson  1968) 


«(»)  -  —2a  7  f  tAtts 

a  d£  n^  -  a^ ) 

(2.1a) 

tr  In  n  (r0)  =  —  /  ^7^  cosh'1  (T)da' 
a  da  a 

(2.1b) 

where  £  =  nr,  and  where  r0  and  a  are  related  by  Bouguer's  rule  (Born  and  Wolf  1975) 
as  in 

a  =  n(r0)r0  (2.2) 

Let  k  (k=  1  k  I  =  2tr/X)  be  the  vector  wavenumber  of  the  radio  wave  after  refraction. 
If  the  planet  is  assumed  stationary  with  respect  to  the  earth,  the  additional  Doppler 
frequency  f  caused  by  atmospheric  refraction  through  an  angle  cv  is 

■Xf  =  k  ■  v  -  vy-tcoso  ll  *  v,cosd  sino  (2.5a) 

where  k  “  k/k.  and  where  0  is  the  (positive)  angle  between  v,  and  the  plane  of 
refraction.  If  v  is  lying  entirely  in  the  plane  of  propagation  1(1=0).  the  simpler  first 
order  form  of  equation  (  2.5a)  is 


Xf=vta  (2.5b) 

From  extrapolation  of  radio  tracking  data  during  times  when  there  is  no  atmospheric 
perturbations  of  the  railio  wave,  the  position  and  velocity  of  the  spacecraft  are 
known.  According  to  equations  (2.5a)  or  (2.3b/,  a  measurement  of  the  utmospherii 
Doppler  f  yields  the  bending  angle  a  and  hence  the  impact  parameter  a.  Front  corres¬ 
ponding  values  of  a  and  a,  the  differential  da/da  may  be  formed  and  hence,  upon 
performing  the  indic  ated  integration  in  equation  ( 2 .  J  b ) ,  the  refractive  index  correspon¬ 
ding  to  a  g  ceil  impact  parameter  is  obtained.  Application  of  the  Bouguer-relation 
(2.2)  finally  yields  the  desired  n(r0 (  profile. 

Calculation  of  a  refractivity  profile  from  stellar  or  radio  occultation  intensity 
measurements  follows  an  essentially  similar  path,  lhe  fundamental  equation  is  still 
(2.1b),  but  now  da/ da  is  expressed  in  terms  of  the  intensity  of  the  probing  signal. 
Indeed,  for  small  bending  angles  and  neglecting  focusing  of  the  wavefront  around  the 
curved  limb  (al)  <  R),  geometrical  optics  flux  conservation  yields  for  the  intensity  0 
of  the  near-limb  ray 

1  0'1  =  Dda/da  (2.4) 

when  0  is  normalized  to  unity  before  refraction.  Thus,  since  a(t)  is  known,  a  measure¬ 
ment  of  0  versus  time  t  establishes  da/da,  and  via  equations  (2.1b)  and  (2.2)  also  the 
refractivity  as  a  function  of  atmospheric  depth. 
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For  an  isothermal  atmosphere,  where  v  varies  exponentially  with  height  with  a  con¬ 
stant  scale  height  11,  equation  (2.4)  may  he  manipulated  to  yield 

^  =  (l-2)+,n(i-1)  (2.5) 

H  <p  <p 

Here  vq  is  the  component  of  the  spacecraft  velocity  in  the  plane  of  refraction  that  is 
perpendicular  to  the  ray  before  refraction.  The  time  t  is  taken  to  be  zero  when 
<P  =  1/2,  i  e,  corresponding  to  the  nominal  occultation  level.  An  essentially  identical 
form  of  equation  (2.5)  was  first  derived  by  Baum  and  Code  (1953)  and  has  since 
been  used  extensively  in  the  interpretation  of  stellar  occultation  measurements.  Thus, 
lor  an  isothermal  atmosphere  the  intensity-time  profile  <p( t)  yields  the  temperature  T 
directly  if  the  mean  molecular  mass  m  is  known,  since  T  =  mgH/K,  where  g  is  the 
acceleration  of  gravity  and  K.  is  Bolt/.man’s  constant. 

Equations  (2.4)  and  (2.5)  strictly  apply  only  to  a  point  source.  For  stellar  occulta- 
tions  the  intensity  <p  should  be  convolved  with  the  brightness  distribution  of  the 
projected  stellar  image  in  the  occulting  atmosphere.  As  shown  by  P  J  Young  (1977), 
however,  the  resulting  corrections  for  the  finite  angular  width  of  the  star  are  normally 
altogether  negligible 

Once  a  refractivity  profile  has  been  obtained  from  either  Doppler  or  intensity 
measurements,  the  hydrostatic  equation  may  be  integrated  to  yield  the  temperature 
profile  T(r0),  and  by  the  perfect  gas  law  also  the  pressure  profile  p(r0)-  In  a  well- 
mixed  atmosphere  of  non-polar  molecules  we  have  (eg,  Fjcldbo  and  Eshleman  1968) 

I  (r0 )  =  T(r0 , )  +  v  'V  \  ^  S(ro )  "(ro )  (2-6) 

v(r0)  kv(r0)r01 


P(r„)  =^---  "M  (2-7) 

tn 

In  these  equations  m  is  the  (constant)  mean  molecular  mass  of  the  atmospheric 
constituents,  v  is  the  mean  "refractivitv  volume”  (/  c.  the  refractivity  at  standard 
conditions  divided  by  Eoschmidt’s  number),  and  K  is  Boltzman’s  constant.  Equation 
(2.6)  assumes  a  boundary  condition  high  in  the  atmosphere  at  r0  =  r0,. 

Having  reviewed  the  formal  occultation  theory  it  is  appropriate  to  add  a  few  com¬ 
ments  relative  to  our  objective,  as  stated  in  section  1.2.  We  note  first  that  the 
procedure  described  above  assumes  the  atmospheric  refractivity  to  be  spherically 
stratified.  Turbulence,  representing  small  inhomogencities  in  v  over  a  continuum  of 
scale  sizes,  clearly  violates  this  provision.  One  might  think  that,  since  occultation 
experiments  arc  sensitive  to  atmospheric  properties  over  a  distance  C  =  (RH)12  in  the 
direction  of  the  ray,  the  calculated  refractivitv  profile  when  turbulence  is  present 
would  contain  an  additional  component  equal  to  the  average  turbulent  refractivity  over 
the  distance  C.  This  is  not  the  case,  however.  As  indicated  in  section  1.2,  turbulent 
refractivity  fluctuations  in  fact  produce  changes  in  signal  strength  that  grow'  with 
distance  and  may  eventually  become  comparable  to  the  mean  occulted  intensity,  even 
though  the  turbulent  refractivity  fluctuations  arc  a  small  fraction  only  of  the  ambient 
refractivity.  The  intensity  fluctuations  also  depend  on  the  radiation  wavelength  X, 
whereas  the  geometrical  optics  intensity  formulae  (2.4)  and  (2.5)  by  their  nature  arc 
independent  of  X.  We  thus  conclude  that  turbulent  refractivity  fluctuations  cause  the 
derived  refractivity  profile  to  depart  from  the  "true"  profile  of  the  ambient  atmo¬ 
sphere  in  ways  that  may  be  largely  independent  of  the  generic  fluctuations  in  refractivity 
at  the  ray  periapsis.  For  these  reasons,  exploring  the  ellects  of  turbulence  on  atmo¬ 
spheric  profiles  obtained  from  occultation  measurements  constitutes  both  a  legitimate 
and  important  scientific  endeavour. 
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2.2  Effects  of  boundary  conditions  and  various  measurement  errors 


For  practical  reasons  the  integration  in  equation  (2.1b)  has  to  be  truncated  at  some 
upper  limit  a’=  a) ,  say.  For  an  isothermal  atmosphere  the  resulting  fractional  error  in 
derived  refractivity,  bv,  is  approximately  (Wasserman  and  Veverka  1973) 

5„  =  1  -  erfl0(rol  -  r0)|  (2.8) 

where  |3'1  =  H  is  the  scale  height,  erf ( • )  is  the  error  function,  and  r0!  is  related  to  U| 
via  Bouguer’s  rule.  The  form  of  equation  (2.8)  implies  that  the  refractivity  error 
decreases  rapidly  with  increasing  penetration  of  the  rav;  for  r0)  r0  =  211,  bv  -  3%, 
while  three  scale  heights  below  the  starting  point  of  the  integration  &v  =  1 .5%.  F'rom 
equations  (2.6)  and  (2.7)  it  is  found  that  also  the  associated  temperature  and  pressure 
errors  decrease  rapidly  below  the  first  few  scale  heights  of  these  profiles. 


An  additional  error  in  F  or  p  arises  from  the  essentially  unknown  boundary  condition 
in  temperature.  From  equations  (2.6)  and  (2.7)  it  is  seen  that  an  error  AT(r0i  )  in  the 
boundary  temperature  yields  errors 


and 


AT(r0)  =  AT(r01) 


Ap(r0) 


_  Fcf(r0i 


m 


(2-9) 

^(r0  ) 

AT(r0l) 

(2.10) 

in  temperature  and  pressure,  respectively,  at  depth  r0.  While  the  temperature  error 
decreases  rapidly  with  depth  (eg,  exponentially  for  an  isothermal  atmosphere),  the 
absolute  pressure  error  is  constant  with  depth.  The  fractional  error  in  pressure 
decreases  rapidly  with  atmospheric  depth,  however. 


The  effect  of  different  types  of  measurement  errors  on  the  determination  of  p  and  T 
has  been  considered  by  Eshleman  ( 1 973),  and  in  part  also  by  Hubbard  ct  al  (1975). 
Following  Eshlcntan  (1975),  the  temperature  scale  height  of  an  assumed  isothermal 
atmosphere  may  be  expressed  in  terms  of  spacecraft  trajectory  data,  the  radio 
I  hippier  f,  and  its  time  derivative  df/dt  as  in 


MD  (^ll__i)M 

V1H  XD  |— | 

dt 

Alternatively  this  relation  may  be  written  as 


(2.11) 


M(S  -  1)  =  1  (2.12) 

where  the  dimensionless  magnitude  factor  M  and  the  sensitivity  factor  S  arc  defined 
by  comparison.  Since  \1  for  planetary  occultations  may  reach  a  hundred  or  more.  S 
only  slightly  exceeds  unity  over  most  of  the  atmospheric  region  probed.  As  pointed 
out  by  Eshleman  (1975),  S  is  identically  unity  for  bending  uiouml  a  fixed  limb 
(11=0),  as  in  classical  knife-edge  diffraction.  All  information  about  atmospheric  struc¬ 
ture  is  therefore  contained  in  the  small  excess  of  S  over  unity. 

Referring  to  equations  (2.11)  and  (2.12).  absolute  errors  AS  in  S  and  A(aD)  = 
A(Xf/V| )  in  the  deflection  of  the  rav  produce  a  fractional  error  in  11  given  b\ 

AH=MAS  +  ^t>»  (2.111 

II  OD  1 

Thus  fractional  errors  in  any  of  the  quantities  vf ,  I),  X  or  Idf/dtl  comprising  S  arc 
magnified  by  the  potentially  large  factor  M  =  ftl)/H,  while  a  fractional  error  in  a  (re, 
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in  f)  produces  only  a  comparable  error  in  H.  If  !df/dt|,  D  or  X  equals  (1— A)  times  its 
true  value,  or  if  the  derived  Vj2  is  (1+A)  times  its  true  value,  then  to  first  order  in  A 
and  for  large  M  (Eshleman  1975) 


Pf 


=  1  +  MA 


r  3/2 


(2.14a) 


(2.14b) 


where  I  f  and  l't  denote  the  false  and  true  temperature  profiles,  respectively,  and  with 
a  similar  notation  for  p. 


Despite  the  potentially  large  error  in  derived  profiles  induced  by  small  errors  in  S, 
temperature/pressure  profiles  derived  from  Doppler  measurements  are  still  more  accu¬ 
rate  than  corresponding  profiles  derived  from  radio  or  stellar  occultation  intensity 
measurements,  provided  that  the  radio  Doppler  can  be  measured  at  least  M  times 
more  accurately  than  the  intensity,  which  is  normally  the  case.  From  equation  (2.5) 
it  is  indeed  found  that  errors  in  either  v^  or  0  yield  comparable  errors  in  the  scale 
height  determination,  and  hence  in  T.  This  is  also  the  reason  why  stellar  occultations 
are  routinely  terminated  when  the  residual  light  flux  is  comparable  to  the  expected 
error  in  baseline  determination,  which  is  typically  a  few  per  cent. 

Digitization  and  subsequent  data  processing  also  introduce  errors  in  the  final  T— p 
profiles.  Such  errors  arc  highly  dependent  on  signal  to  noise  ratio  during  data  trans¬ 
mission,  and  on  the  specifics  of  the  data  processing  algorithms  used,  as  discussed  in 
detail  by  l.ipa  and  Tyler  (1979). 

The  previous  analysis  apparently  allows  the  following  arguments  to  be  made  concern¬ 
ing  the  effect  of  turbulence  on  occultation  profiles.  Turbulence  causes  the  actually 
observed  bending  angle  a  to  deviate  slightly  from  the  "true”  value  of  the  quiescent, 
non-turbulent  atmosphere.  According  to  equation  (2.13).  a  change  in  or  produces  a 
comparable  (t  c.  small)  error  in  II  and  hence  in  T  if  these  quantities  are  obtained 
from  radio  Doppler  measurements.  Atmospheric  profiles  derived  from  intensity 
measurements,  on  the  other  hand,  contain  errors  that  are  comparable  to  those  in 
intensity  according  to  equation  (2.5).  Since  intensity  scintillations  due  to  turbulence 
may  become  comparable  to  0  itself,  it  follows  that  the  resulting  temperature  error 
may  also  be  appreciable,  unless  some  "smoothing”  procedure  is  applied  that  suppresses 
the  scintillations.  We  shall  see  from  the  more  stringent  analysis  in  section  5  that  these 
contentions  are  indeed  true. 


3  ELECTROMAGNETIC  WAVE  PROPAGATION  IN  A  NON-TURBULENT  PLANE¬ 
TARY  ATMOSPHERE 

3.1  General  formulation 

Standard  formulations  of  electromagnetic  wave  propagation  in  a  spherically  sym¬ 
metric,  non-turbulent  atmosphere  arc  based  on  geometrical  optics.  The  rationale  for 
this  is  the  tact  that  the  radius  of  the  first  Fresnel  zone  is  smaller  than  the  atmo¬ 
spheric  scale  height  by  typically  a  factor  on  the  order  of  10,  cf  Table  3.1.  We  depart 
here  from  the  common  geometrical  optics  approach  by  formulating  the  propagation 
problem  in  a  wave-optical  context,  from  which  the  proper  approximations  leading  to 
the  standard  geometrical  optics  results  will  appear  in  a  natural  way.  The  small  wave- 
optical  corrections  to  these  results  arc  also  readily  obtained  from  this  approach, 
which  also  provides  the  notational  framework  for  the  generalization  in  section  4  of 
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MARINER  10/VENUS 

10 

1.5  -104 

1.4 

MARINER  4/MARS 

9 

2.6  •  104 
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30 

5.6  10s 

7.3 

|3Sco  A/JUPITER 

30 

6.5  ■  10s 

0.57 

BD17°43 88/NEPTUNE 

60 

4.5  10y 

1,5 

Table  3.1  Representative  values  for  various  radio  and  stellar  occultations 

The  table  gives  representative  values  for  the  stale  height  11,  approximate  plancl-to- 
spaeeeraft  (or  earth)  distance  If,  radius  of  the  frec-space  Fresnel  zone  up  .  and 
projected  linear  radius  as  p  of  two  occulted  stars.  All  lengths  are  in  kilometers.  The 
two  stars  are  observed  at  A  =  O.a/im.  while  X  =  0.1  3  m  (S -band)  has  been  assumed  in 
the  calculation  of  up  p  for  the  radio  occ  ultations. 


PHASE  SCREEN 


PLANE  WAVE 


Figured.!  .1  phase  changing  screen  oj 
thickness  C  <4  P  represents  the 
inhomogeneous  ambient  atmo¬ 
sphere 

Position  vectors  I,  and  r2  refer  to  the  exit  fate  of 
the  screen  and  an  arbitrary  field  point  a  distance 
i r j i  3  1)  C  from  the  screen,  respei  lively 


//KU, 


exp|ik  lr2 

ir2  r, | 


the  weak  scattering  scintillation  theory 
to  include  the  inhomogeneous  ambient 
atmosphere. 

We  shall  assume  throughout  that  the 
effect  of  the  actual  atmosphere  on 
electromagnetic  waves  may  be  repre¬ 
sented  by  a  thin  slab  of  thickness  C  <  I) 
md  refract ivity  f0,  the  effect  of  which 
is  to  change  the  phase  of  the  wave  but 
leaving  its  amplitude  unchanged,  see 
Figure  3.1.  We  shall  find  that  this  thin 
screen  approximation  is  always  valid 
for  ray  bending  angles  a  <S  (H/R),/J , 
where  II  is  a  typical  scale  of  variation 
of  the  spherically  stratified  atmosphere. 

Let  F.(r, )  be  the  complex  electric  field 
of  the  wave  in  the  plane  x  =  C/2  upon 
exit  from  the  phase  screen.  In  position 
r2  behind  the  screen  the  field  is(<‘£, 
Goodman  1968) 


cos(n.c)  dr, 


where  dru  -  dy,d/.,  is  a  differential  area  on  the  exit  plane  x  =  C/2  of  the  screen.  The 
obliquity  factor,  cos(n.e),  accounts  for  the  angle  between  the  normal  n  to  the  phase 
screen  and  the  direction  e  =  (r2  r,  )/lr2  r,  I  to  the  observer.  If  the  spatial  extent  of 
the  refracting  medium  perpendicular  to  the  ravpath  is  small  compared  to  D.  the 
angle  Ln,c'  0  and  cos(n.e)  '  1,  while  lr2  r,  I  3  1).  This  paraxial  approximation  is 
assumed  throughout.  Consistent  with  the  thin  screen  approximation,  an  incident  plane 
wave  of  unit  amplitude  has  an  electric  field  at  exit  given  by  F(r,  )  =exp(ikS0),  where 
S0  is  the  change  in  phase  path  through  the  screen.  We  shall  see  in  section  3.2  that, 
consistent  with  the  thin  screen  approximation,  the  phase  change  may  be  calculated  by 
integrating  e0  along  a  straight  line  path  through  the  screen,  instead  of  following  the 
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actual  curved  raypath.  Thus  S0  =  fv0{x,y0,Zo  )dx,  where  (y0,z0)  arc  lateral 
coordinates  of  the  ray  upon  entering  the  screen.  If  6  -  (27fRH)|/2  ,  it  also  follows  that 
S0  =  i>06,  where  r’o  is  now  referred  to  the  ray  pcriapsis  of  the  actual  exponential  atmo¬ 
sphere  of  scale  height  11. 

With  these  provisions  equation  (3.1 )  may  he  expressed  in  t he  form 
i  •>  <  —  i j/ 

E(r2)  =  -^-//e  20  dru  (3.2a) 

2ttiD  -« 


where 


*  =  2D(L)+  S0)  +  (y,  y,)2  +  (z2-  z,)2 


(3.2b) 


Following  Haugstad  (1978c)  (sec  also.  Hubbard  and  Jokipii  1 977b),  we  expand  \p 
around  the  point  of  stationary  phase  (y0 .z0 )  of  the  integrand,  i  c 


$  =  -  cy(y ,  ->•„)'(/.,  --z„)-' 

i,j  =  o  'J 


(3.3a) 


where 

=  J_  d(i+i>ii/ 

Sj  i!J'!  d\J  dj  ;T*° 

*1  -t-n 

The  stationary  phase  [joint  is  defined  by  the  conditions  c0 1  (y  i  =y0 , 
cio(Vi=y0.  z,=z0)  =  0,  yielding 

Vo  =  Y:  -  l>0  .  z0  =  z2  -  Oa 


(3.3b) 
Z|=Z0)  = 


where  a  and  (3  are  the  ray  bending  angles  in  the  xz-  and  xy-planc,  respectively,  as 
given  by  a  =  JV0,/  (x,y0.y0)dx  and  (3  =  ft^0  y  (x,y0.z0)dx.  (Subscripts  v.  z,  yy,  yz,  etc 
denote  partial  differentiation  in  the  usual  way.)  Successive  determination  of  the  re¬ 
maining  coefficients  Cy  yields,  through  second  order 

Coo  =  21)2  [  1  +i(a2  +  |32)]+  2DS0 

cot  =  cio  =  0  (3.4) 


c02  1  +  M  ,  C| ,  -  0  ,  c20  -  1  +  N 


I  he  two  quantities  M  and  N  play  a  central  role  in  the  subsequent  analysis;  they  are 
defined  by 


M  - /t'o.zz(x,y0.z0)dx  ,  N  - /t'o,yy(x,yo,z0)dx  (3.5) 

where  the  integration  is  again  understood  to  be  through  the  screen  from  x  =  —6/2  to 
x  =  6/2. 

Equations  (3.2)  (3.5)  constitute  the  basic  mathematical  framework  for  the  sub¬ 
sequent  analysis  in  this  and  also  in  the  next  main  section. 


3.2  Geometrical  optics 

We  shall  assume,  tentatively,  that  the  stationary  phase  expansion  (3.3)  can  be  trun¬ 
cated  after  the  second-order  terms,  i  e,  we  define  Cy  ( i  +  j  >  2)  =  0.  In  this  case  equa¬ 
tion  (5.2a)  becomes 
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,k^-° 

2D 


2  ff  i  D  -~ 


//e 


(0y  (yi  -yo)2  +  <t> 2  (Z,  — Zo  )2  I 


dr 


lx 


(3.6) 


where  we  have  defined 

c20  =  0y  =  1  +  N  ,  c02  h  0Z>  =  1  +  M  (3.7) 

The  integrals  in  equation  (3.6)  are  of  the  Fresnel  type  when  evaluated  between  the 
limits  ±  00  as  indicated.  Apart  from  possible  contributions  from  the  neglected  terms  in 
the  stationary  phase  expansion,  the  very  presence  of  the  planet  itself  requires  that  a 
circular  region  of  radius  R,  centered  at  Z|  =  r0,  be  excluded  from  the  integration. 
However,  due  to  the  rapidly  decreasing  contribution  to  the  integral  from  wavelets 
originating  from  positions  |z,  —  z0l  >  (XD)1/2  0Z'/2  and  lyi  yol  >  (XD/2  0'^  from 
the  stationary  phase  point,  the  presence  of  the  planet  can  be  ignored.  Evaluating 
equation  (3.6)  for  the  intensity  0  =  IE|2  thus  yields 

0  =  0y0z  (3.8) 

This  result,  independent  of  k,  is  the  geometrical  optics  intensity  formula  for  a  spheri¬ 
cally  symmetric  refractivity  field.  It  can  be  derived  from  a  flux  conservation  argument 
directly,  or  from  a  geometrical  optics  formulation  based  on  the  Eikonal  equation 
(Haugstad  1978e).  Truncation  of  the  stationary  phase  expansion  after  the  second- 
order  terms  therefore  restricts  the  validity  of  the  resulting  intensity  formula  to  the 
realm  of  geometrical  optics.  We  calculate  in  the  next  section  the  lowest-order, 
k-dependent  correction  to  equation  (3.8),  and  thereby  also  establish  the  region  of 
validity  of  this  result. 

Equation  (3.6)  provides  an  interesting  wave-optical  interpretation  of  ray  bending  in 
an  inhomogeneous  medium.  Indeed,  the  phase-factor 

0  =  ^l0y  (y  t  ->’o)2  +  4>z  (Z|  -Zo  )2  ] 


has  the  constant  value  tr  on  the  circumference  of  an  ellipse  centered  at  (y0.Zo)  with 
semi-axes 


aFJl  =0|.2(XD),/2  ,  aKv  =  0>/2(XD)''2  (3.9) 

in  the  y-  and  z-directions,  respectively-  Refraction  thus  distorts  the  originally  circular 
shape  of  the  frce-space  Fresnel  zone  of  radius  ap ()  =  (XD)12  into  an  ellipse  with 
semi-axes  given  by  equation  (3.9).  The  property  that  the  intensity  is  proportional  to 
the  square  of  the  area  of  the  first  Fresnel  zone  is  preserved,  however,  since 
(trap  h3 ?,v/ndf2,0  )2  =  0y0z>  *n  agreement  with  equation  (3.8).  We  note  also  that 
diffraction  effects,  through  the  inclusion  of  higher-order  terms  in  the  stationary  phase 
expansion  and  in  0  above,  will  cause  the  atmospheric  Fresnel  zone  to  deviate  slightly 
from  its  strictly  elliptical  shape  in  the  geometrical  optics  limit. 

For  a  spherically  symmetric  and  exponential  atmosphere  the  definitions  (3.5)  may  be 
used  to  show  that 

M  =  aD/H  ,  N  =  -aD/R  (3.10) 

in  this  case.  The  corresponding  result  for  the  intensity  is 

0’ 1  =  (1  +  5(D)  (1  -5®) 

'  H  ’ '  R 


(3.11) 
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which  reduces  to  equation  (2.4)  when 
focusing  is  negligible  (al)  <  R).  In  the 
present  case  refraction  expands  the 
frce-space  Fresnel  zone  by  a  factor 
( 1 -al)/R)',/2  >  1  parallel  to  the  limb, 
while  at  the  same  time  compressing 
the  vertical  scale  by  the  factor 
( 1  +  aD/ll)1'2  <  1,  see  Figure  3.2.  As 
aD  -*  R,  the  spacecraft  approaches  the 
center  of  the  planet,  as  viewed  in  the 
plane  of  the  sky,  and  equation  (3.1 1) 
predicts  Diffraction  effects 

limit,  however,  the  intensity  to  a  finite 
value  whose  magnitude  is  readily  ob¬ 
tained  from  the  physical  picture  de¬ 
veloped  above  (cf,  Kshleman  ct  al 
1979a). 

Fhc  phase  in  the  geometrical  optics 
approximation  is  obtained  from  equa¬ 
tion  (3.6)  by  writing  the  electric  field 
as  F.  =  A  exp  |i(kS-uu)|,  where  A  is 
the  amplitude  and  kS  is  the  phase.  Per¬ 
forming  the  integration  over  the  screen 
it  is  found  that 

kS  ~  k|l)  ‘  S0(y„ .•/.„)  t  *  tf2M  (3-12) 

1  he  sum  kl)|  1  t  l/2(or  t(32)|  is  simple  the  increase  in  phase  by  following  the  geomet¬ 
rical,  refracted  path  to  the  receiver,  while  kS0  is  the  additional  atmospheric  phase 
perturbation  inside  the  refractive  medium. 

The  previous  developments  are  all  based  on  the  thin  screen  approximation,  as  defined 
in  section  3.1.  liquation  (3.8),  or  more  conveniently  equation  (3. 1  1 ),  may  now  be 
used  to  explore  the  validity  of  this  approximation.  For  this  approximation  to  be  valid 
it  is  sufficient  that  the  change  in  intensity  <t>  through  the  medium  be  small.  Thus, 
putting  I)  =  V  5  (Rtl)1  2  in  equation  (3.1  1 ),  we  find  that  the  thin  screen  approxima¬ 
tion  is  satisfied  if 

a  <  11/^5  ( 1 1/ R ) 1  2  (3.13) 

That  is,  the  deflection  of  the  ray  inside  the  refractive  medium  should  be  a  small 
fraction  of  a  scale  height  only.  The  integration  of  v0  along  the  actual  curved  raypath 
to  obtain  the  atmospheric  phase  perturbation  can  therefore  be  approximated  by  inte¬ 
gration  along  a  straight  line  path,  as  already  assumed. 

The  seeminglv  restrictive  condition  (3.13)  is  found  to  have  been  satisfied  for  most 
past  radio  occultations.  It  is  invariably  satisfied  in  stellar  occultations  where  the 
optical  rays  suffer  much  smaller  bending.  Using  a  direct  flux  conservation  argument 
not  restricted  by  the  thin  screen  approximation,  F.shleman  <7  al  ( I  979a)have  derived  a 
generalization  of  our  equation  (3.1 1)  that  reduces  to  this  result  when  the  condition 
(3.13)  is  satisfied. 


o 

1 


Figure  3.2  Schematic  illustration  of  the 
changing  shape  of  the  atmo¬ 
spheric  Fresnel  zone  as  the 
observer  approaches  the  focal 
line  in  the  plane  of  the  sky 
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3.3  Wave-optical  effects 

Since  the  tree-space  Fresnel  zone  is  a  small  hut  finite  fraction  of  atmospheric  scales 
of  variation,  the  true  phase  and  intensity  of  the  radio  wave  deviate  slightly  from  the 
geometrical  optics  results  of  the  previous  section.  Close  to  the  focus,  where  aFh 
increases  without  bound,  wave-optical  effects  eventually  become  dominant  for  any 
finite  value  of  the  wavenumber  k,  as  noted  previously.  We  calculate  below  the  magni¬ 
tude  of  these  diffraction  effects  in  the  two  limits  where,  (i)  a  F  k  =  aj.  0,  but  with  no 
additional  restriction  on  aK  v ,  and  ( ii )  when  at  h  ie,  when  the  intensity  is 

critically  controlled  by  focusing  of  the  wave  around  the  curved  limb. 


The  former  of  these  cases  is  most  easily  treated  by  expanding  the  phase  factor  0 
beyond  the  quadratic  terms  and  retaining  the  lowcst-order  correction  to  the  electric 
field  E.  Writing 


0  -  [0y  (v ,  -y0)2  +  0zl(zi  -z0)2 


coj(zi  -Zo)’  +  C,2(z,  - z0  )2  (y  i  ~yo )  + 


(3.14) 


it  is  found  that  the  leading,  finite  wavelength  correction  to  the  electric  field  involves 
terms  of  the  form  (z,-z0)4  and  (z,  -z0)8.  Putting  E  =  E0(1+6E),  where  the  geomet¬ 
rical  optics  field  E0  =  0V|/20Z1,2  exp(i(kS— cot)]  and  5E  is  the  lowest-order  fractional 
correction  to  E,  we  find  after  considerable  algebra  that 


5E  =  -  -1-  0.,(  1-0., 

32  z  *• 


J32  D 


--  <t>2,(  1 

384 


(?  IT2 
k2 


(3.15) 


where  we  have  defined  j3  -  1/H  as  the  reciprocal  scale  height.  Thus  there  is  a  correc¬ 
tion  kSS  to  the  geometrical  optics  phase  by  an  amount 


d2I) 


k5S=-,^0z(l  -0z)tJ 


(3.16) 


There  is  also  a  fractional  intensity  correction  given  by 


50 

0 


35 

192 


0*(1 


0  )2  ^ 

9/-  k2 


(3.17) 


This  latter  result  is  functionally  similar  to  a  corresponding  result  by  Hubbard  (1979) 
for  propagation  in  two  dimensions. 


Since  the  atmospheric  phase  perturbation  kS0  =  kf0^-  klla,  it  follows  from  equa¬ 
tion  (3.16)  that  the  fractional  change  in  phase  due  to  diffraction  is 


k6S^  1  ,:daI>2 

kS0  32  vz  k2 


(3.18) 


That  is,  there  arc  fractional  corrections  in  phase  and  intensity  that  arc  both  propor¬ 
tional  to  K^j34D2/k2  and  that  decrease  rapidly  with  occultation  depth.  Since 
x  spans  a  range  of  order  of  magnitude  10'4  10’8  for  past  radio  and  stellar  occulta- 

tions,  these  corrections  will  normally  be  altogether  negligible. 


Sufficiently  close  to  the  focal  line  «L)  =  R,  diffraction  effects  eventually  become 
dominant  for  any  Finite  value  of  k.  I.ctting  rx  denote  the  perpendicular  distance  of 
the  spacecraft  from  the  focal  line,  Hubbard  (1977)  has  shown  that  for  an  exponential 
atmosphere  the  intensity  close  to  the  foclis  is  given  by 


0 


4*2  S?J* 


2tr  rx  R 
XI) 


) 


(3.19) 
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where  J0  is  a  Bessel  function  of  zero  order.  The  focal  intensity  maximum  is  thus 

0max  =  47f2  RH/XD  (3-20) 

with  a  characteristic  width 

w  =  \D/R  (3.21) 

For  w  <<  rx  <  R,  the  average  intensity  reduces  to  the  geometrical  optics  result  (3.11) 
when  allowance  is  made  for  a  doubling  of  the  intensity  close  to  the  focal  line  where 
the  two  limbs  yield  nearly  identical  contributions  to  the  average  intensity.  Since  w  is 
typically  a  small  fraction  of  a  meter,  the  ray-optical  intensity  formula  (3.11)  is  valid 
extremely  close  to  the  focal  line.  As  shown  by  Eshleman  et  al  ( 1 979a),  the  results 
(3.20)  and  (3.21)  may  also  be  derived  from  simple  physical  arguments  relating  to  the 
limiting  distortion  of  the  atmospheric  Fresnel  zone  as  the  focal  line  is  approached. 


4  PROPAGATION  IN  A  TURBULENT  PLANETARY  ATMOSPHERE 


4.1  Phase  and  intensity  in  the  weak  scattering  approximation 


The  development  of  the  previous  section  provides  a  framework  for  calculating  the 
effect  of  turbulence  on  electromagnetic  wave  propagation.  In  the  thin  phase  screen 
approximation  turbulence  merely  adds  a  component 


S,  = 


1 12 


"i (x.  y0. z0)  dx 


(4.1) 


to  the  total  optical  path  length  at  the  exit  face  of  the  screen.  It  is  assumed  through¬ 
out  that  the  turbulent  component  of  refractivity,  t>, ,  has  zero  mean  and  satisfies 
<i'2>1'2  <  tv  Equation  (3.1)  takes  the  form 

_  i^(2D2+(y2-yI)2  +  (z2-zI)2  +  2I)(S0+S1)] 

E(r2)  =  -p-  ffc  ~D  dr)1  (4.2) 

2tnD 


In  the  weak  scattering  approximation  kiS ,  1  <  1.  and  the  integrand  in  equation  (4.2) 
may  be  expanded  around  the  point  of  stationary  phase  as  before.  Expanding  the 
electric  field  E  =  Fi0  +  E,  +  ....  where  F>0  =  (9(l  +  t'0),  E,  =  (?(tq),  etc,  we  thus  obtain 


E|  (r2 )  =  k0v0ze 


ik^ft/2 
2D  ' 


/  vi  (x,  v0,  z0)dx  ffc 


*2i)l<^y  (y*  ~‘y°)2+ ^z  (z> 


-Zol 


dr, 


(4.3) 


after  retaining  only  the  quadratic  terms  in  the  stationary  phase  expansion.  Let 


'-.)=///  c 


ijuxx,  +  uv(y,  -y0)+  uz( 


i  —  z0  )  ] 


dZ(u) 


(4.4) 


be  the  Fourier-Stieltjes  representation  of  vt ,  where  dZ(u)  =  Z'(u)du  (du  = 
duxduyduz),  with  ux,  Uy,  uz  being  the  spatial  wavenumber  components  in  the  x-,  y-, 
and  z-dircctions,  respectively.  The  Fourier  amplitude  Z'  of  v,  has  been  referenced  to 
the  stationary  phase  point  (y0,  z0).  For  a  homogeneous  and  isotropic  random  field 
<Z  (u,  )Z’,l(u2)>  =  4>(U|  )  5(u,  u2),  where  ‘I’  is  the  spatial  power  spectrum  of  and 

5(-)  is  the  Dirac  delta-function.  Combination  of  equations  (4.3)  and  (4.4)  yields  after 
algebraic  manipulations 

ik  l  ^ 

E,  (r2 )  --  k0y0ze  21)  fffsinci )  /.'(u)  i2(uA)  du  (4~r>) 
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Merc  sinc(-)  is  the  sinc-function,  ux  =  (uy,  uz)  and  we  have  defined 

2k  lVv  f  1  f  ‘  K'  y“  +  z° )  -  '  9T)^y  y'2  +  0z‘  Z|2 1 

“  ~  °c  ~  ’  ■  f f  e  “  dy'i  dzj 


(4.6) 


The  first-order  perturbation  in  intensity  normalized  by  the  average  intensity,  0|  ,  and 
the  first-order  perturbation  in  phase  path,  S,,  are  related  to  Kt  as  in 

0 i  =2Re(Ko‘K|  i,  S,  =  k'1  Im|Ko  Fi,  I 

respectively.  Since  F0  ~  0V  0Z  exp(ik  c0o/-M),  we  obtain  after  integration  over  the 
phase  screen  in  equation  (4.6) 


0,  =  2kV  /// sinc(Uj^-)  /-'(u)  sin  (P)  du 

2rr 


(4.7] 


u  t  , 

Sj  =  C  fff  sinc(  x  )  X  (u)  cos  (P)  du 

2tr 


(4.8) 


where  we  have  defined 

P  =  ^  (0yuy  +  0zuj)  (4.9) 

For  the  mean  square  of  these  quantities  we  obtain,  assuming  I,0  <  V.  (1.0  is  the  outer 
scale  of  turbulence) 

<0?>  =  87rtk*  //  <l»(ux)  sin2  (P)  dux  (4.10) 

<Sj>  =  27rC  //  <l>(uj  cos:(P)  dux  (4.111 

1  he  quantity  <0:,>  is  called  the  scintillation  index,  the  square  root  of  which  is  often 
termed  the  modulation  index. 


liquations  (4.10)  and  (4.11)  are  the  appropriate  generalizations  of  the  weak  scatter¬ 
ing,  mean  square  fluctuations  in  intensity  and  phase  path  to  turbulence  superimposed 
on  an  inhomogeneous  background.  (For  stellar  ocoultations  the  intensity  should  be 
convolved  with  the  projected  stellar  brightness  distribution  in  the  occulting  atmo¬ 
sphere,  sec  section  4.8.)  The  effect  of  the  ambient  atmosphere  is  seen  to  appear  onlv 
through  the  quantities  0V  and  0Z.  whose  interpretation  in  terms  of  the  mean  signal 
intensity  was  established  in  the  previous  section.  A  sufficient  condition  for  neglecting 
coupling  to  the  inhomogeneous  background  atmosphere  is  therefore  that  0V,  0Z  both 
be  essentially  unity,  or  equivalently,  that  the  change  in  intensity  caused  by  the 
ambient  atmosphere  alone  be  negligible.  This  condition  is  not  satisfied  in  cither  radio 
or  stellar  occultation  experiments,  however,  implying  that  interpretation  of  weak 
scintillation  data  in  such  experiments  should  proceed  from  the  generalized  formulas 
derived  here.  Surveillance  of  the  earth's  upper  atmosphere  using  a  pair  of  satellites 
may  also  involve  sufficient  ray  bending  for  weak  scintillation  data  to  be  analyzed  in 
terms  of  the  results  given  here,  especially  if  the  raypath  traverses  regions  of  thermal 
inversion  where  the  local  value  of  i)‘  v0/dz? ,  contained  in  the  definition  of  0Z>  may 
attain  appreciable  values.  It  may  be  noted  that  although  the  previous  derivation  has 
assumed  e,  to  be  a  strictly  homogeneous  anti  isotropic  random  variable,  the  results 
(4.10)  and  (4.11)  arc  still  valid  if  e,  is  inhomogeneous,  provided  that  the  scales  on 
which  the  inhomogeneitics  occur  arc  large  compared  to  any  other  scale  of  interest, 
i  e,  to  I.0  and  the  Fresnel  scales  a(,  h,  a,  . 

Having  established  the  first-order,  weak  scattering  results  for  phase  and  intensitv 
appropriate  to  radio  and  stellar  occultations,  it  is  important  to  explore  their  validity. 
For  the  ambient  atmosphere  the  thin  screen  approximation  was  found  to  require 


a«(H/R)12.  Since  the  first-order  intensity  and  phase  perturbations  are  both  func¬ 
tions  of  <i>y,  <t>y.  the  condition  on  the  smallness  of  a  also  applies  in  this  case.  1  he  thin 
phase  screen  approximation  also  requires  that  the  scintillations  in  intensity  be  small 
inside  the  screen.  This  condition  is  automatically  satisfied,  however,  since  we  shall 
prove  in  section  4. .3  that  <0J>  is  small  compared  to  unit)  at  any  distance  from  the 
screen  when  klS1(r,)|<^  1,  i  e,  when  the  turbulent  phase  perturbation  at  exit  from 
the  screen  is  small.  We  observe,  finally,  that  since  C  <  1)  the  quantities  0V  and  <t>z  do 
not  change  appreciably  through  the  screen  and  may  hence  be  associated  with  a  fixed 
spacecraft-limb  distance  x  =  1). 


The  stationary  phase  expansion  also  imposes  certain  restrictions  on  the  validity  of  the 
previous  results.  These  limitations  are  most  easily  appreciated  by  noting  from  the 
derivation  of  equation  (4.3)  that  if  Ay  s  y,  y0  and  Az  =  z,  z0.  a  Fourier  compo¬ 
nent  of  wavenumber  ux  =  (uv.uz)  gives  a  contribution  proportional  to 


I0y  I  Ay  t  ^  0yuy)2  +  <*>:}  (Az  +  ^  0zuz)2  ] 


iAv  dAz 


(4.12) 


to  the  perturbed  electric  field  at  the  receiver.  Defining  uv  =  2tr/fcv  and  uz  ~  27t/Cz,  it  is 
seen  front  this  expression  that  an  eddy  of  wavenumber  components  uv,  uz  gives  a 
contributiiin  t > >  the  perturbed  electric  field  that  originates  from  a  position 
(iA\,  -  0V  M)TV,  ,Az!  =  0Z  Xl)/t.  )  relative  to  the  stationary  phase  point.  Such  cle- 
mentary  wave  contributions  are  thus  scattered  through  angles  0y  -  0vX/Cv  and  d7  = 
0Z  A,TZ  in  the  xy  ami  xz  planes,  respectively.  When  0V  -  0Z  =  1,  this  is  essentially 
the  Bragg  condition  for  diffraction  by  a  sinusoidal  spatial  diffraction  grating.  Now. 
for  Kolmogorov  turbulence  <S2>  and  <02>  may  be  shown  to  derive  their  principal 
contributions  from  eddies  of  dimension  (liv,  tz)  ~  1.0  and  (Cv,  5Z)~  (ap  j,,  ap.v),  res¬ 
pective!).  f  he  quadratic  approximation  of  the  stationary  phase  expansion  should 
therefore  be  validated  to  distances  Ay  ~0vXI)/I,o.  ,Azi  ~  0Z  Xl)/L0  and  I  Ay  |  ~ 
ayh,  [Azi  '-ay  v  for  these  two  quantities,  respectively.  F.  valuation  of  the  next-order 
terms  in  the  stationary  phase  expansion  indeed  turn  out  to  be  small  at  these  distances 
if  H  I,„  ^  ay  v  ,  and  provided  also  that  0V  <  (I.0/ay  ,,)  (L0  R/ay  [))'7  in  the  evalua¬ 
tion  of  <Sf  -.and  pN .<  (RID1  2 /ay  n  in  the  evaluation  of  <0(>.  In  shallow  occulta- 
lions  i  hese  constraints  are  alwav  s  met .  Inn  very  close  to  the  focal  line  both  conditions  are 
violated  and  neither  <.$2>  nor  <S2>  will  be  given  by  simple  formulas  of  the  type 
(4.10)  and  (4.1  l). 

We  are  now  prepared  to  disc  uss  another  quantity  of  principal  interest  in  radio  occul- 
tation  experiments:  the  ray  bending  angle.  It  is  readily  shown  from  the  previous 
equations  that  a  major  component  of  the  mean  square  perturbation  in  bending  angle 
involves  an  integral  of  the  form 

f  f  <t>(  Uj )  ii2  cos2  ( 1’)  duy  (4.13) 

For  Kolmogorov  turbulence  this  integral  diverges  as  a  function  of  the  upper  limit,  and 
is  thus  controlled  bv  the  smallest  eddies  of  dimension  -  li0.  the  inner  scale  of  the 
turbulence.  Bv  relation  (4.12),  eddies  of  this  size  yield  scattered  waves  that  originate 
from  distances  cl  --  M)/t0  from  the  stationary  phase  point.  With  C0  1  cm,  d  is 
commensurable  with  the  dimensions  of  the  pi anet.  Fhe  implication  is  therefore  that  a 
calculation  of  the  perturbation  in  ray  bending  angle  requires  the  complete  spatial 
extent  of  the  scattering  medium  to  be  stipulated,  and  the  stationary  phase  expansion 
cannot  lie  truncated  after  the  second-order  terms.  The  problem  is  compounded  by 
the  fact  that  in  going  from  equation  (3.1)  to  (.3.2),  me  exact  kernel 


cxp(|k|r?  r,i| 
It  fi  I 
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of  equation  (3.1)  was  approximated  by  the  Fresnel  kernel 


—  exp( i  — ■  p2 
D  f  2D 


which  requires  that  x  k  p4/ 1)3  «  1,  where  p  -s  ( ( \  2  y,  ):  *  (/.,  y.j  I  ]1  ;  is  the  lateral 
separation  of  the  receiver  and  an  arbitrary  point  on  the  exit  face  of  the  phase  screen 
( Tatarskii  1971,  pp  223).  For  the  perturbations  in  intensity  and  phase  this  condition 
is  strongly  satisfied,  since  for  these  quantities  p  (M))i:  and  <  M)/1.0,  respectively. 
However,  the  fluctuations  it;  ray  bending  angle  derive  their  principal  contribution 
from  waves  scattered  from  distances  1  from  the  stationary  phase  point,  where  fc  is  the 
effective  lateral  extent  of  the  turbulent  medium.  For  turbulence  that  varies  exponen¬ 
tially  with  altitude  on  a  scale  clt  (cf  section  4.2),  the  medium  extends  effectively  to  a 
distance  p  =  t  s  (RHt)'  :  parallel  to  the  limb.  Assuming  llt  -  11  =  10  km.  we  find  that 
for  X-  and  S-band  oecultal ions  of  I'S  spacecraft  typical  values  of  x  arc  such  that 
x  S*  1 .  Hence  a  numerical  approach  is  apparently  required  to  obtain  the  mean  square 
fluctuations  m  bending  angle 


Most  previous  analyses  of  scintillation  data  from  radio  and  stellar  occult  ition  experi¬ 
ments  have  ignored  the  coupling  effect  between  the  turbulence  and  the  ambient 
background  described  and  calculated  here  (Gurvich  1 969,  Golitsyn  and  Gurvich  1971, 
Woo  and  Ishimaru  1973,  Woo  etui  1974,  Woo  and  Ishimaru  1974,  Woo  1973a, 
Jokipii  and  Hubbard  19771.  The  first  attempt  to  account  for  the  coupling  effect  on 

the  scintillation  index  was  apparently  made  by  Woo  etui  (1973).  Later  and  indepen¬ 
dently  the  problem  was  treated  by  Haugstad  (1977,  1978c, e).  by  Hubbard  etui 

(1978)  and  recently  also  by  Woo  et  al  I  1980).  An  interesting  heuristic  approach  to 

the  problem  was  taken  by  A  Young  (1976).  I  lie  problem  of  calculating  the  mean 
square  fluctuations  in  phase  and  bending  angle  does  not  appear  to  have  been  addres¬ 
sed  by  other  authors. 


4.2  An  atmospheric  model 

Explicit  evaluation  of  the  intensity  and  phase  scintillations  requires  a  model  to  be 
stipulated  for  both  u0  and  u, .  Ihe  former  quantity  has  already  been  assumed  to 
correspond  to  a  spherically  symmetric  and  exponential  atmosphere,  for  which 

i\,tr)  =  f'o ( r0  )  e\j>|  p’lr  r„  » |  (4.14) 

where  ,3' 1  =  II  is  the  scale  height,  and  r  r  r  is  measured  from  the  center  ol  the  planet. 
For  the  turbulence  component  we  assume 

lr)  =  <v\  i  r)> 1  :t'?fr)  (4.13a) 

<v]  ( r )>  1  2  -  <tq  (r0  ]>  1  *  exp|  d,(r  r„  )  |  (4.13b) 


1  hat  is,  u,  is  the  product  of  a  spherically  symmetric  component  that  varies  exponen¬ 
tially  with  height  on  a  scale  dj1  ■  Ht.  and  a  strictly  homogeneous  and  isotropic 
random  component  v*  of  zero  mean  and  unit  variance.  Ihe  power  spectrum  ■!’  of  f, 
will  be  assumed  to  be  of  the  Von  Kantian  type,  l  e. 


■l>(u)  -  ^v\ 


F(p;2i  I  f  exp(  V.f, u*  l 

7T3-Fl  S  I  1  •  I.*,u;  I1' 


(4.16) 


Here  u  =  luij.*-  u*-*-  u*)1  2  is  the  three-dimensional  sc  alar  wave  number.  I.„  and  f!„  are 
the  outer  and  inner  scale  of  the  turbulence,  respectively,  and  p  is  a  parameter  control¬ 
ling  the  slope  of  ‘1’  m  the  inertial  subrange,  where  I u  •'  V|,’  .  for  Kolmogorov 
turbulence  p  ~  I  1  /3. 
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For  purposes  for  numerical  evaluation  also  the  magnitude  of  <v]>  should  he  stipula¬ 
ted.  Turbulent  fluctuations  in  refractivity  may  be  related  to  those  in  temperature  and 
pressure  using  the  relation  (e  g,  Tatarskii  1971) 

v  =  7.9  •  10'5  p/T  (4.17) 


which  assumes  that  the  refractivity  from  any  vapors  present  can  be  neglected.  Small 
changes  in  refractivity,  V\  ,  are  thus  related  to  similar  changes  p,  in  pressure  and  I  , 
in  temperature  according  to 

v.  =  7.9  ■  KPl1-  ^  T,  1  (4.18) 

i  j  i :  1 

For  the  centimeter  to  kilometer  eddy  sizes  of  interest  here  the  temperature  relaxation 
is  slow  enough  to  suggest  an  adiabatic  relationship  between  pi  and  lj  (Lawrence  and 
Strohbchn  1970),  so  that  p,/p  =  7/(7  - 1  )•' T,  /  T,  where  7  is  the  ratio  of  specific  heats 
at  constant  pressure  and  volume.  With  this  assumption  equations  (4.1  7)  and  (4.18) 
combine  to  give 


<v\> 


(7-1) 


•  2  <’Ij_> 


(4.19) 


For  an  isothermal  ambient  atmosphere  the  fractional  mean  square  fluctuations  in 
refractivity  thus  vary  with  altitude  as  the  variance  of  the  temperature  fluctuations.  An 
estimate  of  the  latter  quantity  may  be  obtained  from  the  temperature  structure 
function 


I),  (Ar)  -  C\  lArl23  (4.20) 

valid  for  Kolmogorov  turbulence  and  separations  ^  I Ari  <  L0.  An  approximate 
value  of  <  J;  >  is  obtained  from  this  relation  bv  putting  Ar  =  l.„  and  noting  from  the 
definition  of  l.„  that  1)(  (1,0)  2  2  <  1  ]>.  Thus* 

1) !  (!•<,)  =  hi  Lr,1  '  2  <Tj  > 

The  fractional  fluctuations  in  refractivity  are  thus  approximately 


Willies  for  the  outer  scale  may  be  obtained  bv  direct  measurements,  or  its  value  may 
be  inferred  from  the  bulk  How  properties.  The  temperature  structure  constant,  (7 j- , 
may  be  estimated  from  the  mean  flow  characteristics  for  various  types  of  turbulence 
(see,  e  g,  Tatarskii  1971,  pp  72  74).  Alternatively  one  might  use  the  estimate 
2  2- 1  O' 1  K;  tn'2  1  corresponding  to  an  altitude  level  of  '  17)  km  in  the  earth’s  atmo¬ 
sphere  (button  et  al  1972).  Assuming  I.„  Id1  m  and  I’,,  -  200  K  we  obtain 

’  ‘"I 

This  estimate  tnav  be  appropriate  for  Jupiter.  For  the  atmosphere  of  Venus  the 
higher  ambient  temperature  (and  perhaps  also  the  smaller  value  of  L0 )  probably 
makes  equation  (4.22)  an  over-estimate  of  the  turbulence  strength  in  this  case.  As  we 
shall  see  m  set  (ion  7,  the  best  wav  to  estimate  <!']’>  is  probably  bv  direct  analy  sis  of 
radio  or  stellar  scintillation  measurements.  Lacking  such  direct  information  on  the 
turbulence  strength,  we  shall  explore  in  a  forthcoming  section  the  implications  of 
turbulence  of  the  relative  strength  given  bv  relation  (4.22). 


*  tor  Kolmogorov  nirlmli'ii  i-  tin-  more  ,  orrer  I  rrl.itionship  i\  given  bv  C.;|t.^'  l.')|<|  j  - 


(4.22) 


4.3  Radio  and  stellar  intensity  scintillations 


The  turbulence  model  described  in  the  previous  section  may  now  be  combined  with 
the  scintillation  formula  (4. 10)  to  yield  an  explicit  result  for  the  radio  scintillation 
index,  [-'.valuation  of  this  equation  by  transforming  to  polar  coordinates  yields  for 
3  <  p  <  6  (Haugstad  197 8c) 


7TI/2  (27T)3'I,/-  i„  ''T  |  /  Ho  tl)2 

'*■>  =  — — ^ - |scc  4  {~ ~ - 3\VR(^  0z;p)<c  >  (4.23a) 

^  JF.O  JK,0  K  '  '■ 


where 


P 


7T  —  -  1 

U'R((f>v.0z;p)  =/  (0vcos20  +  0zsi»2fl)“  d0 
-  it 


(4.23b) 


In  deriving  this  result  we  have  assumed  that  L0  >  a^j,  >  ay  v,  implying  that  the 
major  contribution  to  the  scintillations  is  from  turbulent  eddy  sizes  in  the  inertial 
subrange.  In  its  present  form  this  result  is  therefore  not  applicable  to  very  deep  occulta- 
tions,  where  a  [.  h  may  exceed  1.0-  We  shall  consider  this  case  separately  in  section  6. 

In  equation  (4.23)  coupling  between  the  turbulence  and  the  inhomogeneous  back¬ 
ground  appears  through  WR  only.  For  a  strictly  homogeneous  background,  or  when 
the  differential  ray  bending  is  sufficiently  small  that  0V  0Z  a  1.  WR  =  2ir  for  anv 
value  of  the  spectral  slope  parameter  p.  In  this  case,  after  manipulations  involving  the 
gamma  function,  equation  (4.23a)  may  be  shown  to  coincide  with  the  classic  Rvtov 
result  in  the  thin  screen  approximation. 


If  we  now  consider  occultations  for  which  the  approximation  0V  ?  1  can  be  made, 
but  with  no  restriction  on  0Z,  evaluation  of  equation  (4.23b)  gives 


WrU^p)  = 


2tr,  0Z  =  1 


■»P-l 


*> 

rip'  i : 


,  0Z *  1 


(4.24) 


Since  for  p  =  11/3  the  lower  va’uc  corresponding  to  0y<  1  is  about  0.6  times  the  upper 
value,  it  follows  that,  (i)  coupling  to  the  inhomogeneous  background  reduces  the 
scintillation  index  by  less  than  a  factor  of  two,  and  (ii)  this  reduction  is  essentially 
independent  of  0  s  0Z  whenever  defocusing  reduces  the  average  intensity  by  a  few  dB  or 
more.  I  his  result  disagrees  with  a  prediction  by  A  Young  (1976),  who  concluded  from  a 
heuristic  calculation  that  the  scintillation  index  should  increase  by  a  factor  0' 23  as 
deeper  atmospheric  regions  were  probed.  Although  neither  the  works  of  Hubbard  ct  al 
(1978)  or  Woo  et  al  (1980)  gave  explicit  results  for  the  coupling  effect,  their  expres¬ 
sions  for  the  scintillation  index  may  be  manipulated  to  yield  a  similar  result  to  that 
derived  here.  We  shall  see  in  the  next  section  that  the  rather  small  effect  of  the 
inhomogeneous  background  on  the  scintillation  index  does  result,  however,  in  sub¬ 
stantial  changes  in  the  shape  of  the  radio  and  stellar  scintillation  power  spectra. 


lor  stellar  occultations  the  point  source  intensity  correlation  function  should  be 
convolved  with  the  normalized  stellar  brightness  distribution  as  projected  above  the 
limb  of  the  planet.  In  the  spectral  domain  this  operation  yields  the  general  expression 
(c  f,  Hubbard  et  al  1978) 


//  ‘IMuJ  !0(UjJidUj_ 

0i 

for  the  scintillation  index  of  a  stellar  source.  In  this  expression  <I>^  (UjJ  is  the  spatial 
power  spectrum  of  <$j>  (i  c,  the  integrand  of  equation  (4.10)),  and 
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when-  ;'ir)  is  the  normalized  brightness  distribution  of  the  projected  stellar  disk.  For  a 
uniformh  illuminated  stellar  disk  of  projec  ted  radius  as  one  obtains 


IvMuJp 


2ii(asvV’ 

Vi 


where  J  ,  is  a  Bessel  function  of  order  one.  Instead  of  the  point  source  result  (4.10) 
we  thus  obtain 


0;>  -  HrrCk sin2  (I* 


asul 


(4.23) 


lbe  additional  factor  in  the  integrand  is  due  to  averaging  over  the  stellar  tlisk  and 
acts  to  suppress  contributions  to  the  scintillation  index  front  eddy  si/cs  smaller  than 
the  projec  ted  stellar  disk.  A  closelv  similar  result  to  ec|uation  (4.23)  has  been  obtained 
bv  1  lubbard  <7  <il  ( 19781. 


If  the  double  inequality  I.„  '  as  aj.  h  is  satisfied,  evaluation  of  equation  (4.25) 
yields  for  3  <.  p  <,  ft 


F  ( — *- )  .  -p 

4  2  ,«•«)  ^.\vs(<,v,0/;pK^>  (4.26a) 

F(l>  =  ) 

•> 


where  in  this  case 

U's(0v,<>/;p)  /  (0vcos :0  t  t>7sin 20)2d0 

■IT  ■ 


(4.26b) 


For  intermediate  and  small  values  of  the  ratio  a  /ay  ,  numerical  integration  of  equa¬ 
tion  (4.23)  is  required.  For  shallow  occultations  equation  (4.26b)  gives 


"s!  I- <>/.-!>) 


2tr,  6  y  1 

r:(3/2i 
1(3)  ’  / 


(4.27) 


Since  the  lower  value  for  0y  <•  1  is  "  1.76.  we  conclude  that  coupling  to  the  inhomo¬ 
geneous  background  reduces  the  stellar  scintillation  index  by  at  most  a  factor  ^  3-6 
under  the  previous  assumptions.  A  Young’s  prediction  for  this  case  was  an  increase  in 
the  scintillation  index  by  a  factor  ',  attending  a  compression  of  the  stellar  image 
by  a  factor  <)/  perpendicular  to  the  planetary  limb. 


It  is  instructive  to  express  the  scintillation  index  as  a  function  of  occultation  depth 
under  the  assumption  that  <fj>/r,o  is  constant.  To  this  end  we  first  note  that 


IF’  M2  <fj> 
2jtR!)2  e,2, 


We  have  here  used  the  fac  t  that  cv  -  (2jtR/1  l)1  2 v0  for  bending  angles  cv<(M/R)12 
and  M  «1)/II  as  before.  Fliminating  <nj>  between  this  relation  and  either  of  equa¬ 
tions  (4.23a)  or  (4.26a)  yields  for  ~  I 


•  > 


|K(l,),!iW)<  l>(ai)' 


V!  <•’!> 
R  i’U 


^Vwd.c^p) 


(4.28) 


For  radio  oc  c  ultations  W  \VRand  a  r  a  j.  0 .  while  for  stellar  occultations  W  =  \VS  and 
a  as.  1  he  numeric  al  factor  K  is  the  co-factor  of  equations  (4.23a)  and  (4.26a),  rrs- 
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Figure  4.1  Variation  with  occultation 
depth  f  of  the  </uantitiex 
FI  and  FI  for  radio  occul¬ 
ta  t  ions 


Figure  4.2  Variation  with  occultation 
depth  0' 1  of  the  quantities 
S2  and  ft  for  stellar  occul- 
t  at  ions 


1  he  solid  curves  display  £2  and  FI  when  coupling  to  the  inhomogeneous  background  is  properly 
accounted  for,  while  the  broken  curves  are  calculated  for  a  model  where  this  effect  is  not  included.  1  he 
curves  are  normalized  such  that  FI  in  the  latter  case  approaches  unity  as  0  1  -*• 


pectivcly.  If  the  spacecraft  (or  the  earth  for  stellar  occultations)  is  moving  perpendi¬ 
cular  to  the  incoming  ray  asymptote,  1)  is  approximately  constant  for  small  a,  and  all 
variation  of  <0:,>  with  depth  is  contained  in  the  two  last  factors  of  equation  (4.28). 
Figures  4.1  and  4.2  show  the  variation  with  occultation  depth  of  the  quantity 


ft  =  ,  1  0 \ *  W(  1,0:1 1/3) 
0  ‘lit 


(4.29) 


for  both  radio  and  stellar  occultations.  along  with  the  result  that  would  have  been 
obtained  if  coupling  to  the  inhomogeneous  background  were  falsely  neglected.  Note 
that  for  constant  1)  and  11,  M  is  proportional  to  atmospheric  density  or  pressure  at  the 
ray  periapsis.  Also  shown  in  the  figures  is  the  corresponding  variation  of  the  absolute 
mean  square  fluctuations  in  intensity,  which  are  proportional  to  FI  ~  f>:  Si.  In  con¬ 
trast  to  the  scintillation  index  which  increases  rapidly  with  depth,  the  absolute  inten- 
sity  scintillations,  <0J>  =  0‘  <0j>,  approaches  a  constant  determined  by  the  asymp¬ 
totic  value  of  W  as  0  -*■  0. 


We  conclude  this  section  by  proving  that  in  the  weak  scattering  regime  <$]  >  <  1  also 
when  coupling  to  the  inhomogeneous  background  is  accounted  for.  For  turbulence 
superimposed  on  a  homogeneous  background  this  result  has  previously  been  proved 
by  Salpeter  (19ti7)  and  others.  Consider  the  mean  square  fluctuations  in  optical  path 
at  exit  from  the  phase  screen.  This  quantity  may  be  obtained  from  equation  (4.1 1 ) 
by  putting  P  =  0  and  evaluating  the  double  integral  over  the  power  spectrum,  yielding 


<S?  (r,  )> 


l’(|,,)  )(|>  2) 


(4.30) 


Kliminating  <fj>  between  this  result  and  equation  (4.23a)  gives 

<$]>  =  K(p)  (i,FJii),-2"y?Y^/.;|»)  (ki  <  si(r,  )>| 

l-o  2tr 


(4.31) 


whore  K|[>I  (27T)''1’  "|sec  jrp'4i/|  4l'(p/2)|  is  of  order  of  magnitude  unity  for  p  -  11/3. 
Since  a  j.  „  1.0  per  assumy  ion,  and  we  have  shown  that\VR<  2tt,  it  follows  that 

1  in  the  weak  s  iltering  approximation  where  k<Sj(r,)>17  <  1.  Kvcn 
though  the  weak  scattering  condition  has  been  used  to  derive  the  previous  results,  it 
is  likely  that  only  the  weaker  condition  <<t>\  >  1  is  necessary,  however.  For  turbu¬ 

lence  superimposed  on  a  homogeneous  background  this  has  been  proved  by  Tatarskii 
(1971,  pp  23t>(.  One-dimensional,  numerical  computations  of  intensity  scintillations 
superimposed  on  an  inhomogeneous  background  indicate  that  <0j  >  <  1  is  indeed 
the  proper  criterion  also  in  this  case  (French  ct  al  1980). 

4  4  Temporal  power  spectra 

4.4.1  theoretical  phase  and  intensity  power  spectra 

I  he  scintillation  index  provides  only  an  averaged  representation  of  the  effect  of  the 
inhomogeneous  background  on  radio  and  optical  wave  propagation  in  turbulent  plane¬ 
tary  atmospheres.  More  detailed  information  is  contained  in  the  temporal  power 
spectrum,  the  shape  of  which  will  be  shown  here  to  depend  critically  on  both  occul- 
tation  depth  and  geometry  when  coupling  to  the  inhomogeneous  background  is 
accounted  for. 

1  he  phase  autocovariance  (unction  in  the  plane  \  -  1)  of  the  observer  is  given  bv 

B  =  k*  <S ,  (l),y  2 ,/;  )S ,  (l),y  2  4  Ay  o.z,  +  Az2)>  (4.32) 

Geometrical  optics  energy  conservation  implies  that  separations  Ay  2 ,  Az2  at  the 
receiver  project  into  ray  separations  Ay,  -  0vAy2  and  A/.,  =  0zAz2  at  ray  periapsis. 
Velocity  components  vv.  vz  of  the  spacecraft  are  therefore  related  to  corresponding 
components  v*  and  v*  at  the  ray  periapsis  according  to 

v*  =  0yvy,  v*  =  </>zvz  (4.33) 

It  is  customary  to  adopt  the  ’’frozen-field"  hypothesis,  whereby  turbulent  inhomo¬ 
geneities  are  assumed  to  be  converted  across  the  line  of  sight  at  constant  velocity, 
which  for  the'  present  case  equals  the  relative  velocity  between  that  of  the  rav  at 
periapsis  and  the  atmospheric  bulk  velocity  ,  due  to  either  planetary  rotation  or  the 
presence  of  atmospheric  winds,  or  both.  We  shall  assume  that  the  components  of  the 
latter  are  small  compared  to  v*.  v*.  In  this  case  ray  separations  Ay,,  Az,  at  the 
periapsis  are  related  to  time  separations  r  according  to  Ay,  =  v*T,  Az,  =  v*t.  Using 
relations  (4.33)  and  drawing  on  tire  analysis  of  section  4.1,  we  find  that 

B  -  k'  ff  <*’s|  (Uj_)  exp|i(uy0yvy  4  U/0zvz)r|duj_  (4.34a) 

where 

<I’S|  (uxl  -  27rUl*(ui)  cos:(1»)  (4.34b) 

is  the  spatial  power  spectrum  of  the  scintillations  in  phase  path,  as  given  implicitly  bv 
equation  (4.1  1 ).  Fourier  transformation  over  r  yields  for  the  temporal  power 
spectrum  of  the  phase  scintillations 

W,  (w)  -  2trk‘  ff  <I'S|  (ux)  6|uv0yv  f  u z</>zvz  to|dux  (4.33) 

Using  the  turbulence  power  spec  trum  (4.1(i)  this  integral  may  be  evaluated  by  trans¬ 
forming  to  polar  coordinates  and  performing  the  integration  over  the  angular  variable. 


After  a  suitable  change  of  integration  variable  in  the  remaining  integral  we  find  that 
fllaugstad  1979a) 

■>  ^  xeos:|(  ):  |0yC°s2  (<F  +  6)  4  0/siir  (dj  +  6i|j 

U’|  (w)  '  « ,  uj  11  f  - - — -  -  — - d\  (4.36) 

1  1 1  |l  m‘jx  )2\''r-(*:  1),: 

CJ„ 

We  have  here  defined  0t  -  0,  =■  cos''(l/x),  to,,  =  toc  =  2JT1”  v*/a|,- q  , 

6  -  tan'1  (v^Vy),  and  k,  -  | 4jt '  -  P( p/2 l/l'f < |»-:4)/2 ) K’k2 1. ?,<!'?>/%•*.  We  observe  from 
this  result  that  the  shape  of  the  phase  power  spectrum  is  controlled  by  the  four 
parameters  cou,  uJ(.,  0V  and  0Z.  In  view  of  their  definition,  u>0  and  c j(.  are  the 
characteristic  frequencies  associated  with  the  ra\  speed  at  periapsis  and  the  outer 
scale  and  the  free  space  Fresnel  /.one,  respectivelv. 

I hrough  entirely  analogous  procedures  the  radio  and  stellar  intensitv  scintillation 
power  spectra  are  obtained  as 

\sin:  {(  ):  lO^.cos2  (dj+6 )  +  dosin’  ( 0j  t  6)|j> 

W :  ( oj )  -  K:^>  £  7  W<7  -  -  ’  cK  i  4.S7) 

1  [if  n>: 


\  ’  |0v  a  i)s:  (d:  •  6)  *  0  sin‘(d:  *  6)1" 

W,(w)  =  11  f  -  -  —  |i(tJX)d\  (1.38) 

1,1  |  1  +  ):  l)1’  uc 

W.» 

respectively.  For  the  former  result  (  16ir"  l'(p/2)/I  ((p-A)/  — ))  fek"  I <1/7 >/\ **' ;  0,, 

02  and  oa„.  0JC  are  defined  as  before.  Flic  stellar  intensity  scintillation  power  spect¬ 
rum  is  subjected  to  the  same  limitations  as  the  expression  (1  .lib)  for  the  scintillation 
index.  In  this  c  ase  u,  v*'as,  and  k  ,  (  1  <»7T  ’  I ' ( p / 12 ) / 1 ' ( ( p - .‘i ) / 12 ) )  I):1,;J  <tq  >/(asv*4  I. 

We  note  that  all  of  the-  power  spec  tra  calculated  here  are  valid  onlv  for  frecpiencies 
much  larger  than  the  characteristic  frecjucncy  v"7II  associated  with  either  of  the  scale 
heights  II  or  1 1( . 

Figures  4.3  to  4.9  show  the  phase  and  intensity  power  spectra  as  a  function  of 
normalized  angular  frequency  w/u>(  or  oj/co,,.  for  different  values  of  m  v*/v*  and  0. 
assuming  0V  1.  In  all  cases  computed  we  have  assumed  u)r/u)„  -  3.  For  the  phase 
spectra  coupling  to  the  inhomogeneous  bac  kground  is  seen  to  reduce  the  fringes  in 
the  spectrum  for  to  >  oJ(  in  a  central  occultation  Int  0)  as  deeper  atmospheric 
regions  are  probed.  Figure  4. displays  the  change  in  spectral  shape  from  a  central  to 
a  more  nearly  grazing  occultation,  for  a  fixed  value  0  0.1  of  the  mean  signal  inten¬ 

sity. 

by  contrast  to  this  behavior,  the  intensity  power  spectra  display  a  strong  dependence 
on  both  occultation  geometry  (tn)  and  depth  (0).  According  to  Figure  1.6.  Un¬ 
characteristic  fringes  in  the  power  spectrum  above  the  Fresnel  frequency  u>(.  alto¬ 
gether  disappear  as  deeper  regions  are  probed  in  .1  central  occultation.  In  a  grazing 
occultation.  however,  these  fringes  become  increasingly  deeper  and  follow  a  characte¬ 
ristic  decrease  in  scintillation  power  below  u)(.  as  denser  regions  are  probed, 
f  igure  4.8  shows  the  progressive  change  in  spectral  shape  from  .1  c  entral  to  a  more 
nearly  grazing  occultation  for  a  fixed  occultation  depth  corresponding  to  0  0.1. 

With  only  minor  modifications  the  changes  in  spec  tral  shape  with  m  and  0  for  a  point 
source  are  reproduced  also  for  an  extended  source  when  the  projected  stellar  radius 
substantially  exceeds  the  free-spare  Fresnel  zone.  Figure  4.9  shows  the  change  in 
spectral  shape  with  occultation  depth  for  a  grazing  occ  ultation.  If  the  projec  ted  stellar 
radius  is  small  compared  to  the-  free-spare  Fresnel  zone,  averaging  over  the  stellar  disk 
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Figure  4.3  Phase  scintillation  power  spec¬ 
trum  at  various  depths  0  in  a 
central  oeeultation.  assuming 


Figure  4.4  Phase  scintillation  power  spec¬ 
trum  at  various  depths  <p  tn  a 
grazing  oeeultation,  assuming 
ojk./co0=3 
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Figure  4.5  Phase  scintillation  power  spec¬ 
trum  for  different  occupation 
geometries  and  a  fixed  depth 
0  0. 1 .  assuming  oj(./ci7()  "  J 


Figure  4.6  Point  source  intensity  scintilla¬ 
tion  power  spectrum  at  vari¬ 
ous  depths  0  in  a  central  oc¬ 
cupation,  assuming  oj^./oJq  -  3 
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will  yield  no  substantial  effect  on  either 
the  scintillation  index  or  the  power 
spectrum.  In  this  case  power  spectra 
will  correspond  closely  to  the  point 
source  scintillation  spectra  of  Fig¬ 
ures  4.6  -4.8. 


0  =  0.33 


The  classical  Rytov  results  for  the 
phase  and  intensity  power  spectra  cal¬ 
culated  for  turbulence  superimposed  on 
a  homogeneous  background  correspond 
to  the  cases  0  =  1.0  in  Figures  4.3,  4.6 
and  4.9. 


6  =  0.1 


0  =  0.01 


W3(w) 


Figure  4.9  Stellar  scintillation  power  spec¬ 
trum  for  various  occupation 
depths  0  in  a  grazing  occupa¬ 
tion,  assuming  co(/u>0=  3  and 

aJah  .0  *  7 


4.4.2  Observed  power  spectra 

Power  spectra  calculated  from  the  Mari¬ 
ner  5  and  10  Venus  encounters  (Woo 
et  al  1974,  Woo  1975a,  Woo  et  al 
1980|  and  from  the  Pioneer-Venus 
entry  probes  (Woo  et  al  1979)  have 
been  reported  in  the  literature.  For  mis¬ 
sions  by  other  planetary  spacecraft  and 
for  stellar  occultations  no  correspond¬ 
ing  results  have  been  published. 

Figures  4.10a-b,  taken  from  Woo  et  al 
(1974),  show  the  recorded  part  of  the 
Mariner  5  S-band  entrance  and  exit 
occultation  as  a  function  of  universal 
time  and  altitude  above  the  mean 
Venusian  surface.  Regions  marked  A  to 
1)  were  suspected  of  containing  struc¬ 
tures  due  to  turbulence,  and  were 
selected  for  further  analysis.  After  re¬ 
moval  of  a  quadratic  trend  function 
from  the  exit  data.  Woo  et  al  (1974) 
calculated  the  temporal  power  spectra 
and  compared  with  predictions  from  a 
layered  turbulence  model  fitted  to  stan¬ 
dard  weak  scintillation  formulas, 
assuming  a  priori  that  F0  <  aj.-  ().  In  a 
later  comparison  Woo  ( 1  975a)  claimed 
improved  theoretical  fit  to  the  Mari¬ 
ner  5  spectra  by  instead  assuming 
L0  >  a j.  () .  Neither  analyses  accounted 
for  coupling  to  the  inhomogeneous 
background,  however.  As  apparent 
from  the  previous  analysis,  accounting 
for  this  effect  is  instrumental  in  the 
present  case  where  0'1  may  reach  a 
hundred  or  more. 
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Figure  4. 1  On 

Time  history  of  signal  intensity  for  the 
recorded  part  of  t he  Mariner  5  I'enus  ent¬ 
rance  occultation,  as  adapted  from  Woo  et 
all  1974) 


Figu  re  4. 1  Oh 

lime  history  of  signal  intensity  for  the 
recorded  part  of  the  Mariner  5  Venus  exit 
occultation,  as  adapted  from  Woo  et  al 
(1914) 


The  upper  abscissa  uives  the  altitude  <>1  the  ray  periapsis  above  the  mean  Venusian  surface.  Regions  A 
to  l)  were  suspected  for  their  turbulence  content  and  selected  for  further  examination. 


In  Figures  4.1  1  and  4.12  we  have  overlaid  experimental  spectra  from  regions  A  and  C 
of  the  Mariner  5  exit  occultation  with  theoretical  spectra  calculated  from  equa¬ 
tion  (4217).  (Kxperimcntal  spectra  from  other  atmospheric  regions  do  not  raise  signifi¬ 
cantly  above  the  noise  level.)  The  relevant  input  parameters  for  this  calculation  were 
obtained  from  Fjcldbo  el  al  (1971),  yielding  aj.  ()  =  1,3-1  ()3  m,  vz  =  7.5-1 03  ms'1  and 
Vy  "  3- 1  ()2  ins'1  for  both  regions  A  and  C.  Since  0V  ?  1.  we  have  v*  s  vv  and 
v*  =  0V,,,  with  0  estimated  from  Figure  4.10b  as  s0.05  for  region  A  and' =0.01  for 
region  C.  Weak  scintillation  theory  is  here  appropriate  since  the  computed  scintillation 
index  for  both  events  is  small  compared  to  unity. 

It  is  apparent  from  Figure  4.1  1  that  the  position  of  the  "knee"  in  the  theoretical  and 
experimental  power  spectra  agree  quite  well,  a  fact  strongly  favoring  turbulence  as  the 
cause  ot  the  observed  fluctuations  in  signal  strength.  For  region  A  there  is  a  clear 
discrepancy  between  the  spectral  shapes,  however.  While  the  theoretical  spectrum 
appears  to  match  the  average  high  frequency  slope  of  the  Mariner  5  spectrum,  the 
strong  fringes  in  the  experimental  spectrum  are  not  reproduced  in  the  theoretical 
prediction.  F.ven  though  both  spectra  display  a  decrease  in  scintillation  power  below 
the  Fresnel  frequency,  they  are  clearly  discrepant  also  in  this  region.  The  reason  for 
this  difference  is  presently  not  clear;  a  contributing  factor  may  be  a  possible  low 
quality  of  the  experimental  spectrum  itself,  however.  Indeed,  the  remarkable  uni¬ 
formity  of  this  spectrum,  as  compared  to  the  other  Mariner  !>  spectrum  in  Figure  4.12 
and  when  contrasted  with  the  90%  confidence  limits  indicated  in  the  figure,  clearly 
cast  doubts  on  the  validity  of  the  spectrum.  It  may  be  that  the  "zero-fill"  operation 
applied  to  the  original  data  in  order  to  evaluate  the  Fourier  transforms  at  closer¬ 
spaced  intervals  than  would  otherwise  lie  possible  (Woo  et  al  1974,  pp  1 700),  has 
introduced  substantial  correlation  between  neighboring  power  estimates  with  a  corres¬ 
ponding  reduction  in  power  variance.  For  a  sharp  spec  tral  knee  this  correlation  would 
leave  essentially  unchanged  the  position  of  the  knee,  but  would  tend  to  smooth  any 
structures,  eg  fringes,  in  the  spectrum.  This  does  not,  however,  explain  the  regular 
fringes  that  still  persist  in  the  experimental  spectrum,  nor  its  strong  decay  at  sub 
Fresnel  frequencies. 
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Figure  4.1 1  Log-amplitude  power  spect¬ 
rum  for  region  A  of  the  Mari¬ 
ner  5  exit  occupation  (dotted 
curve),  as  adapted  from  Woo 
et  al  (1974) 

Superimposed  is  the  theoretical  spectrum  (solid 
curve)  based  on  the  actual  occupation  parameters, 
and  displaced  along  the  ordinate  axis  to  vteld 
maximum  fit.  The  error  bar  represents  a  VO"!,  i  on 
fidence  interval. 


Figure  4.12  Log-amplitude  poieer  spect¬ 
rum  for  region  ('.  of  the  Mari¬ 
ner  5  exit  occupation  (dotted 
curve),  as  adapted  from  Woo 
et  al  ( 1974) 

Superimposed  is  the  theoretical  spectrum  (solid 
curve)  based  on  the  actual  occupation  parameters, 
and  displaced  along  the  ordinate  axis  to  yield 
maximum  fit.  I'he  error  bar  represents  a  00"!.  con 
fidenc  e  interval. 


The  experimental  spectrum  shown  in 
Figure  4. 1  2  displays  scatter  among  indi¬ 
vidual  spectral  estimates  that  appears 
more  consistent  with  the  indicated  90% 
confidence  limits.  There  is  an  indica¬ 
tion  that  both  the  predicted  position  of 
the  spectral  knee  and  a  decrease  in 
power  at  lower  frequencies  is  repro¬ 
duced  in  this  event. 

Power  spectra  calculated  by  Woo 
(1973a)  from  the  Mariner  10  Venus 
occupation  lack  the  low-frequency  part 
below  to  ,  and  also  appear  to  be 
smoothed  at  higher  frequencies  to  an 
extent  that  would  make  a  comparison 
with  the  theoretical  power  spectra  vir¬ 
tually  impossible.  Due  to  the  much 
smaller  turbulent  medium-to-probe  dis¬ 
tances  involved,  the  Pioneer  Venus 
entry  probes  failed  to  detect  any  inten¬ 
sity  scintillations,  establishing  essen¬ 
tially  an  upper  bound  on  the  refracti- 
vitv  fluctuations  in  this  case  (Wo o  et  al 
1979). 


Figure  4. 13  shows  a  power  spectrum 
based  on  the  Mariner  10  Venus  occupa¬ 
tion  profile  as  calculated  by  l.ipa  and 
Tyler  (1979).  The  spectrum  corres¬ 
ponds  in  altitude  approximately  to  the 
region  labeled  A  by  Woo  et  al.  The 
figure  shows  both  the  experimental  and 
theoretical  power  spectra,  together 
with  part  of  the  exit  occupation  inten¬ 
sity  time  profile  and  the  small  segment 
subjected  to  spectral  analysts.  The 
experimental  spectrum  is  calculated  by 
first  fitting  a  polynomial  of  varying 
degree  to  the  intensity  background,  and 
then  subtracting  this  trend  function 
from  the  signal  to  obtain  the  fluctuat¬ 
ing  part.  The  residual  variance  of  the 
signal  relative  to  the  polynomial  trend 
is  shown  as  a  function  of  the  degree  N 
of  the  polynomial  in  Figure  4. 1 4.  Ini¬ 
tially  the  residual  variance  is  seen  to 
decrease  rapidly  with  X,  but  then  be¬ 
coming  virtually  insensitive  to  this  para¬ 
meter  for  N  >  3.  Thus  a  fifth  degree 
polynomial  was  subtracted  from  the 
intensity  time  profile,  and  the  differ¬ 
ence  was  normalized  by  the  trend  func¬ 
tion  and  subjected  to  a  power  spectral 
analysis,  the  result  of  which  is  shown  in 
Figure  4.13.  The  corresponding  theore¬ 
tical  spectrum  has  been  shifted  along 
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Figure  I.  /.>  Experimental  and  theoretical  power 
speetrum  corresponding  to  the 
iipproxmuitelx  Id  s  portion  of  the 
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the  ordinate  axis  to  yield  maxi 
mum  fit-  The  low-frequency  parts 
are  seen  to  agree  fairly  well.  Also 
the  high  frequency  portions  are 
not  in  obvious  disagreement,  the 
large  scatter  of  the  data  points 
taken  into  account.  The  possible 
discrepancy  between  the  position 
of  the  Fresnel  knees  may,  at  least 
in  part,  he  attributed  to  the  limi¬ 
ted  statistical  significance  of  the 
experimental  spectrum. 

Preliminary  analysis  of  other  parts 
of  the  Mariner  10  oeeultation  pro¬ 
file  suspected  of  containing  struc¬ 
tures  due  to  turbulence  has  not 
yielded  a  consistent  result.  Most 
of  the  fluctuations  seen  in  the 
deeper  parts  of  the  oeeultation 
profile  turn  out  to  be  noise,  con¬ 
sistent  with  the  findings  of  Woo 
and  colleagues  for  the  Mariner  5 
oeeultation.  At  least  in  one  other 
region  docs  the  power  spectrum 
raise  above  the  background  noise, 
but  yet  does  not  match  its  theore¬ 
tical  prediction.  It  is  possible  that 
this  region  does  not  correspond  to 
well-developed  turbulence,  or  that 
the  statistical  significance  is 
simply  too  low  to  provide  a  mean¬ 
ingful  comparison  with  the  statis¬ 
tical  mean  spectrum  of  the  theore¬ 
tical  analysis. 


.  ,  ,,  ,  i  ■  ,  Spectra  of  higher  statistical  stgnift- 

inilieateil  in  tin •  ton  er  h  it  part  o,  1  .  ,  . 

,  1  cance.  essential  to  an  experimental 

the  i<’ ure  r  ,  ,  , 

assessment  ol  the  theory  deve¬ 
loped  here,  may  in  principle  be 
obtained  from  the  averaging  of 
power  spectra  from  different  occultations,  or  through  substantial  spatial  averaging  in  a 
single  grazing  oeeultation,  where  the  ray  moves  many  Fresnel  scales  at  essentially 
contain  altitude.  Removal  of  the  background  intensity  from  the  scintillations  poses 
yet  another  fundamental  problem.  When  the  scintillation  spectrum  and  the  trend 
function  have  finite  power  in  a  common  frequency  range,  which  may  normally  be  the 
case,  complete  separation  of  fluctuations  and  trend  is  not  possible  even  in  principle 
(Yaglom  l‘H»2).  It  is  nevertheless  imperative  that  a  common  and  well-defined  proce¬ 
dure  be  used  when  removal  of  the  trend  is  attempted.  The  technique  suggested  here  is 
both  well-defined  and  statistically  sound  (  Draper  and  Smith  1  960.  Chatterjee  and  Price 
71. 


Figure  4. 14  Rest  variance,  (Jp  ,  o\  intensity  fhu  Haitians  around  a  po/x nainial  trend,  as 
function  of  tin  decree  X  o'  the  txdx  nomiui 
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4.5  I  hr  weak  scintillation  approximation  in  planetars  decollations 

Consistent  with  the  contention  made  in  \e<  ti«m  4  b  we  shall  assume  here  that  tin 
theoretical  results  of  the  previous  sections  are  valid  it  ■  1.  a  t  csp.  e  t  tv  e  ot  the 

mannitude  of  the  phase  flucttiatiotis  at  exit  I  root  the  a  t  m.  .heri  below  wc  explore 
to  what  extent  this  weak  scintillation  assumption  m  av  h.oe  :« e$  satis) ted  in  past 
radio  and  stellar  occullations,  examining  both  theoieteal  a  i.t  experimental  evidence. 

He  define  formally  a  critical  ociuliatioit  depth  a",  mu  I.  that 

<0\  tO*  !>  =  1 

1  he  theoretical  results  developed  here  are  thus  valid  onlv  tot  oecultatioti  depths  such 
that  0  0*.  In  terms  of  0* .  the  scintillation  index  at  an  crbttrarv  depth  0  is  from 

equations  ( 4.24  )  and  (4.26l 

pm  Htl,0;i)i  0*  *  1  0  ,  , 

<0i>  -  -  -  J  '  I  i  -  -  i  (4.49 1 

\\(  l,0*;p)  0  1  o* 

w  here  \V  takes  either  of  the  forms  aiven  bv  eiptattons  1 4.24b  i  and  I  4.2(ib)  for  radio 
and  stellar  o<  cultatioiis,  respectively.  In  table  4.1  we  have  calculated  the  critical 
depth  0*  for  two  past  radio  and  stellar  occullations.  lor  all  cases  we  have  assumed 
-  10  s  and  I.„  -  H/4.  where  H  alone;  with  the  free-spare  f  resnel  /.one  and 
the  stellar  annular  radius  are  ejvcn  in  I  able  4.1.  Since  0*  ~  1  for  all  four  occultations. 
ec|uation  14.49)  simplifies  to 


10  ' 0 ) " 


(4.40) 


when  0  1.  It  is  apparent  from  (able  4.1  that  the  weak  scintillation  condition  may 

have  been  violated  in  the  deeper  parts  of  these  occ  libations  if  turbulence  of  the 
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strength  found  at  the  ~15  km  level 
of  the  earth’s  troposphere  is  also  pre¬ 
sent  at  the  higher  altitudes  probed 
on  Venus,  Jupiter  and  Neptune. 

Data  from  the  Mariner  5  Venus 
occ'ultation  have  been  analyzed  by 
Gurvich  (1969),  Golitsyn  and  Gur- 
vich  (1971)  and  by  Woo  vt  al  (1974) 
for  the  variance  of  the  log-amplitude 
fluctuations.  Accounting  fora  factor 
of  4  between  this  quantity  and  the 


fable  4.1  Approximate  critical  occultation 
depths  <p*  and  minimum  occul¬ 
ted  intensity  $mtn  Jor  two  past 
ratlin  and  stellar  occultations 
for  thi*  Ntfllar  intubations  bv  Jupiter  and  Nep¬ 
tune  art*  from  l.illcr  fl  til  (1974)  and  Freeman  and 
l.yni^a  i)970),  respectively . 


scintillation  index,  Golitsyn  and  Cur- 
vich  (1971)  estimated  the  scintilla¬ 
tion  index  to  be  about  2-10"2  for 


both  the  night  side  immersion  and 
day  side  emersion.  This  value  is  an 
average  over  altitude  corresponding 
to  an  approximate  pressure  range  of 
1  -5  bar.  Woo  et  al  (1974)  obtain  a 
similar  figure  (— 4-10'2)  for  the  scin¬ 
tillation  index  in  the  -60  km  altitude  region  of  pronounced  turbulence  identified  in 
their  analysis.  In  his  at  tlx  sis  of  turbulence  in  the  Mariner  10  Venus  occultation,  Woo 
(1975a)  found  a  value  =0.1  for  the  scintillation  index  in  the  same  —55  — 70  km 
altitude  region  previously  probed  by  Mariner  5.  While  the  radio  occultation  measure¬ 
ments.  including  that  of  Vetter  a  9  (Timofeeva  et  al  1978).  appear  mutually  consistent, 
estimates  of  the  refractive  index  structure  constant  that  would  yield  much  higher 
values  for  the  scintillation  index  have  been  derived  from  radio  link  analyses  of  the 
Venera  1  8  landers  t Kolosov  et  al  1970,  Yakovlev  et  al  1974).  While  this  discrepancy 
lias  not  yet  been  sutislat lords  resolved,  the  recent  Pioneer  Veins  entry  probe 
measurements  appear  to  support  the  radio  occultation  data  (Woo  >‘l  ’  1979).  lhe 
spacecraft  measurements  thus  appear  to  predict  a  lower  turbulence  -irength  than 
would  be  predicted  from  extrapolation  of  measurements  in  the  earth’s  troposphere. 


No  determination  of  a  scintillation  index  for  the  Pioneer  Jupiter  missions,  for  the 
Voyager  Jupiter  and  Saturn  missions,  or  for  any  of  the  several  past  stellar  occulta¬ 
tions  have  been  reported.  The  large  excursions  in  intensity  commonly  observed  in 
stellar  occultations  strongly  indicate,  however,  that  non-uniform  atmospheric  struc¬ 
tures.  possibly  turbulence,  would  yield  a  formal  value  for  the  scintillation  index 
comparable  to  or  larger  than  units.  It  is  deemed  important,  after  properly  separating 
the  scintillations  from  the  main  structure  of  the  event,  to  calculate  the  scintillation 
index  also  in  such  cases.  Since  numerical  simulations  indicate  that  the  scintillation 
index  saturates  to  unity  also  when  an  inhomogeneous  atmospheric  background  is 
accounted  for  (French  et  al  1980)’*,  values  of  <$]>  that  were  larger  than  unity,  for 
example,  would  immediately  rule  out  turbulence  as  the  sole  cause  of  these  intensity 
variations.  Additona!  clues  as  to  the  nature  of  the  underlying  atmospheric  structure 
might  then  be  obtained  from  an  examination  of  the  scintillation  power  spectrum,  as 
discussed  more  fuiiv  in  set  t  ion  7. 


*  Kci  ent  results  In  this  author,  based  on  ,i  slight  i>cncrah/uttou  ol  the  st.il imuirv  |diuse  iethnit|ue 
developed  in  seitionti.  shows  that  the  seimillation  index  indeed  saturates  to  unity  also  in  the 
present  e  til  an  inhouio'aeneoiis  hat  Aground. 
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5  PROFILE  ERRORS  DUE  TO  TURBULENCE 

As  previously  stated,  the  inversion  procedure  to  obtain  profiles  of  refractivity,  tempe¬ 
rature  and  pressure  from  occultation  measurements  assumes  the  atmospheric  refracti¬ 
vity  to  have  a  spherical  distribution.  Turbulence  violates  this  condition,  and  since  the 
perturbations  in  signal  intensity  and  Doppler  frequency  due  to  turbulence  arc  non¬ 
linear  functionals  of  the  generic  refractivity  fluctuations,  such  profiles  will  contain 
both  random  (zero  mean)  and  systematic  (finite  mean)  errors.  We  estimate  in  this 
section  the  magnitude  of  these  errors,  and  shall  find  that  the  random  and  systematic 
errors  are  both  small  in  the  limit  of  weak  scintillations,  consistent  with  the  simple 
argument  given  in  section  2.2. 


5.1  Random  errors  in  atmospheric  profiles  from  radio  or  stellar  intensity  measurements 

Let  nt  =  1  +  v*  be  the  refractive  index  calculated  from  the  formal  inversion  given  by 
equation  (2. 1  b).  Let  a  ~  a+  p,  where  p/a  '  p/R  -•  1.  lo  lowest  order  in  v  this  equa¬ 
tion  becomes 

v* 1  ro  >  =  4 *  f  M(a  *  pip1  :  dp  (5.1) 

TTl)  i\  p 

We  have  here  used  the  approximation  cosh"1  (a’/a)  '  ( 2p.  R ) 1  ;,  and  the  argument  a*  p 
of  M  =  cvl)/I I  signifies  that  o  is  the  bending  angle  of  a  ray  with  impact  parameter  a*  p. 
Since  (2trR/lI)'  ’('0  for  small  bending  angles  in  a  sphericalK  symmetric  and  iso¬ 
thermal  atmosphere,  direct  substitution  in  equation  (5.1  I  yields  v  i’(l,  as  re<|uired  in 
this  case. 


When  turbulence  is  present,  the  signal  intensitv  may  be  expanded  in  powers  of  t',  as  in 


0  0Z  *0,  I  ...  0Z(1  t  <5,  *•..)"-  (1  *  M)  1 


(5.2) 


Ihe  last  equality  assumes  that  the  signal  intensity  is  still  related  to  the  ray  bending 
angle  as  if  no  turbulence  were  present  in  a  shallow  occultation.  Combination  of  the 
two  previous  equations  yields  for  the  lowest  order  error  i’*. ,  relative  to  the  "true” 
refractivity  v0  of  the  ambient  atmosphere 


^  -  1  .  -  /  v  *'u‘|,VMp 

Ttl)  R  „  0z(a  +  |>) 


(3.3) 


Ihe  argument  a+ p  of  0,  and  </>z  signifies  again  that  these  quantities  are  calculated 
with  respect  to  a  ray  of  impart  parameter  a  t  j>.  Clearly,  as  <$j  >  -*  0  the  refractivity 
error  vanishes  and  e  tq, . 


Since  is  a  stochastic  quantity,  only  its  mean  and  variance  (and  possibly  also 

higher  statistical  moments)  are  of  interest.  I  he  mean  value  is  zero,  since  clearly 
<0i  >  0.  The  fractional  variance  of  this  quantity  is,  from  equation  (5.5) 


2  <0i  > 

n:  Rl):iv 


//'(  1  •  Mc^1’1 
0  0 


1(1  '  Mt 


(.-dt(p, 


"Pi 


:1>‘  '/“(P: 


Pi  M|’i 


(5.4) 


All  atmospheric  variables  where  ,tti  argument  has  not  been  assigned  tire  referenced  to 
the  rav  impact  parameter  a'  ,e  In  establishing  this  equation  we  have  used  the  fact 
that  <0,  la  +  p,  )0,  (a+p; )  >  <0]  (a)  >exp(-/jt(p,  t-p2))p(p2  p,  ),  where  p  is  the  (nor¬ 

malized)  autocorrelation  function  of  the  (normalized)  intensity  scintillations  at  the 
receiver;  (L1  II  and  |3t'  II,  as  before.  In  the  following  we  shall  not  distinguish 
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between  H  and  1  It .  Let 

I  =  /p(Ap)d(Ap)  (5.5) 

n 

be  the  correlation  length  of  the  intensity  scintillations,  and  let  I-' ( p  |  ,p 2 )  be  a  function 
that  changes  slowly  on  a  length  scale  M,  so  that 

7/ b(p,  ,p3)p(p2  Pi )  dp,  d|>2  5  I  /  F(p, ,p,  )dp,  (5.6) 

on  o 

It  may  be  shown  using  the  definition  (5.5)  together  with  equation  (4.7)  that  I  is  typi¬ 
cally  on  the  order  of  a^.  ()  in  a  shallow  occultation.  Therefore,  in  the  range  0<(Pi. 
P:>  where  e(|iiation  (5.4)  derises  its  principal  contribution,  the  co-factor  of  p  in 
the  integrand  of  this  cipi.it ion  satisfies  the  condition  for  using  relation  (5.6).  For  the 
0  1  region  of  principal  interest  we  thus  obtain 


from  equations  (2.61  and  (2.7)  it  may  be  shown  that  the  corresponding  fractional 
errors  in  temperature  and  pressure  are  approximately 


**5  >  <P*i 


1  <0  > 

11  1 


where  I*,  and  p*,  are  the  lowest -order  errors  in  temperature  and  pressure,  respecti¬ 
vely.  relative  to  the  "true”  ambient  values  T0  and  p0.  Thus  the  fractional  mean 
square  errors  in  derived  refractivity,  temperature  and  pressure  are  all  small  in  the  limit 
ot  weak  scintillations.  For  stellar  occultations  the  same  result  applies,  except  that  the 
correlation  length  1  may  be  different  from  s  ap  ,,  in  this  case.  If  the  projected  stellar 
radius  as  significantly  exceeds  a j.  „ ,  I  =  as>  and  relations  (5.7)  and  (5.8)  are  still  valid 
if  also  as  -v  11. 

It  is  interesting  to  note  that  if  the  Baum  and  Code  equation  (2.5)  were  used  to 
estimate  the  scale  height  and  hence  the  temperature,  the  corresponding  mean  square 
fractional  error  in  temperature  would  be  approximately  <$]>  for  0  <1.  Thus  the 
integral  inversion  technique  is  less  sensitive  to  intensity  scintillations  than  curve-fitting 
to  the  Baum  and  Code  equation  by  a  factor  1/11  <  1. 


5.2  Random  errors  in  atmospheric  profiles  from  Doppler  measurements 

For  a  spherically  stratified  refractivity  field  the  atmospheric  Doppler  frequency  f  is 
related  to  the  small  bending  angle  a  approximately  as  in 


where  \.  is  the  component  ot  the  spacecraft  velocity  in  the  plane  of  refraction  that  is 
perpendicular  to  the  incoming  ray  asymptote.  When  turbulence  is  present,  the 
Doppler  frequency  may  be  expanded  as  f=  fn  +  f, +...,  where  f0  =  ©(i>0),  fi  =  @(v\  ) 
etc.  Since  equation  (5.9)  is  still  used  to  infer  the  bending  angle  also  when  turbulence 
is  present,  the  estimated  angle,  or  -  a  ,  may  be  similarly  expanded  in  powers  of  tq 
according  to  oq  a*o  a*t  +  •  •  •  From  equation  (2.  lb),  and  using  a  similar  expansion 
of  cosh  '(a  /a)  as  in  equation  (5.1),  we  find  that  the  estimated  refractivity  v  is 


n  K  o  dj) 


-(L-),  7<|f1(a,!>)  p'Mp 


R 


(5.11) 


;j  r- 


44 


whore,  dearly,  -  Q){v0),  v* 1  =  ©(^i),  etc.  Using  the  standard,  small  angle  expres¬ 
sion  for  a,  equation  (5.10)  is  readily  shown  to  yield  c*,,  =  v0.  The  lowest-order  error 
to  this  ’’true”  refractivity  profile  is  thus  represented  by  equation  (5.11).  In  terms  of 
the  covariance  function  (assuming  local  homogeneity) 


B(Pi  4>2 ) 


<df,(a+p2_)df1(aljll)>  =  <(df1(a))>^  p(p2  _pi  } 
dp  2  dp  |  da 


the  mean  square  fractional  error  in  refractivitv  is  in  this  case 


X2 


-  -A  ,  <(—)>  (J  e 


,-^tlPi  +  Pa),,,. 


7Tfc  Rv; 


da 


P(l>2  Pi  )P 


I1  d’21  J‘lpi  dp2 


(5.12) 


Ihe  first-order  Doppler  perturbation  is  related  to  the  corresponding  perturbations  in 
bending  angle  as  in 


-Xf,  =  vvfj,  -  v^cv. 


where  J,  and  ft,  are  the  first-order  perturbations  in  ray  bending  angle  in  the  xy-  and 
\ r  plane,  respectively. 

Proceeding  from  this  stage  to  evaluate  <Ulf,  dar>  and  p  via  the  statistical  properties 
p',  and  O;  turns  out  to  be  extremely  difficult,  essentially  for  the  reasons  given  at 
the  end  of  section  4.1 .  As  an  alternative  approach  to  a  direct  theoretical  evaluation,  one 
might  estimate  the  right  hand  side  of  equation  (5.12)  from  occupation  data.  In  so 
doing,  care  should  be  taken  to  properly  subtract  from  the  total  atmospheric  Doppler 
the  lowest-order  component,  f0,  due  to  the  ambient  atmosphere,  as  explained  in 
section  4.4.2.  From  a  record  of  t|  versus  impact  parameter  a,  df,  /da  along  with  its 
mean  square  and  correlation  lunction  p  may  be  estimated,  yielding  even¬ 
tually,  through  the  operations  indicated  in  equation  (5.12),  an  estimate  of  the  mean 
square  fractional  error  in  refractivity. 


We  have  shown  previously  that  in  the  evaluation  of  the  mean  square  fluctuations  in 
••tiding  angle,  the  exact  kernel  in  the  field  integral  (3.1)  cannot  be  approximated  by 
the  quadratic  Fresnel  kernel.  Also  the  complete  spatial  extent  of  the  turbulent 
medium  will  have  to  be  considered  in  such  a  calculation,  as  explained  in  section  4.1. 
Instead  of  attempting  a  rigorous  calculation,  we  shall  make  the  following  heuristic 
arguments,  which  will  hopefully  yield  the  proper  functional  form  and  order  of  magni¬ 
tude  of  the  right  hand  side  of  equation  (5.12).  It  may  be  argued  on  general  grounds 
that  <oq><<d'i>  when  ,<  1.  We  thus  approximate  -f,  -  (Vy/X)/?,,  while  also 
assuming  vv=0(vz).  Due  to  the  assumed  exponential  variation  of  the  rms  turbulent 
refractivity  fluctuations,  the  effective  width  of  the  turbulent  medium  parallel  to  the  limb 
is  iRHjI12.  Perpendicular  to  the  limb  the  effective  width  may  be  shown  to  be  at 
most  on  the  order  of  II  when  0  ''  1.  Ihe  smallest  eddy  sizes  capable  of  scattering 
elementary  waves  within  an  aperture  of  these  dimensions  is  therefore  -  XD/t  and 
kz  "  XD/H  in  the  y-  and  /-directions,  respectively.  The  correlation  length  of  d/3, /da 
can  therefore  be  no  smaller  than  the  least  of  these  dimensions,  i  e,  (Jv.  Also,  when 
'lie  turbulent  phase  perturbation  of  the  wave  is  one  radian,  the  rms  turbulent  bending 
angle  can  be  at  most  on  the  order  of  /3,  -  X/\v.  Omitting  numerical  factors,  we  thus 
lint!  that  an  approximate  upper  bound  on  the  right  hand  side  of  equation  (5.12)  is  such 
that 


( J’U)  ( <L) 

vo  a  |o  [)  R 


(5.13) 


\doptmg  numerical  values  appropriate  for  the  Mariner  10  Venus  and  Pioneer  10 
Jupiter  ex  <  ultations.  we  find  that  the  right  hand  side  of  relation  (5.13)  is  approxi¬ 
mately  O-IO’  and  310'.  respectively.  I  bus  the  refractivity  error  is  small  even  when 
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the  turbulent  phase  perturbation  approaches  the  strong  scattering  regime.  We  also 
observe  from  relation  (5.13)  that  the  refractivity  error  decreases  with  increasing 
planet— spacecraft  distance,  all  other  quantities  remaining  constant.  The  refractivity 
error  incurred  from  the  inversion  of  radio  intensity  measurements  yielded  an  increase  in 
this  error  (proportional  to  D70  for  <t>  <  1;  cf  equations  (4.23)  and  (5.7))  with  pro¬ 
pagation  distance,  and  is  probably  much  larger  than  that  calculated  here  for  Doppler 
measurements  at  any  occultation  depth. 


from  equations  (2.6)  and  (2.7),  and  consistent  with  the  crudeness  of  the  analysis 
presented  here,  we  finally  conclude  that  also  the  fractional  temperature  and  pressure 
errors  are  small,  being  approximately 


(<Ld> 

(  To2 


M  ; 

aF,0  D 


1 

R 


(5.14) 


where,  as  before,  T*,  and  p*i  are  the  lowest-order  errors  in  temperature  and  pres¬ 
sure  relative  to  the  ambient  values  T0  and  p0. 


5.3  Systematic  profile  errors 

The  possibility  that  atmospheric  profiles  obtained  from  occultation  measurements 
were  inflicted  by  systematic  errors  in  addition  to  the  random  errors  discussed  here, 
was  first  raised  by  Hubbard  and  Jokipii  (1975).  As  a  result  of  additional  contribu¬ 
tions  and  subsequent  critical  responses  by  other  authors,  certain  aspects  of  the  theory 
of  these  higher-order  effects,  e  g,  their  wavelength  dependence,  are  still  under  debate 
(cf  discussion  and  references  quoted  in  section  1.1).  However,  both  parties  apparently 
agree  that  the  systematic  effects  arc  too  small  to  be  of  much  practical  concern,  and  a 
discussion  of  this  topic  here  may  therefore  seem  unwarranted.  From  a  theoretical 
point  of  view  these  effects  command  substantial  interest,  however,  and  they  indeed 
appear  to  have  passed  virtually  unnoticed  in  previous  analyses  of  electromagnetic 
wave  propagation  in  randomly  inhomogeneous  media*.  We  shall  therefore  provide 
here  a  brief  discussion  of  these  effects,  developing  first  a  geometrical  optics  result  for 
the  systematic  phase  error  and  proceed  by  generalizing  to  wave  optics.  From  this 
result,  the  form  and  approximate  magnitude  of  the  systematic  errors  in  intensity  and 
bending  angle  will  be  established  along  with  the  resulting  errors  in  the  calculated 
refractivity  profile. 

5.3.1  The  second-order  phase  bias 

Following  the  treatment  of  Eshleman  and  Haugstad  (1977)  and  Haugstad  (1978b), 
consider  a  thin  phase-changing  screen  representing  homogeneous  and  isotropic  turbu¬ 
lence  superimposed  on  a  homogeneous  background,  see  Figure  5.1. 

A  plane  wave  impinging  on  the  screen  as  indicated  is  refracted  through  small  angles 
a,  and  P,  in  the  xz-  and  xy-plancs  of  propagation,  respectively.  The  total  phase 
change  of  the  wave  by  propagating  through  the  screen  to  the  receiver  at  x  =  S/2  +  D  s  I) 
is  kS,  where  in  geometrical  optics  the  optical  path  S  is  simply 
D 

S  =  /  n(x,y(x),  z(x))  ds  (5.15) 


*  Tatarskii  (1971;  pp  2S6)  calculates  the  second-order,  systematic  phase  change  due  to  turbulence, 
but  he  uses  the  result  to  merely  establish  the  conditions  for  convergence  of  the  method  of 
smooth  perturbations.  Keller  (1962)  has  calculated  the  systematic,  second-order  phase  shift  for 
geometrical  optics,  but  his  result  differs  from  that  obtained  here.  Keller  (1977)  has  later  rorrcc 
ted  for  an  error  in  his  calculation,  which  brings  his  result  into  exact  agreement  with  that  given 
by  equation  19.22)  here. 
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Figure  5.1  Radio  ivave  propagation  through  homogeneous  and  isotropic  turbulence 
superimposed  on  a  homogeneous  background 

The  spacecraft  is  located  on  the  x-axis  a  distance  D  >  £  from  the  exit  face  of  the  screen. 

Here  n  =  1  +  f  =  1  +  p0  +  p ,  is  the  refractive  index,  anti 

ds  =  [l  +  (‘i>->2  +  (‘f*)2r<lx  (5.16) 

dx  dx 

is  a  differential  along  the  raypath  represented  by  the  space  curve  y(x),  /.(x). 

The  zero-order  raypath  is  clearly  y0(x)  =  z0(x)  =  0,  since  f0  is  homogeneous.  Denot¬ 
ing  by  Vi(x),  7-i  (x)  the  first-order  perturbations  in  the  raypath  induced  by  the 
turbulence,  we  find  through  second  order  in  small  tpiantitics  by  expanding  f,  and  ds 
relative  to  the  zero-order  raypath 

S  =  /  (1  +  ‘'o  +  v\  4  v\  y  >  i  f  )(1  4  \  /3,2  +  ,7  ay  )dx  (5.17) 

o  - 

All  quantities  are  here  evaluated  on  the  zero-order  raypath,  i  e,  f,  =  tq  (x,  0,0), 
vi  y  =  »>,  (x,0,0)  etc,  and  subscript  y  or  /.  on  tq  denotes  differentiation  with  respect 

to  that  argument.  In  deriving  this  result  we  have  also  used  the  fact  that  0,  =dyi/dx 
and  a,  =  dz,  /dx  for  small  bending  angles  a,  and  0t . 


The 

tions 


integral  in 
V\ 


equation  (5.17)  should  he  evaluated  subject  to  the  boundary  condi- 
=  =V\,L  =0  at  x-0 

(5. 181 


V|  =  z,  -  0  at  x  =  1) 


Expanding  formally  S  =  S0  +  S,  +-S2,  where  S,  =  ©(f,  ),  S2  =  O(i'j),  we  then  obtain 


50  -  /  (1  f  "o)dx 

0 

D 

51  =  /  dx 

o 


s2  -|/  *q  yy,  dx  +  f  Vy iZz,  dx}  +  2  D(0|2  +  0\2 ) 


(5.19a) 

(5.19b) 

(5.19c) 
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In  the  last  result  we  have  used  the  fact  that  Jot]  dx  ^  l)a{  and/(3*  dxs  D(?j  ,  where  Or , 
and  (3,  are  now  the  final,  first-order  bending  angles  in  the  xz-  and  xy  plane,  respecti¬ 
vely.  The  ensemble  averaged  optical  path  length  is  now,  since  <f,>  =  0 

<S>  =  S0+ <S2>  s  S0+ C,  +  C2  (5.20) 

where  we  have  defined,  in  view  of  equation  (5.19c) 

C,  =il)«a,2>  +  </V>) 

(5.21) 

D  i) 

=  /  <"i  vV|>dx  +  /  <y,  /.,>  dx 

o  o 

Integrating  by  parts  the  last  expression  and  applying  the  boundary  conditions  (5.181 
yields 

B2  =  !)(<«, J>  +  < (3 ,2  > ) 

1'hus  there  is  a  net  reduction  in  average  optical  path  by  an  amount 

<S>  S0  =  -  iD(<«,2>  +  <(3,:>)  (5.22) 

relative  to  a  non-turbulent  atmosphere  of  rcfractivity  v0.  An  alternative  statement  of 
this  result  is  that  the  average  phase  velocity  is  increased  in  the  presence  of  turbulence 
(Haugstad  1976,  F.shlcman  and  Haugstad  1977,  Haugstad  1978b).  l'hc  appropriate 
generalization  of  this  result  to  an  inhomogeneous  background  is  (Haugstad  1978b) 

<S2>  =  -  Jj  l)(</)y<(3,-’>  t  0z<a,-’>)  (5.23) 

where  a,  and  (3,  are  the  bending  angles  produced  by  turbulence  on  a  homogeneous 
background. 

The  above  derivation  of  <S2>  apparently  allows  the  following  interpretation  to  be 
made.  The  total  change  in  phase  path  has  two  contributions,  Ki  and  62.  The  former 
component  represents  the  increase  in  phase  path  associated  with  the  longer  geometri¬ 
cal  path  to  the  receiver.  This  increase  is  more  than  offset,  however,  by  the  P2 -compo¬ 
nent,  which  represents  a  twice  as  large  decrease  in  phase  path  associated  with  the  rays 
"success”  in  traversing  the  refractive  medium  in  such  a  way  as  to  encounter,  on  the 
average,  a  lower  than  average  rcfractivity.  We  observe  from  equation  (5.22)  that  the 
bias  in  phase  path  increases  with  propagation  distance  in  vacuum  between  the  turbu¬ 
lent  medium  and  the  receiver. 

A  wave-optical  calculation  shows  that  the  general  weak  scattering  result  for  the  bias 
in  phase  path  is  (Haugstad  1978c) 


<S2>  =  <X,St>  (5.24) 

where  Xi  and  S,  are  the  first-order,  wave  optical  expressions  for  log-amplitude  and 
phase  path.  Employing  equations  (4.7)  and  (4.8)  and  noting  that  2xi  =  .  we  find 

that 

<S2>  =  7rkC  //  ‘f>(ux)sin  (P)dux  (5.25) 

with  I’  given  as  before.  A  comparison  of  this  result  with  the  geometrical  optics 
expression  (5.23)  shows  that  the  two  results  arc  indeed  identical  in  the  limit  k  °°.  A 
further  comparison  of  the  two  results  also  reveals  that  diffraction  strongly  suppresses 
contributions  to  <S2>  from  turbulent  eddy  sizes  smaller  than  the  atmospheric 
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Fresnel  zone.  We  also  find  that  for  Kolmogorov  turbulence  the  bias  in  phase  path 
scales  with  radiation  wavelength  as  A'16,  implying  that  an  initially  non-dispersive 
medium  becomes  slightly  dispersive  by  the  addition  ol  turbulence  (Haugstad  1978c, d, 
Haugstad  and  Eshleman  1979).  A  result  similar  to  equation  (5.25)  has  also  been 
obtained  by  bee  (1976)  from  an  evaluation  of  the  second-order  scattering  produced 
by  simple  diffraction  grating. 

5.3.2  Systematic  errors  in  atmospheric  profiles  from  Doppler  or  intensity  measure- 


If  the  bias  in  optical  path  changes  spatially,  as  a  result  of  a  spatial  variation  in  either 
t>0  or  <t^>,  a  finite  bias  in  both  bending  angle  (and  hence  in  Doppler  frequency)  and 

intensity  may  result,  as  illustrated  in 


ATMOSPHERE 


|v<^>| 


AVERAGE 
PHASE  FRONT 


Figure  5.2  Illustration  of  the  average 
effect  of  turbulence  on  phase 
path,  bending  angle  and  inten¬ 
sity 

The  solid  curve  shows  the  position  of  the  phase 
front  behind  a  homogeneous  ambient  atmosphere 
void  of  turbulence.  The  broken  curve  shows  the 
curved  phase  front  when  turbulence  of  increasing 
strength  in  the  z-direction  is  superimposed. 


intensity  may  result,  as  illustrated  in 
Figure  5.2.  The  turbulence  strength  is 
here  assumed  to  increase  with  atmo¬ 
spheric  depth  so  that  the  average  phase 
advance  is  greater  in  the  denser  region 
of  the  atmosphere  (cf  equations  (5.22) 
or  (5.25)).  It  is  seen  from  the  figure 
that  a  linear  spatial  variation  of 
<v\>  gives  rise  to  a  bias  in  ray  bend¬ 
ing  angle,  whereas  a  non-vanishing 
second  derivative  of  this  quantity  yields 
a  net  focusing  or  dcfocusing  of  the  rays 
and  hence  a  bias  in  signal  intensity. 

Mathematically  the  bias  in  bending 
angle  is  related  to  that  in  phase  path  as 


<os >  =  ,  <B2>  = 

9y  '  9z 


where  a2  and  )32  are  the  second-order 
components  of  the  ray  bending  angle  in 
the  xz-  and  xy  planes  of  propagation, 
respectively.  The  exact  relation  be¬ 
tween  the  intensity  bias  and  that  in 
phase  path  is  still  a  matter  of  contro¬ 
versy,  even  for  the  geometrical  optics 


case  (Eshleman  and  Haugstad  1978, 
Haugstad  1978a,  Eshleman  and  Haugstad  1979,  Hubbard  1979).  However,  for  shallow 
occultations  several  lines  of  reasoning  indicate  that  the  intensity  bias  must  be  of  the 
form 


<02  > 


A  P  9<Q!2>  , 

0/  /  -T  dx 


(5.27) 


in  both  geometrical  optics  and  wave  optics,  since,  for  reasons  of  energy  conservation, 
a  finite  value  of  <02>  is  conditional  on  a  systematic  spatial  variation  of  the  angular 
bias".  As  the  arguments  over  the  precise  form  of  both  the  angular  bias  and  the 
intensity  bias,  along  with  the  question  of  their  wavelength  dependence,  have  become 
extremely  technical  in  nature,  we  shall  make  no  attempt  here  to  provide  a  balanced 
expose  of  these  issues.  Rather,  since  the  resulting  effect  of  these  biases  on  derived 
profiles  is  agreed  by  both  parties  to  be  very  small,  we  shall  instead  use  simple  order 
of  magnitude  arguments  to  merely  establish  the  approximate  magnitude  of  these 


Relation  (5.27)  also  assumes  that  <]S2>  =  0.  This  will  always  be  true  if  both  V0  and  <fj>  have 
strictly  spherical  symmetry,  which  will  be  assumed  here. 


* 
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Considering  the  nominal  occultation  level  (<t>z  -  1/2)  and  disregarding  numerical  fac¬ 
tors,  relations  (5.26)  and  (5.27)  yield 


«x2>  ~  <S2>/H 
<02  >  ~  <S2>/H2 


(5.28) 


where  02  =  02/0z  is  the  fractional  bias  in  intensity,  and  H  is  a  typical  vertical  scale  of 
variation  of  either  <v]>v2  or  p0.  By  virtue  of  the  assumed  spherical  symmetry  of  these 
quantities,  <02>  =  0.  In  equation  (5.26)  the  z-derivative  at  the  receiver  is  related  to  a 
similar  derivative  in  the  atmosphere  according  to  9/dz(receiver)  =  0z9/9z(atmosphere), 
so  that  relations  (5.28),  being  representative  of  (j>z  —  1/2,  actually  overestimates  the 
biases  in  deeper  atmospheric  regions. 


Refractivity  profiles  obtained  from  radio  Doppler  or  intensity  measurements  will  con¬ 
tain  a  bias,  <r.  2>,  that  is  related  to  the  biases  (5.28)  and  to  other  second-order 
effects.  This  refractivity  bias  may  be  obtained  from  equation  (2.1b)  by  formal  expan¬ 
sion  of  v  in  powers  of  v,  and  collecting  all  terms  of  second  order,  yielding 


0*2  > 


_  I 

Doppler  2 

=  I 

Intensity  2 


,  .  9  ?  d<0-,>  ,-2  . 

<’»'>  +  (R)  /“ irp  p 

<„,>  ,  ,,,,  2/ !  |<S.>  - i(i  +  ^)<?t>]p‘'3aP 


(5.29a) 

(5.30a) 


A  closer  examination  of  equation  (5.29a),  using  the  previously  derived  approximate 
upper  bound  on  <v£i>lv o>  reveals  that  the  refractivity  bias  is  probably  dominated  by 
the  second  term  in  this  case,  ie  by  the  angular  bias.  Bv  contrast,  the  right  hand  side 
of  equation  (5.30a)  is  dominated  by  the  second  term  in  the  integrand,  yielding 
approximately  for  the  two  cases 


Doppler  ~  (  b'  ‘  f  P‘  2  dP 


R  o  dp 


0*2  >  intensity  ( 


i  1/2, 


2RD2 


<0,>p'  2  dp 


(5.29b) 

(5.30b) 


Explicit  evaluation  of  these  relations,  using  the  first  of  relations  (5.28)  along  with  an 
explicit  expression  for  <S2>  obtained  from  equation  (5.25),  yields  the  following 
approximate  results  for  the  fractional,  absolute  bias  in  calculated  refractivity 


j  aF,O^P  <'t'i  (5.30) 

"o  Doppler  Raf*0  L0 

(<|iqt2l>)  ~<0]>  (5.31) 

P0  Intensity 

Although  approximate  only,  these  results  are  significant  in  several  respects.  The 
former  result  shows  that  the  fractional  refractivity  bias  is  always  small  in  the  weak 
scattering  limit.  For  the  Mariner  10  Venus  and  the  Pioner  10  Jupiter  occultations,  for 
example,  we  estimate  typical  values  '10'6  — 10'5  if  <v\>/i> \  =  10‘5  is  a  representative 
turbulence  strength.  We  also  observe  that  the  refractivity  bias  decreases  slightly  with 
propagation  distance  (as  D'1*  for  Kolmogorov  turbulence)  and  scales  with  radio  wave¬ 
length  in  the  same  way  as  the  bias  in  phase  path,  i  e,  as  X',<6  for  Kolmogorov 
turbulence.  By  contrast,  the  fractional  systematic  error  in  refractivity  profiles  derived 
from  intensity  measurements  is  controlled  by  the  entirely  different  quantity  <0  2>, 
which  may  approach  unity  in  the  present  theory.  However,  as  long  as  the  scintilla¬ 
tions  are  weak  also  the  fractional  refractivity  bias  will  be  small,  although  being  much 
larger  than  the  corresponding  error  in  profiles  derived  from  Doppler  measurements. 
Thus  the  fundamentally  different  character  and  magnitude  found  for  the  random 
refractivity  error  for  the  two  cases  also  carry  over  to  the  refractivity  bias,  giving  the 
radio  Doppler  technique  a  further  advantage  over  the  less  accurate  intensity  measure¬ 
ments  as  a  tool  for  remote  probing  of  planetary  atmospheres. 
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6  TURBULENCE  IN  DEEP  RADIO  OCCULTA TIONS 
6.1  General  considerations 

Most  previous  spacecraft  occultations  by  planets  have  involved  fairly  shallow  penetra¬ 
tion  of  the  probing  radio  waves  into  their  atmospheres.  Increased  sensitivity  of  the 
radio  links  of  current  (Voyager  1  and  2)  and  future  spacecraft  will  allow  deeper 
atmospheric  regions  to  be  probed  bv  the  radio  occultation  method,  following 
Eshlcman  ct  al  ( 1979a),  we  shall  reserve  the  term  "deep  occultation"  to  denote  situa¬ 
tions  where  atmospheric  refraction  deflects  the  ray  Ivy  about  a  tenth  of  a  planetary 
radius  or  more,  as  viewed  in  the  plane  of  the  sky.  Under  such  circumstances  the 
probing  wave  may  experience  significant  focusing  due  to  the  curved  planetary  limb. 
As  explained  in  sections  3.2  and  3.3,  in  the  central  regions  behind  the  planet  the 
focusing  may  become  comparable  to  and  eventually  substantially  exceed  the  initial 
defocusing  caused  bv  differential  refraction  in  the  plane  of  propagation.  Correspond¬ 
ing  optical  enhancements  have  been  observed  at  earth  during  the  central  occultation 
of  e-Geminorum  by  Mars  (Elliot  ct  al  1977b). 

For  a  general  refractivitv  field  the  focus,  or  more  properly  the  caustic  surface,  is  the 
locus  of  intersection  of  infinitesimally  close  normals  (rays)  to  the  e<|uirefractivity 
contours.  If  these  contours  are  everywhere  parallel  to  the  planetary  limb,  this  locus  is 
also  the  evolute  of  the  planetary  limb.  This  circumstance  has  given  rise  to  the  term 
"evolute  flash"  (Eshleman  ct  al  1979a).  which  refers  to  the  brief  but  very  strong  signal 
enhancement  associated  with  a  crossing  of  the  evolute.  For  an  oblate  planet  characte¬ 
rized  bv  two  orthogonal  radii  of  curvature,  the  locus  resembles  a  four-cusp  cylinder, 
extending  along  a  line  through  the  planetary  center  from  a  certain  minimum  distance 
behind  the  planet  to  infinity.  For  a  strictly  spherical  refractivitv  field  the  evolute 
degenerates  into  a  single  line. 

The  basic  theory  of  stub  evolute  flashes  has  been  worked  out  in  two  recent  publica¬ 
tions.  These  papers  have  either  emphasized  the  potential  information  about  gravitatio¬ 
nal  moments,  planetary  rotation  and  zonal  winds  that  can  be  obtained  from  observa¬ 
tion  of  the  evolute  flash  (Eshleman  ct  al  1979a).  or  the  possibility  of  using  the  gravita¬ 
tional  focus  of  the  sun  for  eavesdropping  and  communication  over  interstellar  dis¬ 
tances  (Eshleman  1979).  I  hc  theoretical  predictions  of  these  papers  were  based, 
however,  on  the  assumption  of  a  spherical  or  near-spherical  distribution  of  refracti¬ 
vitv.  Any  departure  from  such  a  degree  of  symmetry,  as  caused  for  example  by  a 
more  complex  global  distribution  of  refractivitv  or  by  smaller-scale  turbulence,  will  in 
principle  affect  the  predictions  from  this  theory.  The  extent  to  which  this  occurs  is 
investigated  here  by  evaluating  the  general  results  developed  in  section  4  for  the 
scintillation  index  and  intensity  power  spectrum  (lose  to  the  focal  line,  following 
essentially  the  development  of  Haugstad  (1981).  The  general  results  of  this  analysis 
will  then  be  applied  to  the  actual  conditions  of  the  recent  Voyager  1  Jupiter  focal 
evolute  crossing,  and  also  to  examine  the  effect  of  solar  corona  plasma  irregularities 
on  the  performance  of  the  gravitational  lens  of  the  sun. 

For  the  reasons  given  at  the  end  of  section  4.1.  no  attempt  will  be  made  to  evaluate 
the  fluctuations  in  Doppler  frequence  near  the  focus. 


6.2  Scintillation  index  and  power  spectrum  close  to  the  focus 

The  weak  scattering  scintillation  index  associated  with  the  near-limb  ray  is,  from 
equations  (4.10)  and  14.2')) 


87rfk:  //'!>( 


ul)sinJ(n|!2.|l(asu^)|2J*dux 


(6.1) 
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where  the  upper  form  is  for  a  coherent  (radio)  source,  and  the  lower  form  for  an 
extended  (stellar)  source  of  projected  angular  radius  as  above  the  planetary  limb.  For 
the  fur-limb  ray,  the  quantities  0V  and  0Z  contained  in  the  definition  of  P  remain  of 
the  same  form,  except  that  0V -+'|0V|,  since  0V  is  negative  (aD  >  R)  for  the  ray  over 
the  far  limb. 


Partly  because  of  the  ray-optical  nature  of  0y  and  0Z,  and  in  part  because  of  approxi¬ 
mations  made  in  the  derivation  of  equation  (6.1)  itself,  this  result  is  only  valid 
beyond  a  certain  minimum  distance  from  the  focus.  The  former  constraint  coincides 
with  the  distance  w  s  \I3/R  given  in  relation  (3.21)  as  the  width  of  the  focal  maxi¬ 
mum.  I  he  more  constraining  condition  associated  with  the  neglected  terms  in  the 
stationary  phase  expansion  was  established  in  section  4. l  as  0V  =  (R/rx)  <  (RH) 1/2 / 
ay  (),  where  ri  is  the  distance  from  the  spacecraft  to  the  tocai  line  in  the  plane  of 
the  sky.  Thus  equation  (6. 1 )  should  not  be  used  close  to  or  within  a  distance  rL  = 
Uj.  (,  ( R/l i) 1 2  from  the  focal  line. 

Considering  in  the  sequel  the  coherent  source  form  of  equation  (6.1 ),  the  form  of  the 
scintillation  index  in  the  focal  region  is  found  to  depend  strongly  on  which  of  the 
two  inequalities 


aK.h  5  U 


that  are  satisfied.  When  ay  o  <  ay  ^  <  Lo.  the  appropriate  form  of  equation  (6.1)  is 
given  by  equations  (4.23a,b).  Evaluation  of  these  equations  yields  for  0Z  ^  1 


(6.2) 


where  K,  -  (tt1  *  ( 2tr| '-isccjrp/4l2»’- 1  r2  ((p-1  )/2)|  /  |r((p-3)/2)T(p-l)|.  This  result  is 
valid  for  3  <  p  <  6,  and  for  depths  of  occultation  such  that  0Z/0V  —  Hr^R-r^R  <  1, 
while  simultaneously  ay^R/r^1  2  <  l,0.  Thus  focusing  around  the  curved  limb  causes 
the  scintillation  index  to  increase  as  r^5*  (assuming  p  =  11/3)  as  the  focus  is 
approached. 


When  the  occultation  has  proceeded  to  the  level  where  ay  y,  >  L0,  equation  (6.1)  may 
be  evaluated  subject  to  the  stronger  condition  a(.  h  >  I.0,  to  yield  asymptotically 

<0i>  >  K:  l-9-V  <tq>  (6.3) 

A* 

where  K:  1  1 1  (p/2)/  (p-2)I  (3/2)T  ((p-3)/2)].  According  to  this  result  there  is  no 

explicit  ^-dependence  ol  the  scintillation  index  in  this  regime,  the  only  possible 
change  in  <05  >  being  due  to  the  very  slow  change  in  <T{ >  towards  the  focus.  But 
close  to  the  locus,  where  ay  (l  L0  and  0Z  3:  Il/R,  the  change  in  altitude  of  the  ray 
periapsis  as  rA -*  0  is  a  small  fraction  of  a  scale  height  only,  l  hus  <t'j>,  and  hence 
also  <01  >.  is  very  nearly  constant  in  this  region.  This  situation  is  depicted  schemati- 
i  ally  in  Figure  6.1 . 


1  he  immediate  reason  for  the  plateau  in  Figure  6.1  is  the  assumption  imbedded  in  the 
turbulence  power  spectrum,  equation  (4.1 6),  that  d’  is  itself  constant  for  u  I.q  .  A 
constant  power  spectrum  in  this  region  is  merely  a  convenient  way  of  imposing  a 
finite  variance  on  the  fluctuations  in  refractivity,  however,  since  cannot  increase  as 
u  "  all  the  way  towards  tm  0  if  p>  1  and  <iq  >  is  finite.  In  actual  practice  no 
well-defined  outer  scale  exists,  anti  any  structure  in  ‘1*  in  the  range  u<Lq  will,  of 
necessity ,  also  allect  <05  >  close  to  the  focus. 
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Figure  6.1  Illustration  of  the  qualitative  variation  in  scintillation  index  with  fractional 
occultation  depth,  rJR,  for  a  fixed  planet  to  spacecraft  distance 
R  is  the  planetary  radius,  and  rx  is  the  distance  from  the  spacecraft  to  the  focal  line  in 
the  plane  of  the  sky. 


The  point  source  intensity  scintillation  power  spectrum  is,  from  equation  (4.37) 


2 

W(cj)  =  k2cj  Z 
i  =  l 


xsin2{(^-)  [0ycos 2(0y+  6)  +  0zsin2  (0j  +  6)|j- 
7 - - : - dx  ( 6.4) 

i  ll  +  (^X)2F/2(x2_1)h2 


As  before,  oa  =  2tr  1/2v*/a p  0  >s  the  Fresnel  frequency,  co0  =  v*/L0,  and  5  =  tan'1  (v*/ 
v*),  where  Vy  =  0yVy,  v£  =  <Pzv7  are  the  ray  velocity  components  at  the  ray  periapsis 
corresponding  to  the  velocity  components  Vy  and  vz  of  the  spacecraft  in  directions 
parallel  and  perpendicular  to  the  limb,  respectively.  For  the  current  conditions  it  is 
possible  to  derive  an  approximate  analytical  form  of  equation  (6.4).  To  this  end  wc 
first  observe  that  normally  v*/v*  <  (vz/vv)-(10Hrj/R2 )  <  1  in  deep  occultations,  so  that 
6  a  0.  characteristic  of  a  grazing  occultation.  Defining  co^  =  coc0y'2,  and  noting  from 
the  definition  of  coc  that  co'c  in  deep  occultations.  we  find  that 


W(cj)  a  co  sin2  ( )  •  /  — - -  (6.5) 

«c  1  [1  +  (^)2]l>/2(x2  -  l)1  2 

<x)n 

F.xamination  of  the  integral  shows  that 


\V(co)  a 


’sin2(~~)  for 
“c 

w,-psin2  (-^fr) 
to  2 


CO  <  co0 


for  co  >  co0 


(6.6) 


Note  that,  while  initially  co0/co^  ~  (aj.  0/L0)  1  per  assumption,  in  very  deep  occul¬ 

tations  this  ratio  is  reversed,  since  there  co0/co^.  "  (ap  o/LqMR/Tl)1'2  >  1. 


The  relations  (6.6)  imply  that  in  very  deep  occultations  the  scintillation  power 
spectrum  has  the  form  depicted  qualitatively  in  Figure  6.2,  with  the  bulk  scintillation 
power  confined  within  a  band  coj- ~  co  £  co0.  As  the  focus  is  approached,  both  co^. 
and  co0  increase.  The  bandwidth  also  increases,  however,  since  to'c  <*  (R/rx)*  2  while 
co0  ac  R/r±  and  therefore  increases  faster  than  toc  by  the  factor  (R/rx),/2.  The  increase 
in  co0  towards  the  focus  by  the  factor  R/rx  is  caused  alone  by  the  similar  increase  in 
ray  speed  at  the  periapsis.  The  slower  increase  of  to .  is  due,  however,  to  the  compet¬ 
ing  expansion  of  the  horizontal  Fresnel  scale  by  a  factor  (R/rx)12,  which,  if  acting 
alone  would  reduce  tot  by  a  similar  amount,  and  the  increased  ray  speed  which  is 
proportional  to  R/rx. 
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figure  6.2  Qualitative  sketch  of  the  point  source  intensity  scintillation  power  spect¬ 
rum  m  very  deep  occupations 

The  charac  teristic  frequencies  COt  and  to 0  are  defined  in  the  text. 


In  the  temporal  domain,  the  bandlimitcd  power  spectrum  described  above  corresponds 
to  typical  scintillation  time  scales  less  than,  or  at  most  equal  to,  that  associated  with 
the  low  frequency  boundary  of  the  pass-band,  that  is  .  Thus  cj^.  is  also  a 

characteristic  temporal  frequency  in  very  deep  radio  occultations. 

by  way  of  illustration,  it  is  readily  proved  using  the  previous  relations  that  one  scale 
height  from  the  focus,  where  the  background  intensity  has  risen  to  the  unocculted 
lc  el  (0  =  1)  according  to  equation  (3.1 1),  the  lower  boundary  of  the  scintillation 
pass  band  is  at  about  100  11/.,  if  vv  =  1  kins'1 ,  a j.-  ()  =  1  km  and  R  =  10<I  km  arc  repre¬ 
sentative  values.  Such  rapid  scintillations  may  severely  impair  the  ability  to  track  the 
radio  signal  through  an  evolute  Hash,  in  particular  if  the  scintillations  are  also  strong 
(  <0  1  >  2  1). 


6.3  Applications 

6.3.1  l  he  Voyager  1  Jupiter  encounter 

1  he  first  positive  identification  of  the  refractive  focus  of  a  planetary  atmosphere  is 
due  to  Klliot  et  al  ( 1977b),  who  observed  the  "central  flash"  of  e  Geminorum  as  it  was 
occulted  by  Mars  on  8  April  1976.  By  contrast,  the  first  attempted  radio  transmission 
from  the  focus  occurred  by  the  recent  passage  of  Voyager  1  through  the  focal  evolute 
of  Jupiter.  According  to  Martin  et  al  (1980),  the  radio  signal  displayed  several  peaks 
in  intensity  near  the  expected  position  of  the  evolute  crossing,  but  none  could  be 
identified  as  the  evolute  flash.  The  absence  of  a  detectable  evolute  flash  led  Martin 
and  colleagues  to  establish  a  lower  bound  on  the  microwave  absorption,  most  likely 
caused  by  Nil,,  near  the  ^ 4.3  bar  level  of  the  ray  periapsis  when  crossing  of  the 
evolute  occurred.  H'c  estimate  below  the  weak  scattering  scintillation  index  during 
this  evolute  crossing,  and  argue  that  its  high  value  may  have  contributed  to  the 
unsuccessful  attempt  to  identify  the  flash  in  this  case. 

A  formal  calculation  of  the  scintillation  index  near  focus  requires  an  estimate  of 
'  We  will  assume  that 

(6.7) 
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where  i)0  is  the  ambient  refract ivity  at  the  rav  periapsis  during  the  evolute  crossing. 
File  value  in  relation  (6.7)  corresponds  roughlv  to  turbulence  o)  the  strength  found 
immediately  above  the  earth’s  boundary  layer  (liulton  cl  til  1972).  if  also  the  outer 
scale  and  the  ambient  temperature  at  the  ray  periapsis  are  equated  to  1  km  and 
200  K,  respectively.  The  ambient  refractivity  is  related  to  its  value  at  standard  tempe¬ 
rature  and  pressure  (S  I  P)  by 

_  =  IT  lo 
*'o.STP  Po  1 

where  (p0  ,T0  )  refer  to  the  standard  condition  (SIP).  I ’situ;  f,,  s  I  t  0.000  1 10)  as  (or 
a  pure  H: -atmosphere,  we  find  that  tq,  '  S-ll)'4,  which  yields  the  value  '6  -  10'* 
for  <f;>.  Putting  t  =  (RH)1"  and  11  7  30  km,  we  obtain  for  the  S-band  (X  0.12  m) 
occultation  using  equation  (6.2)  and  the  previous  estimates 

<0{ >  a  6  •  10*  ((».«) 

Ibis  high  value  of  the  weak  scattering  scintillation  index  suggests  that  the  scintilla¬ 
tions  have  saturated  to  unity  well  before  the  focus  is  reached.  If  real,  such  strong 
scintillations  would  significantly  change  the  maximum  intensity  relative  to  the  ideal 
prediction  for  an  oblate  planet  (c  /  equation  If)  of  F.shlcman  ct  til  1070a).  further¬ 
more,  the  matched  filter  technique  employed  bv  Martin  <■!  u!  (15)80)  to  detec  t  the 
flash,  requires  signal  coherence  over  a  certain  minimum  time.  In  view  of  the-  discus¬ 
sion  in  section  (i.2  it  is  clear,  however,  that  tins  provision  cannot  be  met  close  to  the 
focus  where  the  scintillation  time  scale  approaches  zero.  1  hits  scintillations  nuv  have 
impaired  the  mere  detection  ol  the  evolute  Hash,  above  and  bevond  the  fundamental 
difficulty  imposed  bv  the  verv  strong  Nil,  absorption  of  the  S-band  radio  signal 
along  the  atmospheric  portion  of  the  lavpath. 

6.2.2  The  gravitational  lens  of  the  sun 

Tshleman  (1  9751)  gives  the  maximum  toe  used  mteusiiv  lot  this  case  as 


0 


max 


2ff  kr„ 


(6.5>l 


where  r,<  is  the  gravitational  radius  ol  the  sun  u..  2‘tVi  no.  I  his  tesult  is  derived  lor  a 
plane,  coherent  wave  of  wavenumber  k  dellected  thioiigb  an  angle'  a  2re.nl  ( a  is  the 
impac  t  parameter)  bv  the  solar  gravitational  held.  \nv  sol.u  olilateness  is  too  small  to 
significantly  degrade  the  maximum  intensific  ation  given  bv  ecpi.uion  (tc.‘)|.  Such 
degradation  max  be  caused,  however.  1  > \  the  sol.u  corona  plasma,  which  imposes 
additional  bending  on  the  wave  tiutix  from  the  locus  bv  an  amount  depending  on  the 
local  plasma  gradient  and  the  radio  wavelength. 


Inference  from  the  continuum  ol  scale  sizes  found  :n  the  solar  wind  plasma  (e  g, 
Jokipii  15)72),  suggests  that  also  the  solar  corona  is  inhomogeneous  on  a  continuum 
of  scale  sizes  ranging  from  the  Dcbvc  length  to  those  associated  with  the  global 
dimensions  of  the  corona. 


Ihe  approximate  effect  of  large-sc  ale  coronal  plasma  stun  lures  on  -.•>  mav  be 

appreciated  bv  finding  the-  condition  under  which  ihe  additional  phase  shift  imposed 
by  the  plasma  is  less  than  X/2.  I  his  is  the  condition  that  the  1  resnel  zone  lie- 
coherent  and  the  signal  intensitv  being  esscnli.dk  unallected  bv  the  plasma  for 
larger  phase  shifts  the  Fresnel  zone  breaks  up  into  two  or  more  non  c  ounce  led  paits. 
with  a  corresponding  drop  in  maximum  intensitv  below  the-  value  given  bv  cqua 
lion  (6-9l.  We  thus  require 

k:<  Sj  >  •»:  I  (t,  10) 

where  <Sj  >  is  the  mean  square  c  hange  in  phase  path  of  the  radio  wave  imposed  bv 
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the  eoronul  plasma.  An  additional  constraint  to  ensure  the  validity  of  equation  (6.9) 
is  that  <0]>  <1.  In  very  deep  occultations,  where  aj.-  h  S>  L0,  this  condition  is,  from 
equation  (6.15 1 

<$]>  '  K,  <tq>  <1  (6.1 1) 

X‘ 

It  is  important  to  observe  that  as  Ion”  as  a  j.  h  >  L0,  the  two  previous  conditions  are 
essentially  identical  regardless  of  the  magnitude  of  L0.  Indeed,  from  an  explicit  evalu¬ 
ation  of  <S^>  it  follows  that  in  this  case 

<01  >  --  2  k:<Sj>  (6.12) 

It  is  thus  immaterial  which  of  conditions  (6. 1  0)  or  (6.11)  that  arc  imposed  to  ensure 
the  validity  of  equation  (6.9).  A  result  similar  to  equation  (6. 1 2)  has  been  obtained 
by  Salpeter  (1967)  for  a  one-dimensional  turbulence  model.  Equation  (6.12) 
extends  this  result  to  three-dimensional  turbulence  on  an  inhomogeneous  back¬ 
ground  when  the  leading  dimension  of  the  atmospheric  Fresnel  zone  substantially 
exceeds  the  outer  scale  of  turbulence. 

Assuming  I,0  '  R  ( R  is  the  solar  radius),  so  that  the  rms  phase  change  is  comparable 
to  that  imposed  by  the  average  corona,  we  calculate  below  the  coronal  change  in 
phase  path  corresponding  to  the  modified  Baumbach-Allen  plasma  model  proposed  by 
lyler  ct  al  ( 1977),  which  accounts  for  latitudinal  variations  and  the  effect  of  the  solar 
wind  plasma.  Denoting  by  p  and  0  the  heliocentric  distance  (in  units  of  the  solar 
radius)  and  solar  latitude,  respectively,  this  model  asserts  that 

\  I  (  +  1  .-55)1  O' 4  4-  9.--A4  10"  I  (COS20  +  ,4Tsin20)1'2  (6.13) 

l'  p16  ph  p2  b4 


Figure  6.  j  1  dines  of  the  phase  shift  x  induced  h\  the  solar  corona  plasma,  for  various 

fractional  heliocentric  distances  p  and  wavelengths  X 
[  hr  corresponding  fractional  distance  D  I)/Dm;n  to  the  focus  is  also  indicated. 


r>6 


where  N'c  is  the  density  of  electrons  per  inJ .  The  associated  refractivity  is 

v  = -4.48- lO'16  X2  Nc  (6.14) 

with  X  in  units  of  meters.  Using  equation  (6.14)  and  the  previous  equations  we  have 
calculated  in  Figure  6.3  the  ratio  k  =  <S2>,;2/X  for  different  values  of  p  and  X.  The 
very  large  values  attained  by  K  for  even  large  heliocentric  distances  (p  2  10)  and  small 
wavelengths  (X  a  1  mm),  strongly  suggest  that  coronal  plasma  irregularities  both 
reduce  the  coherent  size  of  the  Fresnel  zone  (and  thus  0,nax)  and  cause  severe  signal 
fading  as  the  focus  is  approached. 

According  to  Figure  6.3,  near-optimum  conditions  can  be  restored  either  by  moving 
the  ray  periupsis  to  larger  heliocentric  distances,  or  by  reducing  the  wavelength  below 
about  10‘4  m.  The  former  alternative  is  not  particularly  attractive,  since  the  distance 
to  the  focus  behind  the  sun  increases  as  p2,  ant)  is  already  10()l)mjn  at  p=  10 
(l)min~  350  A L'  is  the  distance  to  the  nearest  focus  corresponding  to  p=  1).  The 
apparent  solution  is  therefore  to  use  very  short  wavelengths,  provided  present  techno¬ 
logical  difficulties  of  sub-millimeter  communication  can  be  overcome. 


7  DIAGNOSTICS  OF  SMALL-SCALE  ATMOSPHERIC  STRUCTURES 
7.1  General  considerations 

turbulence  in  the  tree  atmosphere  owes  its  existence  to.  (i)  thermal  instabilities,  (ii) 
wind  shear,  or  (iii)  the  breaking  of  large-scale,  vertically  propagating  waves,  such  as 
internal  gravity  waves.  These  are  physically  distinct  processes  that  mav  lead  to  differ¬ 
ent  spatial  distributions  and  characteristics  of  the  turbulence.  For  example,  a  convec¬ 
tive  instability  may  be  more  likely  to  occur  near  the  equatorial  plane  on  the  dav  side 
where  the  solar  heating  rate  is  a  maximum.  On  the  other  hand,  wind  shear  and 
the  associated  turbulence  may  be  sharply  confined  in  altitude,  but  may  not  be  other¬ 
wise  restricted  to  particular  planetographic  locations.  A  common  feature  of  well- 
developed  turbulence,  regardless  of  the  characteristics  of  the  source,  is  the  universal 
form  of  the  turbulence  power  spectrum  in  the  inertial  subrange  where  LJ,1  <  u  <  CJ,1. 

1  he  magnitude  ol  the  outer  scale,  being  a  characteristic  dimension  of  the  input 
region,  depends  directly  on  the  specifics  of  the  source,  however. 

In  the  turbulence  model  (4.16),  the  two  quantities  <>J>  and  1.0  are  of  principal 
interest.  One  might  also  want  to  estimate  the  slope  parameter  p  from  the  occultation 
data,  as  has  been  done  by  Woo  ct  til  (1976,  1977).  It  should  be  emphasized,  however, 
that  a  statistically  significant  departure  of  p  from  the  universal  value  p=  11/3  in  the 
inertial  subrange,  would  imply  the  existence  of  turbulence  that  is  not  well-developed. 
1  hat  is,  it  is  of  an  intermittent  and  nonstationary  nature.  The  positive  identification 
of  such  intermittent  sources  of  turbulence  would  in  itself  be  important,  but  the 
particular  spectral  shape  recorded  would  be  directly  linked  to  the  specific  event  and 
would  not  bear  any  universal  significance.  The  quantities  <v\>  and  L„  in  well- 
developed  turbulence  would  be  of  much  greater  interest,  since  the  atmospheric  pro¬ 
cesses  generating  the  turbulence  will  be  of  a  more  stationary  character  in  this  case. 
I  he  specifics  of  such  processes,  such  as  the  turbulence  strength  and  characteristic 
dimensions  of  the  source,  may  provide  important  observational  constraints  on  theore¬ 
tical  circ  ulation  models  of  planetary  atmospheres,  as  suggested  initially  by  W'oo  ct  al 
(1974). 

In  order  to  estimate  O',2''  and  L0  from  occultation  data,  one  should  in  principle 
consider  quantities  that  exhibit  the  strongest  possible  dependence  on  the  magnitude  of 
these  parameters.  For  the  former,  the  scintillation  index  and  the  associated  power  spec- 
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trum  is  such  a  quantity,  while  for  lhc  outer  stale  the  phase  stimulations  is  the 
natural  quantity  to  examine.  Also  other  signal  characteristics  may  be  used,  but  will 
not  be  considered  here.* 

We  discuss  below  specifically  how  the  characteristics  of  small-scale  atmospheric  struc¬ 
tures  can  be  obtained  from  occultation  measurements.  It  is  deemed  important  to 
explore  both  the  limitations  and  capabilities  of  the  occultation  technique  in  this 
regard,  since  once  this  potential  is  realized,  such  experiments  may  be  planned  to  take 
1  idlest  possible  advantage  of  their  intrinsic  capabilities.  Parts  of  the  discussion  here 
c  losely  parallel  that  of  Ilaugstad  (1979a). 


7.2  Estimation  of  the  turbulence  strength 

Before  any  "wiggly"  occultation  curve  is  analyzed  lor  its  expected  turbulence  con¬ 
tent.  it  should  be  ascertained  that  the  recorded  small-scale  structures  are  in  fact  due 
to  turbulence  and  not  to  other  atmospheric  struc  tures  (e  g.  "lavers”),  or  to  numerical 
artifacts  introduced  during  the  data  recording  or  the  subsequent  analysis.  The  best 
way  to  make  such  a  check,  would  probably  be  by  examining  the  intensity  power 
spectrum  and  the  scintillation  index.  Assuming  that  the  intensity  background  has 
been  properly  removed,  the  weak  scintillation  power  spectrum  should  conform,  with¬ 
in  proper  limits  ol  statistical  significance,  with  one  of  the  characteristic  shapes  dis¬ 
played  in  Figures  4.5  to  4.9,  for  various  geometries  and  depths  of  occultation.  As  an 
alternative,  one  might  identify  the  Fresnel  "knee”  at  u>c  in  the  intensity  power  spectrum, 
and  determine  if  it  scales  with  wavelength  as  X1  2  for  radio  occultations.  Within  wide 
bounds  on  the-  spectral  slope  parameter  p  (U  <  p<  (i).  this  scaling  is  independent  of 
[>■  In  stellar  occultation  power  spectra,  the  position  and  shape  of  the  corresponding 
"knee"  depend  on  the  relative  size  of  the  projected  stellar  radius  and  the  frec-spaee 
Fresnel  zone.  It  as  a,  ().  the  "knee"  is  determined  by  ns  and  is  therefore  essentially 
independent  of  wavelength.  For  a^  ■  a,  ().  stellar  scintillation  power  spectra  would  in 
all  respects  resemble  closeh  those  of  a  coherent  source. 

A  final  check  on  the  expected  authenticity  of  the  turbulence  is  provided  bv  the 
scintillation  index,  which  scales  with  wavelength  as  A'h  for  values  small  compared  to 
unit  y . 

Once  well-developed  turbulence  lias  been  reliably  identified,  its  strength  and  spatial 
variation  m.t\  in  principle  be  obtained  from  the  observations  via  the  theoretical  rela¬ 
tionships  developed  in  section  4.5.  Thus,  for  shallow  radio  occultations  and  <  1 . 
equation  |  l.L’5al  yields 


<i’f>  j  i.,p  h,  ‘ 2  k ^  ir",x,‘<0IJ>  (7.n 

where  j  is  a  numerical  constant.  If  an  estimate  of  l.«  has  been  obtained,  as  discussed 
in  the  next  section,  equation  ( 7. 1 )  provides  an  estimate  of  <tq2>  in  terms  of  the 
sc  intillation  index  if  also  Ilt  is  known.  A  natural  first  choice  of  Ut,  consistent  with 
the  atmospheric  model  adopted  here,  would  be  11,^11.  A  more  realistic  estimate  can 
be  obtained  from  the  spatial  variation  of  <$2>  itself  if  this  quantity  has  been  com¬ 
puted  over  a  few  turbulence  scale  heights  or  more.  Indeed,  if  good  quality  scintilla¬ 
tion  data  are  available  over  a  range  in  altitude,  the  actual  variation  of  the  scintillation 
index  may  be  used  to  improve  the  exponential  turbulence  model  suggested  here,  and 
thereby  yield  a  more  accurate  estimate  also  of  <»;>■ 


W'iUtrt  al  ( 1977)  have  suggested  using  the  mutual  i  nherence  function  anti  also  the  quantity 
(Woo  I *17. ’•III.  which  is  proportional  to  <S2>-  lhc  latter  quantity  does  not  appear  to  offer  any 
advantage  over  the  s<  mlillation  index.  In  the  weak  Mattering  regime,  incorporation  of  the  inhomo* 
ig  neous  background,  not  considered  bv  Woo  rt  al  (1977).  shows  that  the  mutual  coherence 
function  is  essentially  identical  to  the  phase  scintillations  in  terms  of  its  potential  usefulness  for 
determining  or  b(). 
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For  stellar  occultations  the  analogous  relationship  to  equation  ( 7. 1 )  may  he  obtained 
from  equation  (4. 2.r>).  As  previously  noted,  however,  the  explicit  form  of  this  relation 
ship  depends  critically  on  the  ratio  as/.q  No  general  relationship  ol  the  lornt  (7.1) 
can  therefore  he  given  in  this  case,  except  when  the  ratio  is  either  very  small  or  very 
large  compared  to  unity.  In  the  former  case  equation  (7. 1 )  clearly  applies,  while  in 
the  latter  case  the  appropriate  relationship  is  obtained  from  equation  (4.2(>a)  as 


<*T> 


V  l-o ' 


H,-’ 


R1 2 I)2  <0,2> 


(7.2) 


where  r;  is  a  numerical  constant.  Again  knowledge  of  I.„  and  II,  yields  <c2>  in  terms 
of  the  scintillation  index  for  the  high  altitude  regions  probed  in  stellar  occultation 
experiments. 

Instead  of  <»>2>  one  might  want  to  estimate  the  refractive  index  structure  constant, 
C,2.  Since 


Dn(Ar)  =•  C,2  |Arr  '  ,  H0«Ar«I.„ 


where  I)n  is  the  refractive  index  structure  function,  it  follows  from  the  definitions  of 
I)„  and  I,0  that 
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From  the  previous  results  it  therefore  follows  that 
'krII;'  2R-' 
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(7.4) 


The  upper  form  with  sr  —  2f  is  for  radio  occultations,  and  the  lower  form  with 
ks?z'2ti  for  stellar  occultations.  The  unknown  value  of  l.„  is  here  absorbed  into  tin- 
structure  constant.  1  he  temperature  structure  constant  corresponding  to  (!*  lor  as¬ 
sumed  adiabatic  conditions  is  given  bv  equation  (4.211.  For  radio  occultations  we  find 
that 

/x  , 

C  j  "  2?<7  1 ) 2  1 ,5 1 1, 1  2  R  1  2 1 )  ' l\  ’ "  (7.3) 

1  i'- 

it 

and  with  a  corresponding  result  for  stellar  occultations  when  the  conditions  leading  to 
equation  (7.2)  are  satisfied. 

Thus  equations  (7.1)  and  (7.2)  vield  the  variance  of  the  refractive  index  fluctuations 
in  terms  of  the  scintillation  index  once  l.„  and  IIt  have  been  stipulated.  From  the 
relationships  (7.4)  and  1 7.3).  rough  estimates  of  the  appropriate  structure  constants 
may  also  be  obtained  from  a  measurement  ol  the  scintillation  index,  using  prior 
knowledge  of  11,  only. 

From  the  analysis  of  section  (i  it  may  be  shown  that  formal  relationships  closely 
similar  to  those  developed  here  also  hold  for  deep  occultations. 
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7.3  Estimation  of  the  outer  scale 


As  we  have  seen,  an  approximate  measure  of  L0  is  a  prerequisite  for  estimating  <r,2> 
from  the  scintillation  data.  However,  by  virtue  of  its  character  as  a  typical  dimension 
of  the  macroscopic  flow  generating  the  turbulence,  knowledge  of  the  magnitude  of 
l.„  is  also  important  in  its  own  right. 


The  natural  quantity  to  examine  in  this  case  is  the  phase  scintillation  power  spec¬ 
trum,  as  given  in  analytic  form  and  computed  for  different  occultation  geometries 
and  depths  in  section  4.4.  The  relevant  feature  of  the  spectra  is  the  "knee”  at 
co  s  co0  -  v^/L.o.  Since  v*  is  known,  the  position  of  the  ’’knee”  is  directly  related  to 
the  magnitude  of  1,0.  It  is  moreover  significant  that  the  position  of  the  ’’knee”  is 
independent  of  both  occultation  depth  and  geometry  for  shallow  occupations.  In 
actual  practice  the  turbulence  power  spectrum  will  not  have  the  idealized  form  as¬ 
sumed  here,  with  constant  power  for  u  I/J.  A  significant  departure ,  one  way  or  the 
other,  from  the  average  high  frequency  slope  of  the  phase  spectrum  is  probably  the 
feature  that  in  practice  signifies  the  wavenumber  location  of  the  input  region. 


In  deep  occupations  the  estimation  of  I.0  is  less  easy.  I  he  situation  is  simplest  for 
moderately  deep  occupations,  where  ay  q  <  a yj,  <  L0.  In  this  case  the  spectral  shape 
corresponds  closely  to  that  of  Figure  4.4,  with  the  transition  to  small  fringes  at  co  =  <oc 
approaching  co0  as  ay  h  increases  towards  L„.  At  still  greater  depths  of  penetration  of 
the  ray,  where  ay  h>  I.0,  the  phase  power  spectrum  assumes  the  asymptotic  form 


/  co 2 , 

cos  (  . 

CO  <  CO, 

■> 

-  ft  3  ->  CJ*  . 

COS2  (  ,2/ 

,  cos>  to,, 

COc 

The  power  spectrum  thus  has  the  form  displayed  qualitatively  in  Figure  7.1.  Identifi¬ 
cation  of  I.0  is  seen  to  depend  on  the  possibility  of  detecting  the  onset  of  the  co'8  3 
decrease  in  scintillation  power  above  >co0.  The  large  fringes  in  the  power  spectrum  in 
this  region  may  hamper  this  identification,  a  problem  likely  to  be  compounded  by  the 
expected  rapid  increase  of  co(,  itself  near  the  focus.  On  the  other  hand,  if  the  transi¬ 
tion  in  spectral  shape  at  co  "  co„  has  been  identified,  determination  of  E0  is  readily 
accomplished  by  the  usual  relationship  co0  = 


Figure  7.1  Qualitative  sketch  of  the  phase  scintillation  power  spectrum  in  very  deep 
occultatio  ns 

I  he  i  h.ir.n  U  ristii  frci|ucn<  ics  C0(  and  C00  arc  defined  in  the  text. 
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7.4  Planetary  rotation  and  atmospheric  winds 

We  discuss  here  the  possibility  of  extracting  from  the  scintillation  measurements 
information  about  planetary  rotation  and  atmospheric  winds.  In  the  earth’s  atmo¬ 
sphere  this  technique  has  been  used  extensively,  having  in  fact  reached  the  level  of 
sophistication  where  a  horizontal  wind  profile  has  been  obtained  (Lee  1974).  Scintil¬ 
lation  measurements  have  also  been  used  to  derive  bulk  flow  velocities  of  the  solar 
corona  and  solar  wind  (<  g,  Woo  1978). 

Lite  theoretical  feasibility  of  such  measurements  rests  in  all  cases  on  the  identification 
of  characteristic  features  in  the  scintillation  power  spectrum  whose  location  or  shape 
depends  on  the  atmospheric  velocity  across  the  line  of  sight.  Considering  the  fixed- 
axes  coordinate  system  of  Figure  2.1,  the  ray  velocity  components  at  periapsis  relative 
to  the  atmospheric  motion  are 
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where  a  superscript  s  has  now  been  added  to  the  spacecraft  velocity  components  to 
distinguish  from  the  corresponding  components  of  atmospheric  motion.  v‘‘.  and  va. 
The  latter  may  be  dominated  by  atmospheric  rotation  or  by  winds,  depending  on 
planetary  rotation  rate  and  occultation  geometry. 

While  in  the  previous  sections  we  have  assumed  both  va  and  v‘*  to  be  negligible,  we 
here  ask  for  effects  that  can  be  attributed  to  the  finite  values  of  these  quantities.  For 
a  spherically  symmetric  atmosphere  the  rotational  part  of  vf  is  zero.  Since  vertical 
wind  velocities  are  normally  very  small  compared  to  their  horizontal  counterparts,  we 
shall  assume  in  the  sequel  that  va  =  0.  In  this  case 
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where  the  relative  ray  speed  vK  and  obliquity  factor  m  determine  the  location  and 
shape,  respectively,  of  the  intensity  power  spectrum,  as  depicted  in  Figures  4.6  4.9 

lor  shallow  occultations.  These  spectra  are  also  valid  in  deep  occultations,  provided 
aK  !■<>•  as  explained  in  section  6.2.  In  very  deep  occultations,  where  Uj.-  h  I,0, 
the  fonn  of  the  intensity  power  spectrum  must  in  general  be  established  from  the 
general  result  (4.87).  If  v'  is  such  that  nnl  ?>  1.  the  power  spectrum  has  the  charac¬ 
teristic  shape  display  ed  in  Figure  6.2,  with  the  low  frequency  boundary  related  to  the 
characteristic  frequency  u/  -  27T1  2v*/a j.  (),  and  hence  to  v*  and  va. 

It  is  clear  from  equations  ( 7.8a)  and  (7.8b)  that  in  order  to  relate  v^'.1  to  the  location 
of  the  spectral  "knee"  at  co  =  (to’  for  deep  occultations),  or  to  the  spectral  shape, 
it  is  necessary  that  iva|  be  at  least  comparable  to  the  component  <t>yvy  due  to  the 
spacecraft  alone.  This'  fact  accentuates  in  the  present  context  the' importance  of 
near-central  (/ e,  deep)  occultations.  Indeed,  if  over  some  range  of  the  occultation 
0V |vy|  ^  Ivyi,  then  both  v*  and  m  depend  crucially  on  va  whenever  0zlv^l<  :va|, 
which  is  likely  to  occur  at  some  stage  of  the  occultation'.  In  this  regime,  both  the 
shape  and  location  of  the  power  spectrum  yield  indirect  information  on  the  magni¬ 
tude  of  va.  In  some  cases  the  magnitude  of  v“  may  be  such  as  to  dominate  v*  and  m 
even  for  shallow  occultations.  For  Jupiter,  for  example,  the  rotational  component 
alone  is  about  18  kms'1  for  a  meridional  occultation  over  the  equator.  For  an  occul¬ 
tation  over  the  polar  regions,  the  rotational  component  is  very  small  and  va  may  here 
be  dominated  by  the  wind  component. 


It  is  apparent  from  the  sketchy  discussion  above  that  information  on  wind  speeds  and 
atmospheric  rotation  rates  are  indeed  within  the  capability  of  the  occultation  tech¬ 
nique.  However,  since  such  measurements  depend  critically  on  details  of  the  occulta¬ 
tion  geometry,  it  is  important  to  fully  examine  these  capabilities  in  the  planning  stage 
of  any  spacecraft  mission. 

7.5  Practical  considerations 

In  the  discussion  of  the  previous  sections  we  have  only  briefly  alluded  to  the  several 
practical  problems  involved  in  the  estimation  of  turbulence  parameters  and  atmo¬ 
spheric  motion  from  occultation  measurements. 

The  major  practical  obstacle  in  the  estimation  of  <nj>  appears  to  be  the  limited 
statistical  significance  of  the  scintillation  index.  This  problem  is  most  acute  for  a 
central  or  near-central  occultation,  since  a  meaningful  evaluation  of  <$}>  must 
employ  scintillation  data  covering  a  range  in  altitude  of  less  than  a  turbulence  scale 
height.  In  a  near-central  occultation.  the  number  of  independent  ’’ripples”  in  the 
intensity-time  record  is  at  most  N()  -  1  It /a k.o ,  since  aj.- ()  is  an  approximate  lower 
bound  on  the  correlation  length  of  the  intensity  fluctuations  in  this  case.  Since  N0  is 
typically  on  the  order  of  10.  it  is  clear  that  the  statistical  significance  of  <0j>,  and 
hence  of  <v\>,  may  indeed  be  very  poor  in  this  case.  The  problem  of  low  statistical 
significance  is  compounded  by  the  principal  difficulty  of  subtracting  the  intensity 
background  from  the  fluctuation  part  of  the  signal,  as  discussed  more  fully  in  sec¬ 
tion  4.4.2.  Because  also  1  I/a j.  ,,  =  (0(10),  the  intensity  trend  function  may  have  signifi¬ 
cant  power  at  spatial  frequencies  ~aj.-'(j,  implying  that  irremovable  contamination  of 
the  scintillation  index  and  power  spectrum  from  the  ambient  background  may  be  the 
normal  situation  in  near-central  occultations.  Despite  this  fundamental  difficulty,  it  is 
important  that  an  efficient  and  common  procedure  be  used  to  separate  the  intensity 
background  from  the  fluctuation  part  in  order  to  minimize  this  problem,  as  pre¬ 
viously  argued. 

1  he  situation  is  entirely  different  for  shallow  occultations.  In  this  case  the  ray  has 
only  a  very  slow  vertical  motion,  but  covers  a  large  horizontal  distance  at  essentially 
constant  altitude.  Indeed,  in  the  central  portion  of  the  occultation  the  ray  moves  a 
horizontal  distance  -A  —  (KH,)1*  before  the  altitude  of  the  ray  periapsis  is  changed 
by  one  turbulence  scale  height.  In  this  case  N,  ~C/aF,o^^«<  so  that  considerable 
spatial  averaging  and  corresponding  high  statistical  significance  of  <<?;>  and  </-’(>  is 
possible. 

I  he  practical  estimation  of  l.„  is  probably  most  hampered  by  the  inherent  lack  of  a 
well-defined  value  of  this  quantity.  Indeed,  the  transition  from  spatial  scales  in  the 
inertial  subrange  to  the  input-range  will  for  all  practical  cases  be  gradual,  and  any 
quantity  or  explicit  result  that  depend  on  the  precise  value  of  L0  should  therefore  be 
handled  with  corresponding  caution.  Apart  from  this  inherent  imprecision  of  L0 
itself,  the  statistical  argument  presented  above  apply  for  similar  reasons  also  to  the 
estimation  of  this  quantity. 

I  he  estimation  of  atmospheric  bulk  flow  and  turbulence  parameters  in  deep  occulta- 
tions  is  affected  by  statistical  limitations  in  the  same  way  as  described  above  for 
shallow  occultations.  Kxlrcinely  dose  to  the  focus,  where  both  cu0  and  ui'c  change 
rapidly  due  to  the^ rapid  change  in  tf>v,  the  effective  length  of  the  scintillation  time 
series,  over  which  0,  may  be  considered  stationary,  depends  critically  on  the  various 
parameters  involved  and  must  be  assessed  for  each  particular  case. 

We  have  described  above  some  important  constraints  of  statistical  nature  that  are 
likely  to  affect  the  practical  estimation  ol  turbulence  parameters  and  atmospheric 
motions  from  occultation  measurements.  Neither  the  potential  nor  the  statistical  limi¬ 
tations  of  this  technique  appear  to  have  been  fully  appreciated  in  the  past.  It  is 


ili'cmctl  ini|>‘>Maii!  tli.it  both  aspects  be  .idilrcssrd  .il  an  earliest  possible  stu^e  in  t lit* 
planning  ol  | > Lim-l ,tr\  il\  l>\  missions,  so  .is  to  realize  in  t'.uli  <  .iso  the  maximum 
potential  ol  the  on  nil. il  ion  irclmiipic. 
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ll  is  i  notable  I.,,  l  lli.it  while  elicits  ol  tuilmleni  c  on  i.ulio  .mil  optical  wave  propa- 
, .  1 1 1 .  >n  m  the  earth's  atmosphere  ami  ionosphere  ha\e  .t  ttr.nl  it  1  i  onsidcraMc  attention 
own  the  | i.ist  .let  .nies.  onl\  \ei\  rcccntk  has  this  issue  heeu  [Hirsueil  s\  sterna! ieallv  in 
oiiiilialiou  im  asm vnients  ol  jilanelai-.  at mos|iheres.  I  his  l.u  k  ol  interest  is  natur.tlh 
explained  loi  i.ulio  m  t  iiltations  h\  the  short  hisioix  ol  this  kmil  ol  measurements, 
ami  also  l>\  the  l.u  t  that  the  litst  anaK  ses  ol  swell  ohscrv  at  ions  (/ e.  the  Mariner 
missions  to  Mars  ami  Venusl  prodm  etl  atmos|iherie  pioliles  that  were  lar^eK  accept- 
al>Ie.  I  lie  low  inletest  is  less  expln  ahle  lot  stell.it  oi  iult.it  ions,  the  milch  longer 
histoix  ol  w  fm  ii  is  marreil  !t\  "spikv"  liiminosil \  proliles  imln.iti\i'  ol  tiirlmlenee  or 
i  it  I  ui  small  si  ale  a  t  mospliei  ii  strut  (tires. 


I  he  ptesenl  anaKsis  attempts  to  iate;on/e  am!  anah /e  the  various  aspeets  ol  the 
overall  prohlem  poseil  l>\  turlmlenee  in  plain  laix  atmospluies  suhjei  teil  to  r.ulio  or 
stellat  mi  (illation  measurements.  I  Inis  the  pmn.ux  task  has  heen  to  extern!  existin'; 

llieoix  '>1  elei  ti  oma  ;net  ii  wave  prop. i  ;at ion  in  a  i  atnlom  medium  to  .mount  lot  the 

min  alio  n  neons  ambient  ainiosphete  upon  ulmh  tin  luihuleme  is  siiperinrposeil. 
U  li  i  If  ne.;li;tMe  in  terrestrial  ptop.t  .atioti  i  xpet  intents  inmkin  ;  shot  t  rav,  paths,  in 
i.ulio  am!  stellat  oi  i  nil  at  ion  expel  intent  s  o\  m  mm  h  longer  pat  Its  sin  It  ,  oil  pi  in.;  to  the 
amiuent  atmosphere  bci  omes  \ei\  impott.utl.  I  he  lheoretie.il  results  derived  here  lot 
the  si  ml  ill.it  uni  index  ape.  with  similar  results  obtained  simultaneousk  llliihlt.ud 
I  n  7  x  "t  .i|  I. net  i  V.  oo  ,  ■  I  "ISIH  !>\  »t  he  i  ant  hois.  In  ilte  lorm  expressed  here 
tile  I  o'  a  |  el  lei  i  o|  tin  .  oiiphn  on  '  lie  s.  mi  illal  ton  index  is  i  cadih  assessed.  He  have 
loitm!  that  in  shallow  ,,  ,  nit  a!  io  •  e  ,  oiiplut  ;  to  the  ittiliietit  I  >ai  k  ;t  on  ml  tedmes  the 
si  an  ilia :  ions  :>v  at  most  i  tetoi  .  ,t  two  loi  i  'point  soitree.  while  lor  an  extended 
si'll  ii  so  in a  e  the  red  it.  l  ton  vanes  het  ween  this  limit  when  tin-  pro  jet  t  ed  stellar  t  ad  ills 
.  .initiated  ;  tiie  t a .! 1 1 1 s  >i  tin  Itee  spate  In  sttel  /one.  to  a  maximum  redtte 
on  e.  a;  e,  it xtnt.it  el\  a  l.u  tor  o!  btiwhen  this  ratio  is  larpe.  I  Itese  lank  moderate 

.  hall  -  oe  .  .vnltasled  b\  vet'  uhstanlial  .  hail  ms  in  the  shape  ol  the  illleusitv  power 

spin  t :  w  hit  h  is  now  i  'ti."i  lumtion  ,.|  !>..th  oi  t  ultation  depth  and  peontettx. 

litis,  i  on.  hisio;i.  as  iiitthe]  d.  i  tile,  I  m  tin  appropriate  sections,  express  a  major 
t  lie.  ui  ..  al  :  esiilt  o  |  i  lie  pt  .  sent  w  ■  w  |. 

I  lie  m  i !  v  sis  also  sh,  ,vv  -  ’  b  it  tin  .  mt  ih  .  I  s,  im  illal  ton  I  u  nil  das  .lei  iv  ed  het  e  should 
t  epl.n  .  'I,.-!!  a  n  il.  ms  :ri  tin  s>  im  l.u  >1  Is  v  lot  appr 'o.u  I)  w  Iteitix  .  r  d  il  I  e  rent  ial  ret  t  at  t  ion 

.hat!  .  s  'he  iv. aa  e  int.'usitv  b\  a  si  null,  ant  luionnl  Mink  sis  ol  i.idto  and  optit.il 
s.  t:  1 1 1 ;!  it  lot  is  ,  1 1  u  ;n  1  .  iilt.iiio:!  oi  met  v  at  tons  should  t  Inn  e  I  ore  alvv.tv  s  proceed  I  rout 

tin  .net. ill  'id  lo'i’.iulis  del  iv  e,!  hep  lint  also  othei.  tetreslti.il  appli.  at  tons  mat  lie 
:  Mind  lot  these  i.  stills  x.it.ilite  to  satellite  i  omnium,  at  loll,  where  the  i.imilltltm  .1 
M’i  pad:  is  ini.ai.pt.,:  bv  tin  .  .nth's  itniosiiheie  .u  ionosphere.  ma\  involve  Midi 
■  ■ .  i ! !  <  i .  :  u  i  a !  telt  ■,  I  *  j .  ,pa  n  ton  the  earth’s  at  m  osp  here,  vv  here  the  tu\ 

nav.i'.s  t,  pons  .  •  1  th.umii  imei'ion  mav  also  \  iel.1  siilist.mti.il  values  ol  \|  nD'II. 
xen  ’h  Midi  tin  ell,,  live  distam,  l>  helween  the  medium  and  the  transmitter  or 
I.'.  <  no',  is  .in  ill  in  this  ,  ase  It  is  possihl,  .  Im. ilk,  that  a  version  ol  the  present 


atiou  liieof,  in. iv  I’m I  apphtation  in  uotistii  lot  elastiel  wave  propagation  in 
r  :  -  i'uio.phi  te.  in  lit,  oi  eatt  oi  ui  the  earth's  mtertoi.  w  henever  siillieieitt 
,  ; tents  due  to  tnomahi  s  ut  «  tthet  letup, n  nine  or  d.  itsitv  are  present. 

.null,  ition  in, I  .valuation  ,d  the  wn  olid  older .  svstemalte  propauat  ion  elicits 
.  n  s  1 1 1 1 1 1 1  •  a  pi  t  n  i  ip  al  tesult.  thoii;lt  ol  m. link  theoietiial  interest.  ('.outran  to 
u  ial  !,  '  ill  ol  lluhhaid  ami  |oktpu  tl'ITA),  who  i  om  lulled  that  the  aver,i"C 


phase  was  retarded  by  turbulent e  relative  l < >  a  mm  turbulent  medium  with  simil.tr 
aver.iizc  relrai  tivitv,  we  have  found  from  a  simple  >;eomcl  ri<  .d  tt|>lits  analv  sis  ill. it  till’ 
avcranc  phase  is  m  l.ul  (iilriiiui'il  in  the  presence  ol  /.cro-nic.in  l  m  hulcin  v  .  I.xtcnsiun 
ol  this  result  lo  wave-optics  has  furthermore  demonstrated  a  ucak  iuvi-!i-ni;tli  dcpcu- 
ikiiii'  ol  the  In. is  in  phase  path  lot  phase  speed),  \  leldnr;  the  important  result  that  all 
in il la il\  non  (lispersne  medium  heroines  sightly  dispersive  h\  tin-  addition  ol  turhu 
leme.  When  a  s\  stemat  u  spatial  variation  in  turbulent  c  stien;llt  hi  in  the  amhiellt 
retrartivitv  is  present,  a  bias  in  ra\  heiulini;  alible  and  intensity  mav  also  result. 

While  these  set  mid-order  propagation  ellerts  were  loiind  to  \  ield  only  msiaitil  ir  ant 
errors  in  oeeullation  experiments,  in  other,  more  controlled  experiments  rec |llirit|o 
extreme  at  a  mat  v  ,  such  elicits  mav  have  to  he  act  minted  lot.  A  potential  candidate 
max  lie  relativistic  time  dcl.tv  experiments,  m  ninth  the  ravpallt  is  intercepted  h\  a 
turbulent  neutral  or  loni/etl  medium. 

It  is  a  principal  t  om  Ittsion  In  mi  the  present  work  that  the  propagation  effects  dis¬ 
cussed  above  produce  onlv  veiv  marginal  tutors  in  proldes  ol  leinperature  anti  pres- 
sure  tlerivetl  limn  ottiiliation  measurements,  ptovitletl  that  the  weak  scintillation 
assumption  lundnmcntu!  to  these  t  ah  ulat unis  is  not  viol.tteii.  In  view  ol  the  Iap;c 
•  •riot's  m  the  Ihonei  I11  ami  I  I  I'mitcr  prolihs.  as  t  aitsed  !>v  the  tnazni  I  it  at  it  ui  effect 
discussed  in  section  tins  <  mn  liision  is  in  no  wav  obvious.  1  he  result  is  however, 

lullv  i  i  insistent  v\  it  !i  i.itho  .  n  t  ult  at  ton  t  xpet  imetits,  wltetf  individual  1 -p  proliles 
bom  the  same  atmosplieie  lioth  exhilui  siihsi.uitial  iiiutua!  uv;rccnieiit .  anti  are  also 
lit  t  i.nilv  totisj'tent  with  other  oiiserv  at  tons  white  tomparisons  tan  lit  made.  Stellar 
ot  t  ult  at  ion  pmt  ties  are  i  v  pit  led  I  >  v  a  mu.  h  In -her  iuteilsii  v  find  nation  level,  pn  ihahlv 
violatin',  the  weak  '■  ml  tll.it  n  >n  assumption,  hut  are  also  more  sttst  epl  ible  to  instru¬ 
mental  firms.  It  is  i  lear  that  errors  ol  this  npc  plaved  a  major  role  in  the  several 
independent  ret  nniin  ;s  ol  the  .)St  n||>ii  ot  t  ultatious  liv  |upiter,  which  tlisplavct!  sill) 
si.miial  dm. i  'ti  eniful  1  see .  Ilunteit  ami  Vcvcrhu  1 ' * 7 < • :  lints  stellar  otiiihation  pro- 
tiles.  i!ioii;h  t  learlv  t  om.iniin  ;  M/ablc  errors,  do  not  net  ess.it  ik  <  iiiiiraditt  the  con- 
<  litsimi  licit  the  part  o|  these  errors  ih.it  t  an  lie  at  I  riliutetl  to  turbulence  is  small 
vvhetievet  the  weak  semt ill.it ion  assumption  is  s.ni'lietl. 

I  lie  n.it  1 1 1  a  I  anti  imporlanl  extension  o  I  the  ihtwn  tlcvclopcd  in  ibis  work  will  be  to 
Hit  or;  ioi  ate  stum;  seatterm;  elicits.  I  Ins  extension  should  also  intitule  an  evaluation, 
iv  tbenretii.il  m  iiumei  it  a!  analvsis.  ot  the  mean  s.piuie  I  hit  l  n.it  ions  in  betulitv;  ar.de. 
essential  to  a  more  ri.uirmis  assessment  ot  luiimlente  etletls  ill  ratlin  Doppler  lllt'.i- 
siiiemeiiis  ilium,  iioili  weak  and  stroii;  st  lUriiir;  t  oiithlions.  1  hese  tasks,  prohahlv 
i  null  ;  o!  i  mime!  it  al  nat  me.  nnlet  d  t  on  si  it  ul  t  I  lit  i ni  nu  ipal  remain  ill  ;  problems  in  the 
theno,  o|  •  iii'uiltuu  e  etletls  m  tadto  and  'tell. it  oiiltltation  expei  intents.  I  he  rapid 
■u.  tease  .  ■:  >  i: 1 1 1 i  1  1 1 1  "i  ic  ii  x  w  1 1  h  1 1  tuns: me r n  tit  pth  in  i  nito  m  «  ult.ilious.  .lion; 

wit:  dll  .  .  ||IU!|.  •••!•.  .  I;ise!  V  fil  "spike  Sit  Ut  I  me"  ot  slfji.u  a  aa  t  1 1 1 1  a  t  lot  |  ll  ;!  1 1  t  Hives. 
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.  ai'Ul.ist.  :  ill-  .1  mi  il’.at  loll  III. lex  involves  .,  i..,u  tuple  inti  ua!  even  HI  the  ■ 'lit 
c.Min  HsC'iiu;  a  ■ .  a  - .  whlh  ul  :  w ' '  ilimellsii  Ills  ev  dliatlaili  <  •  I  an  el  hi  In;  >le  Hilt  ua!  Is 
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strong  scintillation  regime  may  be  beyond  feasibility  .  It  is  possible,  however,  that  a 
more  heuristic  approach,  eg  of  the  sort  suggested  by  Salpeter  (1967)  for  turbulence 
on  a  homogeneous  background,  can  be  used  to  treat  strong  scintillations  with  the 
ambient  background  included.  Together  with  selected  numerical  computations,  first  in 
a  one-dimensional  model,  a  combined  approach  of  lids  oii  may  be  the  only  reali¬ 
zable  path  towards  a  better  understanding  of  the  various  aspects  of  strong  scattering 
effects  in  occultation  experiments*. 

An  extension  of  the  present  theory  to  strong  scattering  should  also  encompass  the 
systematic  effects.  In  strong  scattering  they  are  likely  to  involve  all  powers  of  vx ,  and 
should  therefore  exhibit  entirely  different  characteristics  in  this  regime.  Remarkably, 
the  intensity  bias  may  be  the  easiest  quantity  to  generalize  to  strong  scattering.  This 
contention  rests  of  the  result  of  a  preliminary  analysis  (Haugstad  1979b)  which  indi¬ 
cates  that  the  mutual  coherence  function,  which  reduces  to  the  average  intensity 
when  the  two  field  points  of  its  argument  coalesce,  may  be  extended  in  analytical 
form  to  strong  scattering  and  an  inhomogeneous  background.  No  similar  indications 
have  been  found  for  the  phase  and  bending  angle.  These  latter  quantities  are  more 
easily  obtained  from  numerical  calculations,  since  they  both  involve  only  double 
integrals  for  even  a  two-dimensional  propagation  model.  While  a  numerical  evaluation 
may  not  be  satisfactory  for  assessing  the  qualitative  aspects  of  the  systematic  effects 
in  strong  scattering,  this  approach  is  fully  acceptable  for  assessing  their  approximate 
quantitative  impact  on  occultation  experiments. 

The  fundamental  problem  and  question  motivating  the  present  analysis  has  been  to 
what  extent  turbulence  in  planetary  atmospheres  limit  the  accuracy  of  profiles  of 
refractivity,  temperature  and  pressure  obtained  from  occultation  measurements.  Wc 
have  contended,  either  explicitly  from  rigorous  analysis,  or  from  approximate  or 
heuristic  reasoning,  that  in  the  limit  of  weak  scattering  this  influence  is  very  small. 
An  answer  to  the  basic  question  should  be  pursued  also  for  strong  scattering  condi¬ 
tions,  and  we  have  indicated  how  this  might  be  done.  The  fact  that  turbulence 
effects  arc  in  all  probability  small,  at  least  in  weak  scattering,  implies  an  added 
potential  of  the  occultation  technique  to  also  yield  information  on  small-scale  atmo¬ 
spheric  structures.  Thus,  even  though  the  effects  of  turbulence  on  derived  profiles  arc 
small,  its  various  imprints  on  the  occultation  data  may  be  extracted  and  subjected  to 
a  separate  analysis.  Subject  to  the  practical  limitations  discussed  in  section  7.5,  this 
analysis  would  yield  estimates  of  the  turbulence  strength  and  the  outer  scale  in  both 
shallow  and  deep  occultations.  In  very  deep  occultations,  where  atmospheric  rotation 
or  velocities  intrinsic  to  the  atmosphere  dominate  over  the  ray  speed,  information 
also  on  the  two  former  quantities  are  within  the  capabilities  of  the  occultation  tech¬ 
nique.  Since  the  formal,  weak  scattering  scintillation  index  is  typically  very  high  in 
such  experiments,  as  examplified  by  the  estimate  for  the  Voyager  1  Jupiter  evolute 
crossing,  an  extension  to  strong  scattering  conditions  of  the  theory  developed  here  is 
a  prerequisite  for  measurements  of  this  kind.  If  this  has  been  successfully  done,  and  if 
other  sources  of  error  can  be  controlled  and  kept  on  a  low  level,  careful  planning  and 
interpretation  of  occultation  observations  may  provide  a  very  important  source  of 
information  also  on  atmospheric  circulation  and  dynamics. 


*  Recent  results  by  this  author,  based  on  a  slight  generalization  of  the  stationary  phase  technique 
developed  in  section  docs  show,  however,  that  some  analytical,  strong  scattering  results  for  a 
homogeneous  background  (eg,  the  mutual  coherence  function  and  the  asymptotic  scintillation 
index  and  power  spectrum)  may  be  simply  modified  to  accomodate  the  effect  of  an  inhomogeneous 
background. 
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